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RESEARCH NEEDS IN SEDIMENT HYDRAULICS 


Enos J. Carlson,! M. ASCE and Carl R. Miller,2 J.M. ASCE 
(Proc. Paper 953) 


SYNOPSIS 


Considerable progress has been made in the past decade in the field of 
sediment hydraulics, particularly concerning the basic mechanics of the 
phenomena of sediment transport, erosion, and deposition. It is found, how- 
ever, that there are many areas in which the tools available to the engineer 
for solution of sediment problems are inadequate. The Bureau of Reclamation 
in its existing and planned Federal projects is constantly faced with problems 
involving sediment. These problems include sampling techniques and methods 
of computing sediment load, design of stable channels in earth materials, 
channel stability when the natural regime of an existing channel becomes up- 
set, sediment deposits in reservoirs and how they affect storage loss, ag- 
gradation above structures and degradation below structures, and exclusion 
of sediment from diversions into irrigation and powerplant canals. In the 
solution of these problems all the known knowledge and developments of sedi- 
ment transport, scour, and deposition are used. However, in almost every 
phase there is need for research and development, particularly in the applied 
‘sense to increase the reliability and range of application of various formulas, 
and to develop practical procedures that can be applied in the solution of the 
various problems. 

This paper points out some of the problems in which the Bureau of Recla- 
mation in its work has experienced the need for more research in the field of 


sediment hydraulics. 


INTRODUCTION 


Because of the wide scope of the subject matter available, it has been 
necessary to limit the discussion to applying basic research to the solution 
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of certain sediment problems. There are undoubtedly other important sedi- 
ment problems requiring applied research, but which could not be covered in 
this paper. Computation and design methods are only touched upon and de- 
tails on these methods can be obtained from the referenced publications. 

As the title implies, this paper points out areas in which improvement is 
needed and in which inconsistencies exist in the field of sediment hydraulics. 
Actually the science of sediment transport, deposition and design, has ad- 
vanced to a point where it is as accurate as most other hydrologic fields, and 
perhaps more accurate than some. 


Transport and Deposition in Channels 


Suspended Sediment 


In the determination of the suspended sediment load in streams we have 
progressed a long way in recent years with the development of more accurate 
sampling equipment, improved sampling techniques, and more applicable 
procedures for analyzing the data and computing sediment transport. Incor- 
porated with these improvements is the increased basic knowledge of how 
sediment moves in channels.(1)* However, there is still room for improve- 
ment and simplification in all phases of determining the annual suspended 
sediment load of a stream. Those who have made actual sediment measure- 
ments in the field, especially during high flows, and who have been burdened 
with case upon case of pint milk bottles filled with sediment samples, are 
fully aware of the need for technique and equipment improvements. 


Measurement 


Standard samplers. Suspended sediment samplers currently in use by 
most Government agencies include the US D-43, US DH-48, US D-49, and 
US P-46. The first three are depth-integrating samplers only; the DH-48 is 
a wading or hand sampler, and the fourth is a point-integrating sampler that 
can be used for depth integration. The US D-49 is a later and revised model 
of the US D-43. The development of these and other suspended samplers is 
cee in a publication by the Federal Inter-Agency River Basin Com- 
mittee. 

In the determination of the total sediment transport of a stream it is 
usually necessary to combine measured data with empirical computations, 
except when it is possible to sample the total load. Part of the reason for 
this is the physical limitations of the available suspended sediment samplers 
(Figure 1). The distance between the sampling nozzle and the stream bottom 
is given for the various samplers as follows: 


US P-46 0.49 foot 
US D-43 0.43 foot 
US D-49 0.33 foot 


US DH-48 0.30 foot 


The percentage of the total sediment measured, therefore, varies with the 
stream depth and with the sediment concentration curve. Also illustrated in 
Figure 1 is the problem that results when sampling in a stream with moving — 


* Numbers refer to references at the end of the paper. 
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sand dunes. Here it is desirable to stay further away from the bottom than 
the limit of the sampler permits, to avoid obtaining erroneous suspended sedi- 
ment concentrations. Using Imhoff cones in the field, the sand concentrations 

| were determined for the Republican River near Wauneta, Nebraska. It was 

| found that for point samples near the bed, obtained with a DH-48 sampler with 
valve adaptor, concentrations were measured which varied five-fold depend- 

_ ing upon the location of the sampler nozzle with respect to the sand dunes. 


| Photoelectric cell method. Engineers are continuously searching for and 

| devising methods which will do the required sampling and necessary inter- 

| mediate work in a less costly and speedier manner. Use of the sediment 
sampling equipment which has become more or less standardized in the 

| United States requires considerable time and manpower in obtaining the sam- 
ples, analyzing them, and preparing the data in a form that is directly applic- 
able for computing suspended sediment load. As an example, in an intermit- 

tent sampling program it is estimated that to sample a stream of 200 to 300 

' cis at one station and analyze the sediment samples, assembling the data in 

_a form that is applicable to plotting one point on a sediment rating curve 

would take about 1-1/2 man-days, counting travel time. 

| G. Bradeau, Engineer, National Hydraulic Laboratory, Chateau France, 

described a turbidimeter(3) that by using two photoelectric cells with light 

Sources mounted in an instrument having a tunnel through which the sediment- 

| laden water could pass, readings on an indicator would give the concentration 

and size analysis of sediment passing through the tunnel. Although the instru- 

‘ment is in the development stage, Mr. Bradeau indicated it could be utilized 

| to determine sediment sizes and concentrations over a wide range. The use 

of an instrument such as this would cut down very considerably the time and 

: manpower required to obtain and analyze sediment samples. The need for 

research and improvement is evident here. 

The St. Anthony Falls Hydraulic Laboratory in Minneapolis, Minnesota, 
under the sponsorship of the Federal Inter-Agency River Basin Committee, 
is continually improving present samplers and working on the development of 
new ones. Perhaps the samplers have been improved more rapidly than have 
our sampling techniques. 


Equal-transit-rate. Suspended sediment sampling is usually carried out 
by sampling at a given number of verticals across a stream channel or at 
centroids of sections giving equal discharge. A recent improvement in sam- 
pling technique is the equal-transit-rate (E.T.R.) method which has been de- 
veloped by the United States Geological Survey. This method results in a 
more realistic sampling of the sediment movement since it integrates the pro- 
duct of concentration and velocity. It also results in a reduction in laboratory 
costs. Briefly, the procedure is to obtain samples at equidistant verticals 
while lowering and raising the sampler always at the same rate. The rate of 
raising and lowering the sampler is set by the deepest part of the stream 
where it is necessary to move the sampler fast enough to get not more than 
a pint bottle full in the deepest vertical. The samples thus obtained are con- 
solidated for analysis. The E.T.R. procedure may be refined to the point 
where the water discharge could be determined from the sediment sampling 
data. Theoretically, the water obtained in the sampling process is actually a 
measure of discharge if a coefficient for the nozzle can be combined with the 
speed of raising and lowering the sampler to give an accurate integration of 
the velocity distribution. 


a 
+ 
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There are, of course, limitations to the applicability of the E.T.R. proce- 
dure, and these limitations should also be established. It is impractical to set 
down a definite sampling procedure that will fit all field conditions but guide- 
type information is valuable. 


Point sampling. It is sometimes desirable or necessary to obtain point 
samples in measuring suspended sediment. Experience with existing point 
samplers has indicated that improved models are needed. The mechanical 
features of the US P-46,(4) for example, make it a somewhat delicate instru- 
ment and therefore susceptible to operational troubles when field personnel 
are not familiar with the operation. Some improvements have been made to 
the US P-46 sampler recently. 

A point sampling adaptor has been used by the Bureau of Reclamation in 
some instances (Figure 2A). This adaptor can be used on either the US DH-48 
sampler or the heavier depth-integrating samplers. It fits over the sampling 
nozzle and has a spring attachment which opens and closes a flap at the en- 
trance to the nozzle by pulling a line. Where the depths are shallow the 
amount of water entering the exhaust port is so small that it can be neglected. 
For greater depths a flap or plug valve is used on the exhaust port. Although 
some small errors are introduced by the initial inrush conditions and the 
slight interference with the velocity pattern around the nozzle, they are not 
considered serious enough to affect the quantitative answers significantly. 
This adaptor has not received any laboratory analysis to determine its effect 
on accuracy of sampling. It is a field-developed piece of equipment that ful- 
fills the need for a dependable point-sampling device that can be readily 
adapted to a depth-integrating sampler where for some reason the P-46 
sampler cannot be used. 

Another type of point sampler recently developed by the United States 
Geological Survey is operated by an air mechanism. In this sampler a small 
compressed air chamber is attached to the sampler and a valve in the sample 
nozzle is actuated by this pressure. The sampler has been developed pri- 
marily for wading-type samples and is currently receiving field tests that 
thus far indicate good results. 


Sampling near the bottom of a stream. To obtain sediment samples near 
the bed in streams having hard bottoms, a special sampler has been developed 
This sampler was made by the Bureau of Reclamation to obtain Samples near — 
the bottom in a constricted section of Fivemile Creek, Wyoming. It consists 
of a plastic tube 1-3/4 inches inside diameter and 12 inches long with a flap 
gate fitted on each end as shown in Figure 2B. The standard sampler DH-48 
would not allow a sample to be taken close enough to the bed and the nozzles 
were too small for the largest sediment. 


Computation of Suspended load. The length of sediment sampling records 


on most streams in the United States rarely exceeds a few years. It is neces- 
sary in many cases to determine the probable long-time sediment transport _ 
of a stream using the available small amount of data. Suspended sediment 
loads for stations sampled by the United States Geological Survey are pub- 
lished in Quality of Water Supply Papers by River Basins. j 
Practically all reclamation projects require space to be allocated for sedi- 
ment storage in proposed reservoirs. Reservoirs now constructed by the : 
Bureau of Reclamation usually have space for 100 years of sediment accumu- 
lation before any of the water storage allocations are encroached upon. In 
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Providing such space, consideration is given to the cost of the additional in- 

vestment over the period of years before the additional space is needed. Of 

course, many reservoirs on streams with relatively low sediment loads will 
| have no appreciable storage loss for many hundreds of years. 

Several methods are currently in use by various engineers and agencies to 
predict the long-time sediment flow of a stream. The Bureau of Reclamation 
uses the flow-duration, sediment-rating curve method.(5) This procedure, 
as well as most others, is affected by many factors such as the amount of 
sampling data available to define a sediment-rating curve, the seasonal vari- 
ation of the data, and the flow-duration curve for the years of sampling which 
may be wet or dry years. An example of this method, including seasonal vari- 
ation, is shown in Figures 3 to 5 for the San Juan River at Bluff, Utah. 

Generally speaking, most of the sediment transported by a stream, in any 
given period, is carried by discharges that exceed the median, see Figure 5. 
Sampling programs could perhaps be modified to reduce the number of sam- 
ples obtained during low flows and increase the number obtained during higher 
| flows. Sediment studies consistently indicate that there is a narrow range of 
discharges in any stream which, if sampled, would give a measure of the 
average sediment load carried by that stream. These factors all point to the 
: possibility of using selective sampling in arriving at the average sediment 
: load. Examination and analysis of flow-duration and sediment-rating curves 
would lead to a definition of the range in discharges that, if sampled, would 
give the average annual sediment load for various stream types. 


Pump sampling. Often the sediment concentration in a particular stream 
or channel is so small it cannot be accurately determined with the ordinary 
suspended sediment samplers because it is not possible to integrate the flow 
over a sufficient period of time. In most cases this low sediment concentra- 
tion would not be important, but there are times when its determination both 
in quantity and size is important. An example is the diversion of water into 
a powerplant. One method of sampling this type of flow is to use a water 
pump connected to a pipe so arranged as to sample the flows and discharge 
into a container of known volume such as a 50-gallon oil drum. The velocity 
of the water at the intake is adjusted to stream velocity by means of a valve. 
In this way the flow can be sampled over a long period of time, the necessary 
information obtained to determine concentration, and sufficient material ob- 
tained to analyze for size and quantity. This procedure, however, is cumber- 
some, and a simpler procedure is desirable. Undoubtedly others have been 
and will be confronted with this sampling problem. 


Size analysis of suspended sediments. One of the greatest costs in a sedi- 
ment sampling program is the size analysis. Therefore, any sampling pro- 
cedure that will result in fewer samples for laboratory analysis without de- 
creasing data accuracy or increasing field costs would be desirable. Present 
procedures require from three to seven verticals with two samples at each 
vertical to arrive at one river station sediment concentration and size analy- 
sis. If a single sample could be obtained for a given station in a stream ata 
certain discharge that would give a representative answer, only one concen- 
tration and size distribution analysis would be necessary. The E.T.R. sam- 
pling procedure mentioned previously accomplishes this to some degree. 

Reduction in size analysis costs could also be accomplished by improved 
laboratory techniques. The visual accumulation tube (6) method of size analy- 
sis for sands which is becoming more and more widely used in this country 
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is one result of searching for faster and less costly methods. 

Accuracy of determining the size distribution in a given sample is of prim 
importance since this size distribution reflects directly into the methods of 
determining sediment load. Sediment transport determination by such pro- 
cedures as those advocated by Einstein(7) and Schroeder and Raitt(8) depends 
upon the size distribution. The above points are covered in more detail else- 
where but are mentioned here to illustrate the need for accurate size analysi 
determinations. 

Most sediment sampling programs are necessarily restrictive because of 
high sample analysis costs. A single size analysis costs approximately $15, 
and a single concentration analysis costs $2. Any method that can be de- 
veloped to reduce size analysis costs would be advantageous. 


Unit yield method of determining sediment transport. From the sediment 
sampling that has been carried out in various sections of the country and fror 
data obtained from the resurvey of existing reservoirs, the unit yield or sedi- 
ment runoff per square mile of drainage area, with time, is available for 
many streams and in many locations. When it is necessary to make a 
reconnaissance-type estimate of the sediment transport at a particular loca- 
tion on which no data are available, the unit yield from a similar type of 
drainage basin is often used. It would therefore be helpful to all agencies in- 
volved in making sediment estimates if all available data on sediment yields 
by type of drainage basin were readily available. The Federal Inter-Agency 
River Basin Committee has accomplished much along this line by its recent 
publications “Summary of Reservoir Sedimentation Surveys Made in the 
United States, June 19507(9) and “An Inventory on Published and Unpublished 
Sediment Load Data.”(10) It is expected these reports will be kept up-to-date 
Most agencies have carried out sediment studies for many streams by analyz 
ing sampling records, and availability of these data would be helpful. Althoug 
some disagreement is inevitable in certain instances, a compilation of unit 
yield data from these studies would be helpful to all agencies and would reduc 
duplication of effort. In some drainage basins a large percentage of the total 
sediment is contributed by a few tributaries that produce only a small portio: 
of the basin water runoff. To obtain an overall picture of the sediment and 
water yield distribution for the same drainage areas, a basin balance is pre- 
pared which shows water and sediment for the various sub-basins as a per- 
centage of the total basin runoff (Figure 6). 


Bedload 


Involved in the determination of unmeasured or bedload is the definition of 
channel hydraulics and bed material size distribution. In certain streams the 
unmeasured load is insignificant; in others it embraces a large percentage of 
the total sediment in transport. For example, on the Middle Loup River at 
Dunning, Nebraska, more than 50 percent of the total load is unmeasured 
whereas on the Milk River near Nashua, Montana, the unmeasured loadis _ 
estimated at 5 to 10 percent of the total transport. q 


Sampling bed material. An important part of the field data needed in de- 
termining the bed material transport and total transport of a particular 
stream is the bed material size distribution. The importance of obtaining 
representative samples of this type should not be underestimated since most 
formulas for computing bed sediment transport require a determination of 
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jbed material size distribution. An error in the mean sediment size can re- 
sult in an erroneous answer on the unmeasured load. For example, in the 
Schoklitsch bedload formula(11) a change in mean bed size (50 percent finer) 
}from 1.1 mm to 0.8 mm results in a 100 percent increase in the bedload 
{transport, all other factors remaining the same. 

There are many types of equipment available for obtaining underwater bed 
material samples in normal alluvial stream channels. However, there are 
certain phases of the procedure that are not yet adequately defined, for 
example: 


a) How many samples are necessary to define the mean? 

b) Is there a change in the mean bed size with discharge? 

c) How is a representative sample obtained in a cobble-bedded stream? 
d) How are representative samples obtained under deep flowing water ? 
e) How deep should a sample be taken to give particle size information 
\applicable to the various bedload formulas. 

| 


Typical bed material size analysis curves for the Middle Rio Grande River 
jat Section F in the Bernalillo Study reach are shown in Figure 7. Consider- 
able variation in size across the river bed at this section can be noted. The 
studies being carried out by the Bureau of Reclamation in cooperation with 
the United States Geological Survey utilizing the data obtained in this study 
reach include the application of the Einstein bedload function and the Modified 
Einstein Procedure to an alluvial-type stream having a large change in the 
stage-discharge relationship. 

The problem of obtaining a representative sample from the bed of a chan- 
nel containing a preponderance of cobbles or gravel has become more promi- 
nent in recent years in this country as additional developments are planned 
and constructed in mountainous areas. The particular nature of some pro- 
jects, such as diversions into powerplants, makes it necessary to remove the 
material in transport. Various methods are employed at present such as ob- 
taining several sacks of material by use of a shovel, at low flow stream bank 
water surface level, or perhaps by a pickup method such as described by 
Wolman. (12) 


Bedload samplers. Bedload problems have not been as numerous in the 
United States as those of suspended load because most of the important rivers 
have beds composed of fine material and therefore carry much more sediment 
in suspension than along the bed. Because of this, the development and use of 
bedload samplers in this country have not kept up with that of suspended sedi- 
ment samplers. In the 1940 Report No. 2 “Equipment Used for Sampling Bed 
Load and Bed Material,”(13) various bedload samplers are described. Except 
for a few experimental cases the samplers discussed were developed in 
foreign countries and used there. 

It is believed that a bedload sampler for sand bed material could be de- 
veloped which would give a better indication of actual bed movement than is 
indicated by use of bedload formulas. 

One of the difficulties in developing an adequate bedload sampler has been 
preventing scour around the apparatus and having the sampler conform to the 
sand waves of the stream bed. Possibly using one of the new tough pliable 
plastics for a bottom would work. For sands a comparatively small opening 
ipstream with a Venturi-type pressure reducer to extract the water and 
leave the sediment in a reservoir of the sampler may work effectively. The 
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Venturi principle is used on some suspended load samplers to obtain the 
proper intake velocity. This same principle could be used on bedload sam- 
plers to obtain the proper intake velocity. Any bedload sampler could be used 
if calibrated for conditions similar to those in which it will be used. Obtain- 
ing a few bedload samples would give greater confidence in using bedload 
formulas that give comparable answers. 

In streams where the bed material is very coarse, the question frequently 
arises: At what flow does the coarse material begin to move? Georges 
Labaye(14) describes the use of a hydrophonic detector in which the beginning 
of transport of coarse material was indicated by sound against an immersed 
plate and amplified by a microphone pickup. For wide streams this method 
could be used to determine the fraction of the total width over which the bed- 
load was moving. 


Computing bedload. In computing the sediment loads of streams for plan- 
ning reservoirs, an estimate based on empirical formulas is usually made 
for the unmeasured load. Most bedload formulas are based on drag flume 
experiments, and when applied to river problems may not give accurate re- 
sults. Different formulas give widely different bedloads for the same stream, 
as shown in Figure 8. This is probably due to the fact that each formula 
works well for the conditions under which it was developed, but these condi- 
tions were different for different formulas, and when applied to river prob- 
lems they give widely varied answers. Each of the various formulas that 
have been developed requires essentially the same basic data; channel cross 
section, water temperature, water discharge, energy gradient, and size dis- 
tribution of bed material. Some of these factors can be defined fairly accu- 
rately by field measurement, while others are not so easily obtained. Some 
recent procedures for determining the bed material load have been developed 
that do narrow down the deviation in answers obtained. (7) (15) (16) 

It should be the ultimate goal to produce sediment records that evaluate the 
total transport of sediment and its size distribution for a particular river and 
section. Most suspended sampling is currently being carried out by the 
United States Geological Survey, often at their regular stream gaging stations. 
and records are published in Water Supply Papers. If the total sediment load 
was published or some ratio or percentage of measured to unmeasured load, 
it would be much more valuable to those concerned with using sediment re- 
cords. The suspended sediment records are only a partial measure of the 
total transport. 

Since there does seem to be a definite relationship between the coarser 
fractions of suspended sediment load and bedload for varying stream types, 
it might be possible to improve or elaborate on relationships such as de- 
veloped by Lane and Borland(17) by evaluating existing data from streams in 
this country and throughout the world. This method of arriving at bedload 
transport is considered good only for reconnaissance-type estimates. It can, 
however, point out the relative magnitude of the unmeasured load and help 
determine the need for any further more detailed investigation. 


Channel Stability 


Cohesive materials and apparent cohesive materials. The reclamation 
engineer is continually faced with the problem of choosing proper design 


criteria in designing canals in earth materials. The problem is to select the 
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smallest cross section for economy but which is large enough to carry the de- 
sign discharge on a slope that will not cause the earth material to scour. A 
method for design of canals in non-cohesive materials(18) has been developed 
in the Hydraulic Laboratory of the Bureau of Reclamation. It is based on the 
| critical tractive force principle, i.e., a given size of sediment on a stream 
bed requires a certain tractive force T = (wds) to move it. The use of this 
principle to design the shape of the canal (19) in conjunction with the slope is 

| not difficult for noncohesive materials. However, for clays and other plastic 
soils where cohesion is a predominant force (greater than resistance due to 
weight of individual particles) this process cannot be used directly because 
the resistance to erosion due to cohesive forces has not been determined. 
There is possibly a soil mechanics test or group of tests that can be corre- 
lated with resistance to scour or critical tractive force. The plasticity index, 
which is the difference in percent of moisture between plastic limit and liquid 
limit in Atterberg soils tests,(20) has been suggested as a soil characteristic 
that can be used to indicate resistance to scour. For canal design it has been 
_ Mentioned that a plasticity index (P.I.) of 7 is the critical value, with scour 
occurring for moderate tractive forces below this value. Cases frequently 
occur where the P.I. is above 7, but we still get scour. Some fine silty type 
soils such as volcanic ash and glacial flour may have an apparent cohesion 
and even show a P.I. higher than 7. However, if observed under water, their 
soil structure breaks down readily when subjected to a hydraulic tractive 
force. 

This emphasizes the need for laboratory tests and field data correlating 
soils characteristics with critical tractive force and resistance to shear. 
Plasticity index in combination with consolidated shear tests may be charac- 
teristics that can be used. 


Vegetation. In many waste channels and some canals it is desirable to in- 
crease the stability of the soil by planting or encouraging growth of vegetation. 
For slow velocities, grass causes considerable resistance to flow, but for high 
velocities, where grass is forced to lie down, the resistance to flow may even 
be decreased. Certainly on ephemeral drains grass would be desirable, but 
if a complete grass cover cannot be assured over the full length of the waste- 
‘way, scour will result from high velocities in the denuded areas.(21) It may 
be desirable to allow enough water to escape into a drain to assure a vegeta- 
tive growth. Research in the field of channel stability should include an 
evaluation of the influence of vegetation on the resistance to erosive forces. 


Channels carrying suspended sediment. There is a definite effect on chan- 


nel stability where suspended sediment is carried in the stream, but just to 
what degree has not been determined on a scientific basis. Fortier and 
Scobey(22) observed the stability of several canals some of which carried 
clear water and others carried suspended sediment. Their table shows allow- 
able or critical velocities for both conditions. Lane(18) plotted allowable 
tractive forces for clear water, water with low content of fine sediment, and 
for water with high content of fine sediment. The data of these researchers 
indicate a higher allowable tractive force for water-carrying suspended sedi- 
ment. Fortier and Scobey say, for water-carrying colloids, there should be 
a relationship relating characteristics of bed sediment, suspended sediment 
and hydraulic characteristics of the stream that will show just how the allow- 
able tractive force can be raised for a given concentration of suspended sedi- 
ment. The concentration of suspended sediment very likely affects the 
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stability of cohesive soils as well as bedload movement and formation of 
dunes in noncohesive soils. 


Ephemeral and perennial streams. The degree of stability of stream chan 
nels, such as proposed project drains, is dependent to a large extent on wheth 
er the channel is or will be ephemeral or perennial. For the purpose of this 
paper, an ephemeral stream is one in which the flow is discontinuous and is 
directly related to storm runoff in rivers and to irrigation return flow in 
drains and occurs for short periods of time, seldom exceeding a few days. 

A perennial stream is one in which the flow is continuous for long periods of 
time each year, even though the stream may become dry at the end of the 
runoff season. It is known that an ephemeral channel can withstand higher 
tractive forces than a perennial channel. Since in ephemeral channels vege- 
tation can develop and damaging high flows are of short duration, the bed and 
the banks do not become saturated and therefore do not deteriorate under the 
same tractive forces as a perennial channel. Vegetation, in most cases, can- 
not develop under constantly flowing water conditions, and this is one reason 
why perennial streams are not as resistant to forces of flowing water as 
ephemeral streams. It is necessary in the design of a surface drain to deter- 
mine the probable nature of the flow and then design the channel so that the 
occurring tractive forces will not produce excessive channel deterioration. 

The most striking examples of the effect of the flow duration pattern on 
channel stability are those in which flows introduced into a channel are dif- 
ferent from the flow duration pattern for which the channel was developed 
and became stable. The new flow in a channel that has alluvial-type materials 
will cause rapid deterioration as the channel attempts to adjust to the new 
conditions—often resulting in expensive damage and costly repair. (23) 

Research is therefore needed to define adequately the limiting tractive 
forces for various types of channel conditions. This research can be a com- 
bination of laboratory and field data. There are numerous channels on irri- 
gation projects in this country where field data can be obtained for a wide 
variety of conditions. Some work has been done along this line, and additional 
studies are currently underway by the Bureau of Reclamation and others. 


Stabilizing methods. The need for research in the field of channel stability 
is not only one of determining what flow characteristics and sediment load 
result in a stable channel but also what method of control should be applied 
in stabilizing channels already deteriorating and in maintaining the channel 
shape determined by some applicable method. Assuming that the channel 
shape has been adequately defined, the type of control installations to hold 
this shape must be selected. Many types of channel control are used. These 
include jacks (wooden, concrete, and steel), groins (brush or rock-filled), 
riprap, check dams, pervious fencing, and vegetation. The engineer must 
determine the applicable type of protection that will result in the most eco- 
nomical design, minimum initial expenditure, and low future maintenance 
cost. It has been found that a single row of jacks is most effective for flows 
in which they are almost completely inundated and on streams that carry con- 
siderable debris and fine sediment. Wooden jacks or a more permanent type | 
can be used. If wooden jacks are used it must be assumed that their life is 
only between 5 and 10 years. If they deteriorate, the question arises as to 
whether planted and/or natural vegetation will hold the shape created by the 
jacks. Usually vegetation is employed in conjunction with jack installations 
so the job of holding the channel will fall upon the vegetation that has 


ASCE CARLSON - MILLER 953-11 


developed. This being the case, it is necessary to evaluate the ability of the 
| vegetation to resist erosion and the conditions under which it can no longer 
| withstand the erosive forces. 

Groins are found to be most effective in protecting banks at lower flows or 
flows at which they are not inundated. The correct spacing and the angle to 
the flow are all-important. Often groins and jacks are used in combination. 
if rock riprap is used the required thickness, the points in the channel layout 
where it is needed, and the quality of rock required must be established. 


| Riprap below structures. A special case of stable channel design is that 
of channels below hydraulic structures where turbulence and nonuniform 
velocity distribution occur. To take care of the extra forces which are 

| created by drops, chutes, canal headworks, sluiceways, and other similar 

_ Structures, stilling basins supplemented by riprap on earth slopes are pro- 
vided. Choosing size and thickness of riprap and length of canal on which the 
riprap should be placed has been a matter of judgment and experience. Com- 
_ paring scour patterns in coarse sand for model studies below spillways, still- 
| ing basins, canal drops, etc., has been used extensively to establish uniform 
distribution of erosive forces. 


: Canal beach slope. On practically all earth canals built of gravelly mate- 
rial that are operated at a constant depth a beach slope can be noticed at the 
water surface which is flatter than the side slope. At the water surface, the 
Canal banks are subjected to waves caused by wind and other disturbances, 
and these waves exert forces which cause the earth material to wash down 
and take a flatter slope than the constructed side slope. On earth canals 
that are constructed with side slopes of 1-1/2 to 2:1, it is not uncommon to 
observe beaching slopes at the water surface of 5 or 6 tol. 

In conducting some wave studies in the Hydraulic Laboratory (24) to de- 
termine the required cover blanket to prevent fine base material from leach- 
ing due to wave action on Kennewick Main Canal, Yakima Project, Washington, 
several tests were made on the pit-run gravel shipped from the canal site. 
Figure 9 shows a beaching slope of 5.14 to 1 on the gravel material which had 
been subjected to 42 waves per minute with a wave height of 0.43-foot trough 
to crest for a period of 19 hours. 

Canal design which could include a flatter slope in the water surface area 
would approach more closely the natural tendencies in canals constructed in 
gravelly materials. Because of this tendency to take a beach slope near the 
water surface, a design which included a flatter side slope near the water 
surface would be approaching more closely the natural tendency of the canal. 


Gravel blanket over fine-grained noncohesive material. Many canals are 


built in areas where the earth material is fine-grained and noncohesive. The 
material is usually of rather uniform size which means once erosion begins 
it does not stop from the effect of larger particles being exposed. The canal 
is usually stabilized by artificial means. The use of a gravel blanket has 
proved quite satisfactory in this case. Especially where seepage is not a 
problem, the gravel blanket is very useful in providing resistance to higher 
critical tractive forces, resistance to freezing and thawing forces, and the 
‘cost is low compared to concrete and asphalt concrete. 

A typical example of a canal under construction in fine-grained material 
is that of Kennewick Main Canal referred to under canal beach slope above. 
The base material has a median grain size of 0.05 mm with 20 percent 


| 
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coarser than 0.12 mm and 5 percent finer than 0.005 mm. This base material 
has an apparent cohesion when compacted at the best moisture content, but 
when placed under water, the cohesion readily breaks down. Hydraulic Lab- 
oratory Report No. Hyd-381(24) describes wave studies made on the base 
material, a pit-run gravel, and an angular talus material to determine the 
best method of placing a protective blanket to give most resistance to leach- 
ing of the base material. The model studies showed that the angular talus 
material which had been separated on a 3/4-inch screen with the coarse frac- 
tion placed over the fine fraction gave the best protection to the base materi- 
al. 

The problem that is very important is how thick a gravel blanket to use 
and what restrictions must be made for the preparation and placing of the 
gravel blanket. For practical and economical reasons, a pit-run gravel of- 
fers the best solution, but from an operation and maintenance standpoint the 
pit-run gravel available may have a large preponderance of fines which will 
wash out, causing the coarser particles to settle or roll down the slope. The 
result may be a cover so thin it would not be adequate to give protection, and 
also, the canal prism may have a tendency to be filled up. The ideal gravel 
blanket is material that is placed so there are large enough particles on the 
surface to prevent movement from tractive forces and surface waves with 
enough sand between the coarse gravel and original bank material to prevent 
leaching of the bank material. 


| 


Methods of design of stable channels. The research needs in the design of 
stable canal sections in earth materials may be stated in this question, “What 
are the limiting tractive forces that can be successfully resisted by the mate- 
rial of which the canal perimeter is constructed, and what other factors may 
be present in the future to cause this critical tractive force to be higher or 
lower?” The other factors may include vegetative growth on the canal bed 
and banks, canal structures such as drops, chutes, bridges, and checks, and 
the number and degree of bends. These factors would tend to increase or de- 
crease the forces of resistance to the forces which cause scour. 

Designing a channel that will be stable in an alluvial-type material, or 
stabilizing a channel that is already in a state of deterioration, poses many 
problems. There are several methods of approach that can be utilized, but 
the selection and application of any method is dependent upon the engineers’ 
judgment and the existing conditions. The procedures advocated by various 
engineers can serve only as guides and none of them are completely applic- 
able to every situation. P. W. Terrell and W. M. Borland(25) have presented 
some of the requirements and Bureau experience in the design of earth canals 
and other channels. Some of the more familiar publications on stable channel 
design are those by Lane, (18) Lacey,(26) Blench,(27) and Maddock and 
Leopold. (28 

The Lacey and Blench equations for stable channel design are developed 
from conditions in canals and channels in India where conditions of flow, bed 
and bank materials, and channel hydraulics differ from those encountered in 
the western United States. Their equations also require considerable judg- 
ment or basic assumptions on the part of the designing engineer. This may 
be difficult where insufficient data or knowledge of the particular channel are 
available to influence adequately the engineers’ judgment or assumptions 
within a reasonable degree of accuracy. Lane’s design criteria consider the 
tractive forces acting on the periphery of a canal. The Maddock- Leopold 
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relationships require knowledge of certain factors such as sediment transport 
and settling velocity of the sediments in transport that may not be available 

| or are difficult to determine. These relationships are based on analyses of 
many streams in the western United States and are overall concepts that may 
_ hot always be applicable to a particular reach of a given stream. The 
Maddock-Leopold equation currently in use by the Bureau of Reclamation is 


| as follows: 
gl /2 3 
) wa n 0-555 
_ . 3 
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where 


Q = discharge in cfs 

Vs = weighted mean settling velocity for total sediment load 
S = stream gradient 

Cg = total load sediment concentration 

n = Manning roughness coefficient 

W =channel width 

d =channél mean depth 


In this equation the sediment concentration is required and is therefore un- 
satisfactory for clear water conditions such as may occur below a dam or 
where the sediment load is eliminated by channel rectification. 

Research is currently being carried out by the Bureau of Reclamation us- 
ing data obtained in the field from existing canals, drains, and other natural 
and artificial channels, to evaluate more clearly the requirements for stability. 
Of particular concern are the requirements for stability in materials within 
the coarse silt and fine sand bracket where there is a transition from cohesive 
to non-cohesive conditions. 

In the design of a stable channel, one of the first factors that must be de- 
termined is the so-called “dominant discharge.” This might be defined as 

‘that discharge of a natural channel which determines the characteristics and 
principal dimensions of the channel or is most influential in dictating the 
general condition of stability. In some cases it is easily determined, where- 
as in others it is not. Many rules have been presented to arrive at this dis- 
charge, none of which necessarily give the same answer. We find that the 
dominant discharge varies with the sediment characteristics, with the rela- 
tionship between maximum and mean discharge, with the shape of the flow 
duration curve, and with the flood frequency. The establishment of the domi- 
nant discharge is still one of the most difficult problems in channel design 
work. 

Another factor that is not always easily evaluated is the size distribution 
of the bed materials under design conditions. Will it be the same as for the 
channel before rectification, or will it have an entirely different size distri- 
bution curve? 

Still another factor that enters into the design is the probable sediment 
transport under stabilized conditions. For example, when the sediment in 
transport is being derived from the bed and banks what will be the transport 
conditions and the channel requirements for such transport after rectification 
reduces or eliminates the source of sediment as compared with the design 
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based on existing conditions of sediment transport? 

The above discussion points out only some of the elusive factors in chan- 
nel design work. Engineers working on channel design problems do the best 
they can with the tools available and attain a surprising degree of success 
considering the inadequacy of these tools, but the need for more exacting 
methods of design is evident. 


Deposition Behind Structures 


Distribution of sediment in reservoirs. One of the more common problems 
with which the sedimentologist is faced is the prediction of how the sediment 
will be deposited within a proposed reservoir. This prediction is necessary 
to aid the designers in the location and type of outlet works and to determine 
the probable loss of storage capacity within the various storage allocations 
over a given period of time. Most of the present procedures used to predict 
sediment distribution are based on information obtained from resurvey of 
other reservoirs. The old assumption that all the sediment will deposit in 
the bottom of the reservoir has been refuted. 

The sediment distribution in reservoirs is primarily dependent upon the 
following: 


a) Type of reservoir operation 

b) Reservoir configuration 

c) Size distribution of the inflowing sediments 
d) Flood detention period 


There are, of course, other factors that influence the distribution to varying 
degrees such as the trap efficiently, narrow necks within the reservoir area, 
vegetative growth in the delta area, heavy sediment contributing tributaries 
within the reservoir area, and density currents. 

The publications available on sediment distribution do not completely 
evaluate the above-listed major items that influence distribution. Herein 
lies the need for research that can result in better procedures for predicting 
how the sediment will be deposited. A publication by the Bureau of Reclama- 
tion which is in the form of a progress report(29) has attempted to evaluate 
some of the influencing factors. Evaluation of other factors is still depend- 
ent to a large degree upon the engineers’ judgment. Work is being continued 
to improve the procedures presented in these reports. Of course, the predic- 
tion of sediment distribution is not an important item in all reservoirs. In 
those where it is important, but not a major concern, rapid methods of ap- 
proximating the distribution can be utilized. 

We believe that there is sufficient data available on numerous existing 
reservoirs to correlate the influencing factors and arrive at more accurate 
methods of predicting sediment distribution. Assembling and evaluating these 
data would be a noteworthy contribution to the field of sediment hydraulics. 


Deposition above spillway level. In computing a backwater curve above an_ 
obstruction in a river channel, one must assume a discharge, river cross 
sections, and value of channel resistance which can be used in conjunction 
with some stream flow formula. If the stream is a heavy sediment carrier 
the backwater curve will change with time as sediment deposition takes place. 
The rate of sediment deposition is dependent on the reduction in velocity. 
Locally this can be caused by a wide variety of obstructions that occur or are 
built into the stream channel, see Figure 10. 
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Density currents. Much literature has been devoted to the subject of densi- 
ty currents as indicated in the Bibliography of the Bureau of Reclamation Hy- 
draulic Laboratory Report No. Hyd-373.(30) Experience has shown, in the 
United States, that very little sediment has passed through reservoirs as den- 
sity currents, particularly as the reservoir becomes older. The reason for 
this seems to be the fact that the surface width of the sediment deposit be- 
comes wider and wider as more sediment deposits in the reservoir. The 
turbidity current decreases its velocity and less sediment reaches the dam. 
Also, in cases where density currents do occur it may not be possible to re- 
lease water to take advantage of the density currents. 


Effect on intakes at diversions. Sediment deposition behind diversion 


| structures on alluvial streams occurs rapidly as soon as the diversion dam 


is built. In at least two cases on the Republican River—Superior-Courtland 


_ and Cambridge Diversion Dams—the first floods occurring after the dams 


were completed practically filled their pools with sediment. It is common 
practice in the Bureau of Reclamation to include a gated sluiceway providing 
a low-level opening through the dam near the canal intake. One of the vital 
purposes of this sluiceway is to maintain an open channel leading to the canal 


_ headworks during the recession of a flood. Instructions are issued to opera- 


tors in the field suggesting the sluice gates be opened during a receding flood 
in order to insure an open channel to the headworks. For low diversion dams, 


sediment deposited uring a flood could very easily block off a canal head- 


works, causing the normal river flow to go over the spillway which may be a 
few feet higher than the headworks crest. After diversions are constructed 
the best operation criteria are developed by testing and utilizing the operation 
criteria supplied by the designers. 

Where the dam is of low-head type and the sluiceway crest is only a few 
feet below the headworks crest the sediment concentration in the sluiceway 
and headworks is practically the same for continuous sluicing. Model studies 
on the Bartley Diversion Dam(31) showed intermittent sluicing increased the 
total sediment load through the sluiceway many times over that for continuous 


sluicing. 

Channelization through depositional areas. One of the most perplexing 
problems that is associated with sediment deposition behind structures, par- 
ticularly in the western United States, is delta development. The channel in 
the backwater area of a reservoir becomes plugged with a combination of 
sediment deposits and vegetation which is generally in the form of water- 
loving plants such as phreatophytes and tules. This in turn causes troubles 
upstream, especially during flood stages, such as interference with upstream 
structures and diversions, clogging of drain and sewer outlets, water loss to 
the water-loving plants by transpiration, and increased evaporation losses. 
The problem has become so acute in some delta areas that it has become 
necessary to take steps to relieve the situation. Typical examples are the 
Pecos River above MacMillan Reservoir and at various locations on the 


‘Colorado River. So far no completely effective method has been devised to 


annihilate the water-loving vegetation and there is no practical way to pre- 


vent the development of deltas. ) 
The problem then is to convey the water and the sediment through the delta 


area without sediment deposition and with a minimum loss of precious water. 


The procedure generally followed to accomplish this is to construct a de- 


signed channel and then maintain that channel so that it will accomplish the 
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job for which it was created. The methods used in the design of a channel 
such as this have been discussed elsewhere in this paper. In the Middle Rio 
Grande Valley the Bureau of Reclamation has utilized a separate conveyance 
channel and floodway and on the Colorado River a combined floodway and con- 
veyance channel, (32)(33) These channelization works have not been observed 
under flood flow conditions as yet but the ability of a river channelization to 
withstand flood flow forces is the ultimate test of their adequacy. The need 
for adequate tools with which to approach these delta problems is evident. 
How can we effectively eliminate the phreatophyte and tule growth which is 
taking such a large toll in water in an area of the country where all the water 
is needed for beneficial use? How can we create a channel through the delta 
that will prevent the adverse developments upstream and at the same time be 
economically and engineeringly sound? Much progress has been made in 
finding the solution to these problems in recent years but considerable re- 
search and development remain to be undertaken. 


Channel Degradation 


Release of clear water from a dam or introduction of irrigation return 
flows to a channel usually result in a degradation cycle if the channel consists 
of unconsolidated materials. Degradation below Davis Dam is shown on Fig- 
ure 10. It is necessary for the engineer to predict when and to what degree 
this degradation, with accompanying water surface lowering, will occur and 
how it affects the tail water curve for a dam site and installations down- 
stream.(34) No clear cut procedures are available at the present time on 
which to base such a prediction although many cases have occurred in which 
there has been degradation below hydraulic structures. The amount of de- 
gradation is dependent upon several factors which include the variation in 
size of bed material with depth, channel hydraulics, flow duration pattern, 
and the existence of downstream controls. Degradation may be limited in 
some cases only by the amount of storage behind a structure or the limiting 
slope for transport of the materials in the channel. In other cases it may be 
limited by the creation of bed armor. All present methods of computing de- 
gradation are affected by the empirical procedures for determining transport 
and limiting size of transportable materials. 

Of particular concern in the western United States is the effect of irriga- 
tion wastes and return flows on a natural channel which is being used as a 
project wasteway or drain. This usually creates a different flow duration 
pattern than the one on which the channel was developed and a greatly in- 
creased “dominant discharge.” Many such channels are changed from 
ephemeral to perennial(23) by injecting return flow from project wasteways 
or drains. 

Research that will result in more reliable methods of sampling bed mate- 
rial, computing bed material transport, and determining limiting tractive 
forces will reflect directly into the accuracy in predicting degradation. More 
research is also needed to develop procedures for computing the channel de- 
gradation. Perhaps a correlation of records of actual degradation with em- 
pirical computations for determining bed material transport could be accom- | 
plished to devise procedures for predicting degradation below structures. 
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Sediment Control 


In the development of Bureau of Reclamation projects, sediment studies 
based on sediment and water measurements are made to give the best esti- 
mate of sediment quantities involved. Then with the aid of model studies, 
the designer includes a sediment control device in his design. For the coars- 
er bed sediments these sediment control devices have included curved guide 
walls, vortex tubes, and short tunnels.(35) For finer sediments a settling 
basin type of control is usually used in conjunction with mechanical or hy- 
draulic sediment removal equipment. Two factors that influence the decision 
as to the type and extent of sediment control facilities included in the design 
of a diversion are: (1) the size and amount of sediment needing control and 
(2) the cost of sediment control versus benefits derived. 

A typical study is that of Milburn Diversion Dam on the Middle Loup River 
in Nebraska.(36) A sediment study showed the diverted sediment load could 
be as much as 600 acre feet per year without sediment control. A hydraulic 
model study developed a short tunnel type of sediment control device, based 
upon which the designers recommended a diversion including a short tunnel 
and a settling basin for desired sediment removal. The short tunnel is very 
effective because it divides the water on a horizontal plane, allowing the water 

with only light suspended load to enter the canal, while the bedload and lower 
part of the suspended load pass through the short tunnel to the river channel 
downstream. 

In those cases where excess sediment is carried into canals, various de- 
vices have been used to remove the coarse bed sediments. In the countries 
of India and Pakistan these are called excluders. Ralph L. Parshall(37) de- 
veloped the vortex tube for removing sand from canals. His studies and later 
studies by Koonsman(38) indicated that to get good action in the vortex tube, 

a velocity near critical velocity V = Ved was required. To overcome this dif- 
ficulty of requiring such a high average stream velocity to get good action, an 
activator to work in conjunction with the vortex tube was developed in the 
Hydraulic Laboratory of the Bureau of Reclamation as a part of the model 
studies for Republic Diversion Dam. This activating vane is placed directly 
over the vortex tube with tapered upstream and downstream portions so that 


the velocity is increased as it approaches the vortex tube. The necessary 


drag force is thereby created on the channel bottom to create good vortex 
action in the tube keeping the bed material in suspension so it can be drawn 
off. Where the drag forces are insufficient, the vortex tube will fill with bed 
material and the tube will become ineffective as an excluding device. A typi- 
cal design for a vortex tube with activator is shown on Figure 11. Many other 
types of sediment control devices on canals have been tried in this country 
and particularly in India and Pakistan. The success to which these work 
varies a great deal and depends upon the hydraulic characteristics of each 
case. The philosophy in this country has been to remove most of the unwanted 
sediment load with the least amount of water loss. Research along this line 
of thinking is needed to improve the design and operation of our diversions 


and canal systems. 
SUMMARY 


In the foregoing discussion some of the sediment problems encountered by 


Bureau of Reclamation engineers in the planning, design, and maintenance of 
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irrigation projects have been pointed out and some of the methods of solution 
have been indicated. The need for research in sediment hydraulics to im- 
prove and simplify present methods of acquiring basic data, computation, and 
design is apparent. There is a wealth of existing data available in the United 
States and abroad that could be combined with additional field and laboratory 
research to advance the science of sedimentology. 

In some sections of the world, where irrigation developments are centuries 
old, sediment hydraulics has been and still is one of the major engineering 
problems. This points up the need for recognizing these problems early and 
developing corrective measures to reduce sediment problems before they be- 
come enormous. 

The following is a summary of the phases of sediment hydraulics discussed 
herein where research is needed: 


Sediment sampling—Improvement and simplification of present sampling 
equipment, procedures, and techniques will reduce the cost of a sediment 
sampling program and at the same time produce more representative data. 
Work is continually being done by such organizations as the Federal Inter- 
Agency River Basin Committee to improve sampling equipment. This work 
should be continued. 


Bed material transport computations—The need for accurate and versatile 
methods of deriving unmeasured sediment transport is ever present. New 
formulas can be developed and present formulas improved so that they can be 
used with confidence within defined limitations. In developing any sediment 
transport formula, the point must be kept in mind that its usefulness depends 
upon the accuracy of the answers given when computing transport of natural 
streams. 


Transport in nonalluvial streams— Existing bedload formulas were devel- 
oped for alluvial-type channels with bed material of sand and fine gravel sizes. 
There are no such formulas for determining the transport in mountainous- 
type streams that are not strictly alluvial and have bed material ranging from 
sands to large boulders. Can a method of computation be developed to fit this 
case or should efforts be directed toward actual bedload sampling? If a bed- 
load formula can be developed, how should a representative sample of the bed 
material existing under various river stages be obtained? 


Total sediment transport—Because of the mechanical limitations of our 
present day sampling equipment and the early day developments of formulas 
to compute transport, we tend to think of sediment transport as segments 
rather than the overall total load. Total load and its size distribution is what 
must be determined in sediment studies. Inclusion of total sediment load 
figures in the published records would be a notable improvement and would 
be of much greater value to the sedimentologist. 


Stable canals and channels—The need to define the influencing variables 
that create a stable condition in a canal or channel is one of the most urgent 
needs in the field of sediment hydraulics. Various formulas have been de- 
veloped for computing stability but application of these formulas to a variety 
of field problems has indicated that further research is needed. Research 
along this line is being carried out by some organizations and agencies both 
in the field and in the laboratory. 
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Design for stability—The definition of the correct combination of variables 
that result in canal or channel stability will culminate in reliable methods for 
designing for stability. Although current design formulas and procedures ap- 
pear to conflict or are based on a somewhat different approach it is hoped 
that continued investigations and research will reduce these differences and 
result in procedures that can be combined with definable field data to give 
more accurate solutions. The criteria for the type of channel control to use 
under different conditions and the limitations of these control works should 
also be improved and be more firmly established. 


Aggradation and degradation—Very little progress has been made in de- 
veloping methods for predicting channel aggradation and degradation. It is 


believed that a combination or correlation of records on actual happenings 
with sediment transport theories could provide some procedures for predict- 
ing aggradation and degradation. These procedures will, of course, still be 
limited in accuracy by the adequacy of the bedload formulas. 


Deposition in reservoirs— Methods of predicting the distribution of the 
sediment that is to be deposited in a reservoir are still in need of improve- 
ment. A vast amount of data are available on existing reservoirs that could 
be sorted out and combined, perhaps by multiple correlation, to derive pro- 
cedures or formulas for predicting how the sediment will deposit behind a 
structure. 


Diversion structure design—The design of a diversion structure on a 
stream carrying unwanted sediments usually includes a sediment control de- 
vice. The method employed to exclude the sediments can best be determined 
by a laboratory model study. Model study results, combined with the results 
of the prototype operation, will yield even better designs in the future. 


Stabilizing canals below structures—The use of riprap below structures in 


canals is a common design practice to prevent or repair cases of erosion. 
The size and amount of riprap used are usually chosen on the basis of ex- 
perience. Research which would assist in this field is needed. 


Sediment exclusion from canals—Reducing the cost of cleaning canals is a 


- goal towards which operating engineers on all irrigation projects are continu- 


ally striving. Hydraulic sediment excluding devices are in use in some loca- 
tions but for the most part canals are usually cleaned with mechanical equip- 
ment. Development of less costly means of keeping undesirable sediment out 
of canals would be very welcome research. 


Design of canals in cohesive soils—The design of stable canals in earth 
materials to prevent scour necessitates the balancing of the scouring forces 
with the resistive forces. For noncohesive materials considerable informa- 
tion is available. There is a definite need to determine a correlation between 
some soils tests for cohesive materials with their resistance to scour by 
flowing water which can be used as a guide for designing canals in cohesive 


materials. 
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A. DH-48 Sampler with nozzle adaptor for point sampling. 


B. Plastic Tube Sampler for sampling close to hard stream 
bottoms and for large sediment. 


Figure 2, Samplers adapted for special conditions. 


ASCE CARLSON - MILLER 953-25 


| 309 998 99 98 95 90 80 70 60504030 20 10 5 2 05 0.201005 001 
PER CENT OF TIME = INDICATED AMOUNT 


April | to July 15 
(Spring Runoff) 


DAILY DISCHARGE IN C.F.S. 


PER CENT OF TIME & INDICATED AMOUNT 
001 005010205! 2 5 10 20 30405060 70 80 9 9 2 99 998 9999 


Figure 3 -Seasonal Flow-Duration Curves, 
Son Juan River, Bluff, Utah 


HY 2 April, 1956 


953-26 


000'000") 000 ‘oo! 


‘ONITMWYS LNIWIG3S 
O3QN4adSNS ATIVG 4O SYVIA 6I JO 
SLO1d dNOYS NO G3SvVE SaAuND 


YOIN “$3NIg “YeAIY UONP UDS 
“S@AIN) Bulyoy-juawipas papuadsns jouosoas— % ainBiy 


000'o! 000"! 
AVG Y3d SNOL NI QVO7 LN3WIG]aS 


° 
> 
= 
= 
° 
o 
° 
=x 
> 
2 
o 
m 
z 
2 
aid 
om 


=|P®r on Tt -m 


° 


| ASCE 


CARLSON - MILLER 953-27 
Station: Bluff, Utah Period: April 1 to July 15 River: San Juan 
1 2 3 h 5 6 7 
4 % % 2x 2x5 
limits interval midord Ww Qs Qv disch %s disch 
0.00-0.02 0.02 0.01 40,000 3,500,000 8.0 700 
0.02-0.1 0.08 0.06 32,000 2,400,000 25.6 1,920 
0.1 -0.5 0.4 0.3 25,000 1,550,000 101.2 6,200 
0.5 -1.5 1.0 1.0 20,700 1,130,000 207.0 11,300 
1.5 -5.0 3.5 3.25 16,300 740,000 570.5 20,450 
5-15 10 10 11,700 00,000 1,170.0 40,000 
15-25 10 20 8,900 245,000 890.0 2h,500 
25-35 10 30 7,000 157,000 700.0 15,700 
35-45 10 ho 5,650 107,000 565.0 10,700 
45-55 10 50 4,720 79 ,000 472.0 7,900 
55-65 10 60 3,750 52,000 375.0 5,200 
65-75 10 70 2,940 33,500 294.0 3,350 
75-85 10 80 1,950 17,000 195.0 1,700 
85-95 10 90 1,100 7,000 110.0 700 
95-98.5 3.5 96.75 540 2,100 18.9 Th 
98.5-99.5 1.0 99.0 25 12 0.3 
99.5-99.9 0.4 99.7 8 5 
99.9-99.98 0.08 99.9% 2 2 
99.98-100 0.02 99.99 ) () 
TOTAL 5,702.5 150,394 
QvP.D. = 5,702.5 D.D. x 106 x 1.9835 = 1,200,000 (AF)/period 
QgP-D. = 150,398  D.D. x 106 = 15,950,000 tons/period 
20 percent correction for bedload = 200,000 tons/period 
Total sediment discharge = 151504000 tons/period 
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SYNOPSIS 


| Economic considerations may limit the degree of flood protection that may 
ibe provided and thus leave some areas subject to residual flood damages. 
‘Unless this limitation is fully understood by all concerned, progressive ex- 
pansion in the flood zone may eventually reach such proportions that addi- 
tional flood control measures will be required to protect the added wealth 
placed in the path of the flood waters. Flood plain zoning offers a means of 
controlling development within the flood zone to the degree that can be toler- 
ated without increasing the frequency of excessive flood damages. Zoning 
may be accomplished by local ordinance, by purchase of fee title to some 
areas and flowage easements over others, or by general education of the 
public regarding flood dangers that still exist in the area. 


INTRODUCTION 


Flood plain zoning provides a basis for coordinating man’s use of the flood 
zone with that of the river. It seeks to limit the intensity of use to those ac- 
tivities and facilities that by their very nature must be located in that region, 
and to limit other developments to the types and extent that can tolerate the 
flood potentiality that may be expected at their respective locations. The 
orocedure may be used as a Substitute for flood control measures, or as a 
supplement to such facilities. In this paper, attention will be directed pri- 
marily toward the latter use. Attention furthermore will be directed toward 
-onsideration of urban areas, rather than rural or agricultural areas, al- 
hough it is significant to note that the same principles would be applicable in 
sither case. 

In a broad sense, occupation and use of the flood zone represent an en- 
roachment on the river’s domain. This encroachment may be considered as 


Jote: Discussion open until September 1, 1956. Paper 954 is part of the copyrighted 
Journal of the Highway Division of the American Society of Civil Engineers, Vol. 
82, No. HY 2, April, 1956. 

. Prepared for presentation at the New York Convention of the ASCE, Octo- 
ber 24-28, 1955. 

. Chief of Eng. Div., Pittsburgh Dist., Corps of Engrs., U. S. Dept. of the 

Army, Pittsburgh, Pa. 
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of two types. The first is that ordinarily envisioned when the term flood plai 
is used in connection with the placement of structures and other facilities in 
the region between the top of the river bank and the landward limit of the arez 
subject to flooding. Since occupation of this type is outside of the normal flow 
channel of the stream, it generally has little effect on the flood itself, except 
insofar as the volume of overbank storage normally available for flood dis- 
charge is eliminated by structures. The principal effect of this type of occu- 
pation is that it places material wealth in a position where flooding of the 
stream will cause loss or damage to property. 

The second type of encroachment is that where the stream bank or bed is 
occupied by structures. Examples would be where the river bank is extended 
by filling in order to gain area at or above the top-of-bank level, or where 
structures, such as bridge abutments, and even buildings are so placed as to 
encroach seriously on the waterway capacity of the stream. It should be note 
that when structures, particularly buildings, are placed over or encroaching 
on the flow channel of the stream, not only is wealth placed in the pathway of 
the flood water, but the structures themselves may actually impede the move- 
ment of flood water and thus aggravate the flood condition upstream. 


History of Development 


As a background to consideration of the problem of flood plain zoning, it 
may be well to review briefly the history of development within the flood zone 
of our streams. The need for use of the stream for various functions has, 
since the earliest days of settlement, resulted in a concentration of develop- 
ment along the stream banks in urban areas. These essential uses of the 
stream are water supply, both domestic and industrial; waste disposal; trans- 
portation, both on the waterway itself, and by railroads and highways, the 
latter due to the more favorable grades available on the flood plains, particu- 
larly in the more rugged areas of the Country, such as western Pennsylvania 
and West Virginia. A fourth use, which was of considerable importance when 
the Country was young, was that of water power for saw and grist mills, the 
location of which served to concentrate community development in the area 
of the mill. 

The primary cause of encroachment on the flood channel itself, has been 
the progressive need for expansion of properties and the general congestion 
of the urban areas, which would permit expansion only in the direction of the 
stream. 

Several other factors have influenced growth of development in the flood 
zones. The pioneer, of course, lacked information regarding the flood poten- 
tialities of the stream, other than could be discerned from physical evidences 
of past floods. Areas that were subject to more or less frequent flooding 
naturally were developed with full knowledge of the flood danger. The inter- 
vals between floods of greater magnitude, however, have been the cause of 
much development in the flood zone that might not otherwise have taken place 
For example, there has always been a tendency to forget the flood as the tims 
Since the last experience lengthens. Furthermore, where properties have 
changed ownership since the last severe flood, the new owner quite often has 
acquired the property and developed it without full knowledge of the flood 
danger. In general, therefore, past history has indicated that following a se- 
vere flood, there is a period of depressed values and caution in use and 
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development of the flood zone. However, as time passes without severe flood, 
the tendency toward caution lessens and property values gradually increase. 
Through the action of these various factors and the general growth and 
expansion of the Nation, there now are immense values in physical develop- 
iment within the areas along the stream banks that in the past have been sub- 
i jected periodically to floods of various magnitudes. With the flood of 1927 in 
the lower Mississippi River, the Nation was alarmed sufficiently that Con- 
|gress, by the Flood Control Act of May 1928, authorized a project of expanded 
# Scope to provide adequate flood protection for the lower Mississippi River. 
, Approximately nine years later, in March 1936, there occurred floods in the 
j northeastern section of the United States, including the upper Ohio Valley, of 
Such magnitude as to exceed any known flood that had previously occurred in 
the respective streams. This, in turn, spurred Congress to pass the Flood 
| Control Act of June 1936, in which it recognized the flood control problem in 
|} general as one of national interest. Subsequent to these two flood control 
Acts, there have been various others and various flood control measures have 
een passed by States and local governments, with the result that the provi- 
sion of flood control facilities in areas of need has been well advanced, a 
‘number of such projects having already demonstrated their value. 


Limitations of Flood Control Facilities 


Many of these flood control projects, however, are limited in effectiveness, 
a condition of which the general public is not fully aware. Where levees and 
flood walls have been constructed to protect urban areas along the major 
streams, they generally have been designed to provide a high degree of pro- 
tection, in many instances protection against the maximum possible flood or 
at least against a flood of such proportion as would approach that limit. Where 
measures have been constructed to provide protection of that degree, it may 
be stated that for all practical purposes, the flood zone within the protected 
area has been eliminated. In certain instances, however, particularly in the 
headwater areas, it is not economically possible to provide levees or flood 
walls of such a height as would protect against a flood approaching the maxi- 
mum possible. In these cases, the flood zone remains, but it is not subject to 
inundation, except during very extreme floods. 

Where channel improvements are provided, they are effective, at least to 
some degree, on floods of all magnitudes, but it is seldom economically pos- 
sible to provide protection of this nature to the extent that the maximum pos- 
sible flood will be contained within banks. Such a limited channel improve- 
ment, however, will eliminate floods along the extreme outer edge of the flood 
zone and make less frequent the floods in the remainder of the flood zone. It 
is evident, therefore, that there remains in this case also, a residual flood 
potentiality. 

_ Where storage reservoirs are provided as a means of flood control, it is 
seldom economical to provide storage capacity sufficient to control completely 
the runoff from the maximum storm that could occur over the tributary area. 
However, the reservoir will provide complete protection immediately below 
the location of the dam for all except those flood runoffs which exceed the 
capacity of the reservoir. The effectiveness of a reservoir becomes progres- 
sively less with distance downstream from the location of the dam. This is 
due primarily to the lessening in the percentage of area controlled and to 
Jeviation in the synchronization of the maximum effect from the reservoir 
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storage with the flood runoff from the drainage area downstream from the 
dam. A reservoir provides some reduction in flood height for a considerable 
distance downstream, but due to the economic limitation of the structure, 
there remains as in the case of the limited channel improvement, a residual — 
flood potentiality. This potentiality may, of course, be further reduced by 

the addition of other reservoirs at strategic locations and by supplementation 
with flood walls, levees, and channel improvements where such can be found 
to be economically justified. 

It should be evident, however, that even when all of the flood control meas- 
ures now planned are completed, there will still be areas where complete 
flood control or flood protection has not been provided. It is in connection 
with these areas in particular, that this paper is concerned. If the previously 
mentioned factors contributing toward the occupation of the flood zone are 
reconsidered, it will be evident that many of them will still be effective. 
While domestic water supply in a modern community may require only limitec 
facilities along the immediate river bank, the need for industries to be near 
the stream as a source of water supply will continue to be of consequence 
and even more so than in the past. With the growth of modern industry, the 
demands for water supply have multiplied many fold during recent years. 
These needs may be expected to continue and even increase. On the other 
hand, modern waste disposal facilities to the extent that they must occupy the 
flood zone, may be designed in such a manner as to minimize damage from 
flood. In the field of transportation, however, particularly on the navigable 
streams, traffic has increased tremendously over recent years and may be 
expected to continue. Factors contributing toward the occupation and use of 
the flood channel will also continue to be effective if not controlled. 

In addition to these factors, the general public and even some public offi- 
cials are not fully aware of the limitations in effectiveness of such flood con- 
trol measures as have been provided. There has grown up a feeling that with 
flood control as now provided, areas that were formerly subject to flooding 
are now flood-free. Furthermore, civic organizations in their zeal to attract 
new industries and commercial enterprises to their respective areas, have 
often over-stressed the safety feature against floods. 


Need for Supplemental Flood Plain Zoning 


With these various factors still operating, there is danger that unless con- 
trolled, continued development in the flood zone will eventually result in the 
need for even further flood control measures to eliminate damages created 
by the additional encroachment within the flood zone. It may be well to stress 
at this point, that it is not intended to imply that further development within 
the flood zone should be completely prohibited. Certainly, the provision of 
flood control measures has made more useful, land which formerly had been 
of limited potentiality. For example, the development of the Point Park and 
Gateway Center at the junction of the Allegheny and Monongahela Rivers in 
downtown Pittsburgh is considered a wise investment. This development was 
not undertaken before the local officials had carefully investigated the effec- 
tiveness of the flood control reservoir system that was being provided in the 
streams above Pittsburgh. There is danger, however, that in other localities 
local officials may not be as exhaustive in their investigations when consider 
ation is being given to rehabilitation of areas that have been blighted by flood: 
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Some form of control of development along the streams where flood pro- 
jtection has been provided is therefore considered essential. As a basis for 
uch control, it will be necessary that the entire area still subject to flooding 
e classified with respect to frequency and depth of flooding. Recognition will 
ave to be given to the necessity for certain facilities and developments being 
fat or in close proximity to the river bank. The advantages of being near the 
tream will have to be weighed against the probable damages that would re- 
lt from such location; for example, in the case of industry, the additional 
icost of water pumping and additional cost of utilizing navigation facilities will 
ave to be balanced against the potential damages to the plant as a result of 
ood. In effect, it will amount to establishing the types and degrees of devel- 
pment that can be permitted on each classification of property without in- 
reasing the frequency of excessive flood damages. For example, former 
aste land might be made useful for agricultural purposes or for storage of 
i seriais not subject to severe damage by flood. Areas formerly used for 
e latter purpose might now be made available as warehouse or mill sites. 
Former low-class commercial areas may become high-class commercial 
gees with modern stores and office buildings. 


: Methods of Provision 
'_ There are several methods by which flood plain zoning may be attained. 

t may be done by means of local ordinance, which probably would not involve 
payment of damages to the property owner. Without such payment, however, 
it is doubtful that it could be made effective in the removal of any existing 
developments, but might serve as a guide to the future. Another possibility 
would be the purchase of fee title to some areas and easements over others. 
A third possibility would be the general education of the public and its en- 
couragement to follow an adopted community plan for any further develop- 
ments in the area. This would at least make certain that the owner is ap- 
prised of the flood danger before he proceeds with development. 


CONC LUSION 


The need for some action, such as flood plain zoning, in areas that have 
been benefited by a limited degree of flood control, is considered essential 
before expansion in the flood zones and flood channels has progressed to the 
extent that zoning becomes impractical. This need has been at least partially 
recognized in some of the recent reports on flood control measures prepared 
oy the Corps of Engineers, wherein it has been specified that zoning or some 
method of assuring retention of the flood discharge capacity of the stream ibe 
included as an item of local cooperation. er 

In closing, the writer wishes to emphasize that the opinions expressed 
1erein are his own and not necessarily those of the Corps of Engineers. 
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Discussion of 
“THE IMPORTANCE OF FLUVIAL MORPHOLOGY 
IN HYDRAULIC ENGINEERING” 


by E. W. Lane 
(Proc. Paper 745) 


STAFFORD C., HAPP.!_ Fluvial morphology is a field in which greater 


| coordination between geomorphology and engineering should be mutually ad- 
| Vantageous, and the author is to be commended for his accomplishments in 
_ furthering such coordination. His summaries and examples of engineering 


applications of geomorphic concepts should be helpful both to engineers and 
to geologists. This discussion is devoted to some differences between the 
engineering and geologic viewpoints, which have confused previous attempts 
at coordination and hence deserve attention by the author’s readers. 
Difference in time factors is a major problem. The geologic concept of 


erosion cycles involves time measured in millions of years. The average 


rate of stream bed lowering, in a cycle of such duration, might be of the or- 
der of one foot, or less, in a thousand years. Similar magnitudes are sug- 
gested by the amount of erosion accomplished since the earlier stages of 
Pleistocene glaciation. Engineering interpretations of fluvial morphology, on 
the other hand, are usually based upon relatively short records in which the 
slow changes involved in normal geologic cycles are obscured by shorter- 
period fluctuations resulting from normal variations in stream discharge, 
lateral channel shifting, cyclic or other climatic changes, or the effects of 
human activities. Normal shifting of sediment by stream meandering, and 
thick deposits of unconsolidated alluvium which fill many valleys as a result 
of Pleistocene glaciation, are other factors contributing to the complexity of 
engineering problems of fluvial morphology. The geologic concepts of youth- 


‘ful, mature and old age stages in stream development are therefore only of 


limited descriptive value from the engineering viewpoint, and do not indicate 
a sequence of changes within time intervals of engineering significance. 

The author’s classification of stream profile changes illustrates the dif- 
ference between engineering and geologic viewpoints, for it omits diastrophic 


‘uplift or crustal tilting, the primary factor in erosion cycles of geologic 


terminology. Aside from this omission, the classification appears somewhat 
inconsistent by grouping into 2 classes the changes due to 3 factors of the 
basic equation (Q,d ~Qys), whereas 4 classes all depend on changes in the 
fourth factor of the equation. Changes in distribution, as well as volume of 
water discharge, may also deserve mention as factors in the classification. 
It is suggested that grouping degradational and aggradational effects, but 
itemizing changes due to individual factors separately, would be more con- 
sistent with the form of the author’s basic equation. This would facilitate 
brief descriptive identifications of the classes. Inclusion of the primary geo- 
logic factor of crustal changes would make the classification more complete 


1. Geologist, U. S. Atomic Energy Comm., Grand Junction, Colo. 
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from the geologic viewpoint. Such an alternative classification might be as 
follows: 


I. Degradation, due to: 

Increased slope caused by diastrophic uplift or crustal tilting. 
Increased slope caused by lowered base level. 

Increased slope caused by upstream shift of base level. 

Increased volume or more effective distribution of water discharge. 
Decreased sediment quantity. 

Decreased sediment size. 


AmOOW > 


it; Aggradation, due to: 

Decreased slope caused by diastrophic subsidence or crustal tilting. 
Decreased slope caused by rise in base level. 

Decreased slope caused by downstream shift of base level. 
Decreased volume or less effective distribution of water discharge. 
Increased sediment quantity. 

Increased sediment size. 


SAHOO W Pp 


If the revised classification here suggested does not meet with the author’s 
approval, it would undoubtedly be helpful to many other readers if he would 
offer, in his closure, brief descriptive designations for his various classes of 
profile changes. 

The premise of a “graded stream in equilibrium,” on which the author’s 
classification is based, might be improved by substitution of the term “poised” 
in place of “graded.” “Poised” was proposed(!) for streams which satisfy the 
author’s requirement that “the net amount of change is not sufficiently large 
to be detected by quantitative measurements.” It is an appropriate and con- 
venient term, lacking the connotation of confusion long associated with geo- 
logic use of the graded stream concept. 

The term “graded stream” was proposed by W. M. Davis(2) to identify the 
condition of balance between degrading and aggrading. The suggested termi- 
nology was credited to G. K. Gilbert(3) who had previously described a stream 
in equilibrium as one having a supply of debris equal to its capacity, and 
therefore tending to establish a uniform slope by building up the gentler slopes 
of its bed and cutting away the steeper slopes. Davis elaborated the concept 
to include maintenance of an “essential balance” between degradation and ag- 
gradation throughout the long stages of maturity and old age in the erosion 
cycle. The apparent anomaly of long-continued degradation by a stream in 
essential balance between degradation and aggradation was explained by postu- 
lating that declining sediment contribution from the drainage area would per- 
mit renewal of degradation despite the concurrent slope flattening. This ar- 
gument is unconvincing, however, for why should the separate processes of 
drainage area erosion and river bed degradation proceed at equivalent rates? 
If sediment production continued increasing, aggradation might continue for a 
long time; and if sediment production should later decline, why should this 
decline cease, or become slower, when the resulting degradation had just 
balanced the previous aggradation? Or, if sediment production should soon 
decline, after attaining a balanced condition, why might equilibrium not be 
destroyed by renewed degradation? 

The difficulties in geologic usage of the term “graded stream,” appear to 
have arisen from the unfortunate combination of two distinct concepts in the 
definition by Davis. One concept is that of equilibrium between aggradation 
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and degradation; the other is of equilibrium between the capacities of adjacent 
stream segments, resulting in uniformly smooth profile slopes. The latter 
| concept involves grading in the common engineering sense, such as grading 
| of a roadway by smoothing of abrupt irregularities or excessively steep 
| slopes. According to this concept equilibrium must be a precise condition of 
balance and any deviation from such balance initiates a reaction tending to 
| restore the balance; but it is a balance between rates of either aggradation 
or degradation, not a balance between aggradation and degradation. 
Equilibrium between aggradation and degradation is not dependent ona 
| precise or precarious balancing between capacity and load. Aggradation oc- 
curs if load exceeds transporting capacity; degradation occurs if erosive 
capacity exceeds resistance of the bed. The capacity required for degrada- 
| tion is greater than that required to prevent aggradation, however, because 
bed erosion is more difficult than transport of the resulting sediment after 
entrainment. Streams seldom carry a capacity load of bed sediment because 
_ they are unable to produce such material as rapidly as they can carry it away. 
_ The bed load deposits which normally occur along a stream channel do not 
usually indicate aggradation, but are essentially transitory deposits formed 
because sediment transportation is an intermittent process. The magnitude 
of these transitory bed load deposits depends only partly on the rate of bed 
erosion, and more on rates of sediment production by weathering, rain wash, 
creep and other mass movements, and especially stream bank erosion. Ifa 
stream, in which load and transporting capacity were just balanced, should 
have its capacity increased, the amount of bed load deposits would be reduced, 
but degradation would not necessarily follow. There would be some range of 
slope and volume—the controlling factors—within which the capacity would be 
more than adequate to prevent aggradation, yet not adequate to accomplish 
appreciable degradation of the bed. 
Equilibrium between aggradation and degradation may persist over a com- 
paratively wide range of slopes and volumes if a stream is flowing on hard, 
resistant bedrock, but the range will be much narrower if the stream bed is 
underlain by thick, unconsolidated sands. Engineering projects and model 
studies are largely concerned with the latter type of stream conditions, and 
perhaps for this reason engineering literature commonly treats the equilib- 
rium between aggradation and degradation as a precarious, transitory condi- 
tion. Geologic erosion concepts are concerned chiefly with uplifted mountains 
or plateaus in which hard, resistant bedrock is common, and hence equilib- 
rium is a relatively more stable condition from the geologic viewpoint. The 
author has apparently overlooked this factor in the conclusion, stated on 

page 745-5, that “neither in very short nor in very long periods may natural 
streams be considered to be in equilibrium.” 

Stream gradation by development of smooth slopes, as contrasted to 
equilibrium between aggradation and degradation, appears to be a natural re- 
sult of normal stream development. Youthful streams have excess energy 
due to high gradients, and erode their beds actively. Sections of the streams 
flowing over weak rocks are eroded deeper and reduced to lower slopes, 
while those sections flowing on hard rocks retain steeper slopes or even 
waterfalls. In this stage differential erosion tends to increase the irregulari- 
ty in slopes, but the trend toward disparity in slopes is eventually reversed 
by accumulation of bed load materials. 

The coarser and heavier bed load particles are moved only at higher 
stages; during falling stages they are deposited first and concentrated near 
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the base of the sedimentary accumulation. The larger and coarser particles 
are also moved less frequently and less rapidly in the flattened reaches on 
weak rock, and accumulate more thickly in such places. As this process con- 
tinues, greater and greater floods are required to entrain the bed load ac- 
cumulation and expose the underlying bedrock; and the weak rocks are subject 
to erosion less frequently and for shorter periods. Thus erosion of the 
weaker rocks will be retarded relative to that of hard rocks, and the bedrock 
surface will tend to become more uniform in slope. This is a process of 
stream gradation, by which the transporting capacity of adjacent stream seg- 
ments is balanced. Such a graded stream may also attain an approximation, 
or at least an apparent approximation, to equilibrium between aggradation and 
degradation, because stream bed erosion is normally so slow. The bed load 
accumulation does not greatly retard lateral meandering, and a stream may 
effect many changes in alinement, and perhaps even widen its valley con- 
siderably during the time required for relatively little degradation. 

The writer would prefer to restrict the term “graded” to a stream which 
had attained a smooth slope by elimination of irregularities in its profile, and 
which has thereby attained approximate equilibrium of transporting capacity 
between adjacent segments. This restricted concept would eliminate some of 
the difficulties which the author has recognized in his description of the 
equilibrium condition, but it would not be in accord with either the definition 
originally introduced by Davis, or the apparent intent of the author. Substitu- 
tion of “poised” for “graded” would, however, free the author’s classification 
from the ambiguities associated with the latter term. This may appear to be 
a minor detail, but it should be helpful in promoting clearer thinking. Those 
responsible for major river engineering should never forget that a poised 
condition may persist for a long time, yet be no indication of the effects of a 
flood greater than any of previous human record. Wider appreciation of this 
limitation of present data might even be helpful toward more adequate studies 
of long-term river trends, for which much of the evidence must be obtained 
from alluvial valley deposits. 
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Discussion of 
“GROUND WATER PHENOMENA RELATED TO BASIN RECHARGE” 
: 


by Paul Baumann 
| (Proc. Paper 806) 


| NORMAN H. BROOKS,* J.M. ASCE.—Mr. Baumann has presented an ex- 

| cellent explanation of various physical phenomena encountered in ground 

| water recharge. In this new and important phase of ground water hydraulics, 
| Mr. Baumann has again made a significant contribution. 

The writer is particularly interested in the section “Effect of Spreading 
Mounds on Channel Mounds.” The writer concurs with the author in the basic 
_ conclusion that the recharge from a stream channel does not prevent the re- 
| charge from an adjacent spreading ground from crossing under the stream. 

In fact, as Mr. Baumann states, the recharge from the combination of stream 
channel and spreading ground is essentially the same as it would be from a 
stream of the same width as the combination. 

However, the writer would like to point out some limitations in the author’s 
analysis leading up to this conclusion and his corresponding equation (1). In 
the first place, each of the author’s flow nets (Fig. 3a and 3b) shows a zone of 
no flow which is roughly triangular in shape in the lower central part of the 
previous stratum. With the usual assumption that ground water flow follows 
Darcy’s law, there is no theoretical justification for departure from conven- 
tional flow-net analysis; hence, the center streamline should go all the way 
to the bottom of the pervious stratum before dividing at the stagnation point. 
The entire region should be included in the flow net. 

It might be noted that Mr. Baumann’s flow nets would certainly be correct 
if the wedge-shaped sections, which he has in effect excluded from the flow 
net, were impervious. But since they are not, there will be flow in them and 

the limiting streamlines cannot possibly curve downward away from the cen- 
ter line to the impervious base as the author has shown them. 

The overall effect of this approximation may not have a very significant 
effect on the results inasmuch as the flow velocities in the region he has ex- 
cluded are undoubtedly low. The approximation certainly makes the drawing 
of the flow net considerably easier. 

The writer has still another misgiving about the flow nets shown in Figs. 
3a and 3b. The author states on page 806-7: “Actually the flow trajectories 
are curved because of the sloping stream channel and aquifer, but the flow 
net resembles that shown in Figure 3a if developed into a plane.” Theoreti- 
cally, flow nets apply only for two-dimensional flow,* and the developing of a 
curved flow pattern into a plane flow net is an expedient without theoretical 
basis. Although developing a curved flow pattern into a plane flow net is very 
helpful in visualizing the flow, it is nevertheless a risky procedure for quanti- 


tative solutions of problems. 


* Asst. Prof. of Civ. Eng., California Inst. of Technology, Pasadena, Calif. 
1. It is also possible to draw a special type of flow net for axially symmetric 
flow in three dimensions, but there is no axial symmetry here. 
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Furthermore, the way in which the streamlines curve in a plan view de- 
pends upon the length of the reach of the stream from which water flows into 
the ground. If the recharging reach is very short, the streamlines can turn 
into the direction parallel to the stream quite close to the stream. If, on the 
other hand, the recharging reach is infinitely long, then the streamlines for 
the ground water flow can never become parallel to the stream. 

In making the analysis which leads to equation (1), Mr. Baumann has 
avoided the last difficulty by assuming the underlying impervious stratum to 
be horizontal and the flow to be two-dimensional. However, in drawing the 
various flow nets which led the author to equation (1) he must have assumed 

_some kind of downstream control. Without some control point at a finite 
distance, it is impossible to establish steady flow. 

To make an approximate analysis of the effect of the distance to the con- 
trol it will be assumed that the control consists of a horizontal discharge face. 
The flow net in the vicinity of this control, shown in Fig. A, is simply the 
well-known Kozeny solution“ consisting of confocal parabolas. The line of 
seepage in this case is often called the “basic parabola.” This theoretical 
solution is quite good except in the immediate vicinity of the spreading ground. 
There the curvature will be reversed because the water must flow downward 
initially from the spreading basin, and then turn either to the right or the left. 
To get an estimate of the maximum possible rate of discharge which can flow 
laterally through the ground toward the control, an arbitrary point P on the 
line of seepage in Fig. A may be assumed to be at the level of the water sur- 
face at the edge of the spreading ground. 

By the properties of a parabola the length of the discharge face, a, is given 
by the formula 


iT 
Do| 


(VIF +D' - 1), 


a 


where D is the height of point P above the impervious base and L is the hori- 
zontal distance from P to the discharge face. Since the discharge, q’, is 
simply 2Ka (K is the permeability), the rate of flow can be expressed in terms 
of L and D as: 


q=K (VL? +D* -'L) 


Now the acceptance rate q as used by Baumann is equivalent to 2q/W, wherein 
the factor of 2 accounts for flow to either side of the spreading ground. Con- 
sequently, the maximum possible acceptance rate may be given dimension- 


lessly as 
4_5D V\ z) L 2 
= 22( is bli yy fs 
Perhaps it would be more precise to use an inequality sign in equation (2), | 


inasmuch as the quantity of flow represented by the right-hand side is 


2. See “Seepage Through Dams” by A. Casagrande, Contributions to Soil 
Mechanics, 1925-40, Boston Soc. of Civil Engineers, p. 310. 


DISCUSSION 


Fig. A - Kozeny's basic parabola solution for horizontal discharge face. 
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undoubtedly somewhat larger than what could possibly actually occur, because 
the point P is taken to be at the same height (D) as the top of the spreading 
ground. This condition is impossible because of the initial downward flow, 
but the use of D still gives an upper bound for the lateral rate of flow and 
hence the acceptance rate for the spreading ground. 

In Figure B a comparison is made between the author’s equation (1) and 
the writer’s equation (2). It is apparent that for many values of L/D the 
author’s equation (1) gives a value of the acceptance rate which exceeds the 
capacity of the pervious stratum for transmitting the flow horizontally to the 
control. For example, if L/D = 10 and D/W = 2, then the upper bound by 
Eq. (2) is q/K = 0.20; on the other hand, the author’s Eq. (1) gives q/K = 0.63, 
a value which is at least three times too large for this case. In fact, the dis- 
tance to the control (L) must not be more than about twice the depth (D) if the 
value of q found by the author’s equation (1) is to be consistently less than the 
maximum possible as indicated by equation (2). It appears then that the dis- 
tance from the spreading ground to the control is a very important variable 
in the problem which must be taken into account. 

It may be noted from equation (2) that for some combinations of variables, 


the calculated value of z will exceed 1. This is physically impossible, be- 


cause the acceptance rate can never exceed K, which is the rate correspond- 
ing to downward percolation without hydrostatic pressure. Consequently, in 
these cases where q/K > 1 by Eq. (2), the recharge mound cannot possibly 
rise up to the spreading ground level, but reaches its ultimate height some 
distance below the spreading ground, as pointed out by Mr. Baumann on page 


806-3. 


FREDERICK L. HOTEs,! M. ASCE.—Ground water problems of all types 
are increasing throughout the United States. Although there are many com- 
plex phenomena connected therewith that engineers and hydrologists would 
like to know more about, none are more important than that of basin recharge. 
Recharge, both natural and artificial, places definite limits on the safe yield 
of ground water basins on a sustained basis. The author has set forth some 


- of the phenomena connected with recharge problems for both free and con- 


fined aquifers, and has indicated some of the practical limitations learned 
from personal experience. 

This discussion will be limited primarily to topics associated with sea 
water intrusion into pressure aquifers. Several research projects conducted 
at the University of California, Berkeley, California, during recent years 
have contributed information of considerable assistance in the interpretation 
of ground water behavior in coastal regions, and on the effectiveness of in- 
jection wells in repelling sea water intrusion. 


Model Studies 


In addition to the contract with the Los Angeles Flood Control District 
mentioned by the author, the California State Water Resources Board also 
contracted with the University of California to conduct laboratory and model 
studies of sea water intrusion. It was intended that these studies would sup- 
plement and aid in the prototype experiments being performed by the Flood 


Control District. 
1. Assoc. Prof. of Irrig. Eng., Univ. of California, Berkeley, Calif. 


955-12 HY 2 April, 1956 
Model Aquifer 


Following a thorough search of the literature, (2) a sand and lucite model 
was constructed. A detailed description of the model and the experiments, 
together with the results of research on impermeable cutoff wall materials, 
has been previously published, (3) and hence only a summary will be pre- 
sented here. Figure 1 is a schematic sketch of the apparatus. 

The main chamber, representing an idealized pressure aquifer was 4 feet 
long, 6 inches high, and 3 inches wide. It was packed with a commercially 
available washed quartz sand, Del Monte 30 mesh. This sand has an effective 
size (10 percent) of 0.22 mm, a 50 percent size of 0.28 mm, and a uniformity 
coefficient (60/10) of 1.3. When packed within the model its coefficient of 
permeability averaged 0.032 cm/sec (680 gallons/square foot/day). The salt 
water was introduced through the chamber on the left and the fresh water 
through the chamber on the right. Both fluids were supplied through constant 
head flow regulators. 

The piezometric tubes were used to establish the piezometric surface un- 
der any desired condition flow. The model injection wells were located only 
along the near side of the aquifer. They were made of split glass tubing one 
fourth inch in diameter. The near face represents a plane of symmetry, nor- 
mal to the coast line, and through the centerline of one of a line of injection 
wells of equal spacing parallel to the coast. The far face represents the 
plane of symmetry normal to the coast and midway between two injection 
wells. No flow crosses these planes of symmetry in a uniform injection pro- 
cess, so they are represented by solid lucite walls in the model. 


Model Fluids 


The salt water used in the model was adjusted to a specific gravity of 
1.100. Both calcium chloride and sodium chloride solutions were successfully 
used. This procedure increased the model flow rate and generally speeded up 
an otherwise very slow action. The viscosity of these solutions was about 30 
percent greater than that of sea water. This discrepancy did not affect mea- 
surements when the sea water wedge was stationary, but did have some effect 
on rates of wedge movement. 

A novel idea that produced startling and successful results was the use of 
fluorescent dyes to aid in distinguishing the fresh from the salt water. A 
convenient concentration of ordinary dyes does not show up well in the model 
because of the thin layer of solution visible between the transparent sides and 
the first layer of sand grains. An ultraviolet light source was also used to 
produce a strong, positive response from the thin solution layer next to the 
wall. Rhodamine B in the saline water and Fluorosol in the fresh water gave 
good contrasting colors of red sea water and blue fresh water. 


Results 


Harder,(3) by an adaption of Muskat’s (4) approach to the solution of seep- 
age through an earth dam on an impervious foundation, derived the following 


expression for the fresh water flow and length of intruded saline wedge under 
equilibrium conditions. 


a =(1/2)(s-1) 2 . 
Eq. (1) 


| 
| 


: 


: 
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in which 

q = seaward fresh water flow per foot of ocean front 

S = specific gravity of sea water 

M-= thickness of aquifer, down to the lowest depth which must be protected 


T = aquifer transmissibility for a 100 percent hydraulic gradient 
L = length of sea water wedge, from ocean outlet to the toe 


It is to be noted that the flow is independent of the depth of the aquifer below 
sea level. 

Figure 2 shows the basic relationships involved and the confirmation of the 
equation by experimental results, shown plotted as circles. Fresh water is 
flowing from right to left. The same equation was obtained from potential 
theory by assuming a parabolic interface between fresh and salt water. It 
applies only to an idealized aquifer that would never be found under natural 
conditions. However, when applied with good engineering. judgment it can be 
quite useful in the study and solution of practical seawater intrusion problems. 

Figure 3 is a composite photograph showing the postion of the wedge for 
three different values of fresh water flow. The dark vertical lines are caused 
by the metal bracing of the model. The shapes of the interface are typical. 
The distances shown are for a hypothetical prototype aquifer 100 ft thick, but 
are not limited to this one case. One of the great advantages of a model such 
as this is that the results can be extended by means of model laws to any pro- 
totype aquifer. A new model does not have to be constructed for each aquifer— 
the single study serves all. Of course, models of aquifers with known varia- 
tions in permeability might justify individual models for the study of peculiar 
conditions. 

Figure 4 shows the behavior of the sea water wedge when water is injected 
through a line of wells 1000 feet inland, and spaced at 500 foot intervals. In 
this series there is no flow to or from the inland area. Hence all the injected 
water escapes to sea, and resists the intrusion. In the series on the left the 
injection rate is such that it would be sufficient to hold the wedge steady at 
2000 feet. Since the flow is not sufficient to hold the wedge seaward from the 
well, however, sea water moves between the wells and enters the inland re- 


- gion as a density current. On the series on the right the injection rate has 


been increased threefold. This rate is sufficient to hold the wedge seaward 
from the well at about 667 feet from the ocean outlet, so that part of the sea 
water wedge inland from the well is completely cut off. 

An approximation of the prototype time involved in these photographs will 
indicate the value of the model in speeding up an inherently slow process. 
The Silverado aquifer is approximately 100 feet thick. Its permeability as 
obtained from an informal conversation with the author is 25 feet per hour 
(4500 gallons/square foot/day). The prototype—model time scale ratio is 
therefore approximately 2900 to 1; or, one minute of model time represents 


2900 prototype minutes (approximately two days). The time after start of in- 


jection for each photograph in Figure 4 is as follows. 


Frame No. Left Right 
1 er (Before injection) 
2 6 months 19 days 
3 2.5 years 4.7 months 
4 5.5 years 3.5 years 
5 8.1 years 10.8 years 
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Figure 5 more nearly represents conditions at the Manhattan Beach test 
site described by the author. The series of Figure 5 (a) shows the sea water 
wedge intruding from left to right under an inland overdraft condition. Fresh 
water is being withdrawn from the right at a prototype rate of 900 gallons/ 
day/ft. It should be noted that the wedge interface is no longer concave down- 
ward, but, once inland flow is established, closely approximates a straight 
line, with localized variations at top and bottom. The interface lengthens as 
it advances inland. The toe velocity is rapid at first, then gradually de- 
creases, and asymptotically approaches the average velocity determined by 
the rate of withdrawal and the effective cross sectional area. 

This straight line interface between two fluids with different densities and 
moving through porous media has also been demonstrated by Naor\”/ ina 
Hele-Shaw model. Naor hypothesized that the interface is actually the result 
of two superimposed linear motions. One is the horizontal movement of a 
fluid of single density and gives a vertical front. The other motion is that 
investigated and described by Keulegan.(6) Keulegan observed the motion of 
the interface between two fluids of different densities in a porous media 
originally separated by a vertical gate. Upon removal of the gate a tilting 
motion occurred. The interface appeared to rotate about the mid-height, and 
to remain a straight line during the motion. The superposition of this linear 
rotation on the linear horizontal displacement results in the straight line in- 
terface observed by Naor and seen in Figure 5 (a). 

Figure 5 (b) shows the effect of injection through a well located 2000 feet 
inland at a prototype rate of 415 gpm for wells spaced at 500 foot intervals. 
This rate is sufficient to provide for the overdraft of 900 gals/day/ft, and 
the seaward flow of 293 gals/day/ft required to stabilize the wedge at 2000 
feet inland from the coast. 

It should be emphasized at this point that the toe of the wedge would prob- 
ably be about 2250 feet inland midway between the wells. The result of this 
tendency of the wedge to work inland between wells is clearly seen in the 
author’s Figure 8. It is recommended that the line of injection wells be lo- 
cated at least one half well spacing inland from the allowed maximum intru- 
sion of the wedge. The flow from the injection wells tend to equalize at this 
distance from the wells and a relatively uniform flow rate per foot of coast 
is obtained. If the wells remain at 2000 feet, a greater injection rate is re- 
quired to keep the wedge at the midpoint between wells from passing beyond 
the line of wells. This means greater loss of injected water, at no benefit. 

One of the most important results of the test was to show that the wedge 
would be stabilized at the same point, for the same rate of seaward flow, no 
matter how far the wells were spaced inland, as long as they were at least 
one half well spacing beyond the maximum allowable intrusion distance. 
Tests were also made on wells penetrating only the upper and/or the lower 
half of the aquifer, and from distant inland areas. All results confirmed 
theory—the location of the source of the seaward flow was immaterial, within 
the limits previously mentioned. 

The distortion of the vertical scale is tenfold in Figure 5 as compared 
with the Silverado aquifer. This vertical exaggeration should give a fair in- 
dication of the relative importance of the diffusion zone at the interface. The 
interface between the sea water wedge and the overlying fresh water re- 
mained quite sharp in the model and such should be the case for the prototype. 
Even the fresh water cone from the injection well showed little mixing with 


the fresh water originally in the aquifer, but tended to move out by a replace- 
ment process. 
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In addition to the still photographs, a 16 mm technicolor motion picture, 
with titles, was made by means of lapse time photography. This film is 
available on loan from either the California State Water Resources Board or 
the University of California. 


Injection Wells for Fresh Water Barrier 


The use of injection wells for the formation of a “fresh water barrier” 
against sea water intrusion is a tempting concept; and perhaps a practical 
and economical alleviating solution to the problem in some cases. However, 
certain basic principles must be born in mind if this method is to be the most 
economical solution. There are four broad questions to be answered. 


1) How much water should be injected? 

2) Where and how should the wells be located? 

3) How can the required amount be injected into the aquifer? 
4) Is this the most economical solution? 


How Much Water? 


The answer to this question lies in the principles of the Conservation of 
Matter, and of Continuity. Sea water intrusion is caused by an inland over- 
draft condition. The laws of nature require that a small amount of fresh 
water continually flow out to sea to repel the sea water. This amount can be 
estimated from equation (1) for any specified allowable location of the toe of 
the wedge. Whenever this seaward flow is reduced the wedge will move in- 
land. If the flow is not only reduced, but eliminated, and perhaps even made 
a negative quantity by inland pumping, sea water will flow inland. Wedge 
stability can be achieved only by furnishing a seaward flow of fresh water. 
Water in sufficient quantity will not flow seaward unless the overdraft causing 
the invasion is eliminated. In simple language, the total amount of water re- 


quired is equal to the sum of the seaward stability flow plus the existing over- 


draft. This amount of fresh water can be supplied by several methods or 


combinations that will be discussed later. 


Location of Wells 


_ The writer is not as optimistic as the author regarding the status of the 
sea water trapped inland from the line of injection wells, and shown in 
Figure 6 (b) of the author. The writer believes that this salt water will con- 
tinue to move inland and will begin to appear in wells that have previously 
been in a fresh water zone only. Figure 8 of the author indicates that this is 
exactly what is taking place at Manhattan Beach. The iso-chlors of Figure 5 
that are inland from the injection wells are even further inland in Figure 8. 
Any wells pumping near the wedge will induce some of the salt water to flow 
to them. This salt water will mix in the well with fresh water. It cannot be 
removed except by pumping, as long as the injection wells remain as shown. 

The writer recommends that injection wells to be used for sea water re- 
pulsion be located inland from the most advanced intrusion of the wedge toe. 
The same amount of injected water is required. This location has the great 
advantage of immediately causing all the wedge to begin a retreat to the sea. 
All pumping wells should improve in quality and none should become more 
saline. All the invaded sea water will be flushed out into the ocean and the 


aquifer will once again be fresh. 
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Problems of Injection 


It is well recognized that it is much more difficult to get water into the 
ground by pumping than out of it.. The large scale tests reported by the author 
are certainly some of the most significant yet undertaken. Experiments con- 
ducted on a well field (Figure 6) at the University of California have yielded 
interesting results on injection rates for both clear water and water contain- 
ing 20 percent of primary settled sewage.(7, 8) Clear water injection rates 
of 8.4 gallons per foot of aquifer with a permeability of 1900 gallons per 
square foot per day were successfully obtained. The sewage injection trials 
were also successful, but a redevelopment technique, reported in detail in 
the references, is vital. Of especial importance was the finding that no coli- 
form organism appeared in the observation well some 190 feet away. Some 
coliform organisms did appear at a distance of 100 feet. Increasing the in- 
jection rate beyond a certain limiting value did not increase the travel dis- 
tance of coliform organisms. These studies confirmed again the inherent 
natural filtering action of an aquifer. If sewage rather than domestic water 
can be used for recharge by injection, considerable savings in cost can pos- 
sibly be realized. Further studies of this aspect are believed warranted. 


Economical Solutions 


The quantity of water required to repel sea water intrusion can be esti- 
mated. Usually the total water use in a critical area cannot be expected to 
decrease, but on the contrary, present estimates usually predict even greater 
water use in such areas. One part of the solution is obviously then either ad- 
ditional importation of water, reclamation of waste water, or both. Once this 
supplemental supply is brought into a overdraft region, its disposition should 
become a subject of economic analysis. The sea water intrusion can be 
eliminated by either of the following methods. 


1) Injecting the water into the aquifer. 

2) Supplying the water directly to the users through a surface distribution 
system. This reduces the draft on the aquifer and allows the required 
seaward flow to take place. 


The total amount required is the same in each case. However, Method 
One requires the water to be pumped in, and then pumped out again. Obvious- 
ly this costs money, and probably should not be justified unless this is a 
cheaper method than the construction of a surface distribution system. Thus 
the value of the pressure aquifer is as a conduit. It provides no storage, but 
provides a “spigot (well)” for each overlying land owner. 

In the case of an unconfined aquifer, there are additional benefits to be 
credited to Method One by its keeping sea water out of the reservoir during 
a period of dry years, when the water users may lower the ground water level 
below sea level. | 

The value of the aquifer in acting as a purifying filter for injected sewage 
must also be considered in cases where waste water is available in sufficient 
quantities. 

Ground water problems are inherently more difficult to grasp than surface 
water problems, because in most instances they cannot be actually seen, and 
must be dealt with in abstract terms, mathematics and drawings. It is only 
through studies and tests such as the author has described, that can practical 
solutions be found. His contributions to the literature are appreciated. 
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JAMES A. HARDER, !—Mr. Baumann’s description of the ground water re- 
charge operations being undertaken by the Los Angeles County Flood Control 
District is a welcome contribution to the literature dealing with water con- 
servation—a subject not only of present vital importance to such semi-arid 
regions as Southern California, but one sure to be of increasing importance 
in other parts of the United States. This writer would beg to take exception, 
however, to certain of his statements minimizing the importance of the 
“saline wave.” boas 

As shown by the iso-chlors in Figure 5 of the paper, at the time injection 
was begun the injection site was underlain by water having a chloride ion con- 
sentration about 90% of that of sea water. This fact lead to a fear that if this 
jarticular site were to be used, the sea water thus trapped would cause 
rouble as it moved inland. Accordingly, the State of California (which pro- 
vided the initial financial support) required the District either to use a site 
nland from the intrusion or to obtain waivers from major water users 
xxempting the state from responsibility for damage done. The waivers were 
btained, and the program proceeded. é 

Returning to the iso-chlors drawn in Figure 5 of the paper, an estimate 
nay be made from them of the total chloride content of the water which lay 
yehind the injection well site before injection. This can be done by making a 
rofile of the concentration vs. distance inland from the injection well site 
nd applying graphical integration. Based on profiles through the end points 
nd the center of the site, and assuming an aquifer thickness of 100 feet and 
. porosity of 30%, the average quantity of chloride ion existing behind the 


. Asst. Research Eng., Univ. of California, Berkeley, Calif. 


955-24 HY 2 April, 1956 : 


4,000 foot long site is thus computed as being about 20 tons of Cl” per foot of 
front, or about 80,000 tons altogether. This amount of Cl~, by simple arith- 
metic, can contaminate 240,000 acre feet of water to a concentration of 250 
parts per million, a generally accepted maximum for domestic use. 

Of course there is a large amount of fresh water yet remaining in the in- 
land portion of the aquifer which might be available for dilution of the saline 
wave. From Figure 4 of the paper the distance from the injection well site 
to the Inglewood Fault, the inland boundary of the basin, is about six miles. 
Assuming that the aquifer thickens to 150 feet and continues to the fault, 
there are 150,000 acre feet of fresh water within the 4,000 foot wide strip be- 
tween the injection well site and the fault. The above computed weight of 
chloride could contaminate this water to a concentration of nearly 400 parts 
per million if it were uniformly distributed. All of this fresh water will not 
be available for dilution-duty, however, as inland water users may be ex- 
pected to liberate as much as possible of it from the burden of acting as a 
diluting agent. 

Furthermore it does not appear that the natural recharge will contribute 
sufficient fresh water to be significant in the diluting effort. On the basis of 
its area relative to the basin, this strip behind the injection site should expect 
to receive no more than 2,000 acre feet per year of new water from natural 
recharge. Artificial recharge from the injection well site is contributing a 
like amount, from the authors data. 

There remains the hope that the fresh water which overlies saline water 
in the “wave” or “trapped in depressions along the bottom of the aquifer” 
could still be recovered free of salt by means of partially penetrating wells. 
This is possible, but only under severe restrictions on the average pa 
rate. For example, a calculation based on experimental data from Muskat 
will show the maximum drawdown at which a well penetrating the top fifty 
percent of the fresh water zone of an aquifer can operate without “coning.” 
Assuming that there is a layer of sea water lying at the bottom of a 100 foot 
aquifer, a drawdown of only 3.5 feet will induce the salt water to rise ina 
“cone” and enter the bottom of the well. Obviously this is a much smaller 
drawdown than is usually encountered in pumping wells supplying a normal 
yield. 

An idea of the past, present, and future status of sea water intrusion in 
Manhattan Beach may be had from Figure 5(b) of the discussion by Professor 
Hotes inthis issue. This is a series of photographs taken by the writer of 
sea water intrusion into a model aquifer constructed at the University of 
California, Berkeley.“ Overdraft and injection rates were set up to approxi- 
mate conditions at Manhattan Beach; and with scale factors adjusted so that 
the four foot long model represents a section of the aquifer 8,000 feet long 
perpendicular to the coast, the total time elapsed between the first and last 
frames represents ten years. The present situation would be similar to a 
frame between numbers four and five. In frames five and six the saline wave 
is shown moving inland together with the fresh water from the injection wells 
at a rate corresponding to about 500 feet per year in the prototype. The shift 


2. Muskat, Flow of Homogeneous Fluids Through Porous Media, first edition, 
pp 480-496, McGraw-Hill, 1937. 
3. Report on Laboratory and Model Studies of Sea Water Intrusion, Technical 


Bulletin 11, Sanitary Engineering Research Laboratory, University of 
California, Berkeley 4, California. 
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| 
in the iso-chlors between Figures 5 and 8 in Mr. Baumann’s paper seem to 
indicate the same approximate rate. The wave tends to flatten with the pas- 
| sage of time. 
| CHARLES H. LEE,! M. ASCE.—A very interesting equation is presented in 
| this paper (Equation 1) expressing the relation between the average accept- 
| ance (infiltration) rate per unit of wetted area and the parameter 4 where D 
is the depth of aquifer to impervious stratum and W the width of stream chan- 
nel (or stream channel plus spreading ground where the latter is adjacent to 


| the stream). The author well points out that this parameter is of economic 


| importance. The writer also believes that it is of physical significance. In 

order to illustrate this graphically the writer has plotted the full range of -- 

for K = 1000 Meinzer Units (gals. per sq. ft. per 24 hours for gradient of 
unity (Writer’s Diag. 1). This value is equivalent to a flow of 67.5 cubic feet 


| per second per acre for gradient of unity. The curve approaches closely to 
its asymptote, q = 67.5, for values of ~ exceeding 8. The depth of unconfined 


| water-bearing alluvial deposits as found in nature seldom exceed 500 ft. and 


: 


is ordinarily much less. Spreading basins of 500 to 1500 ft. width, and more, 


: 


b=ihe) 


are in common use. For a basin of only 500-ft. width the maximum value of = 


in actual practice would, therefore, seldom exceed unity. It is usually very 


much smaller. For a! = 0.1 the value of q is 3.6 sec. ft. for gradient unity, 


WwW 


as indicated by writer’s diagram 1. 
Referring to the model tests which the author has performed to supplement 


his mathematical analysis it is found that the value of the parameter 2 in the 


model was 12.8 ins. = 10.7.2 This extreme divergence in value of the param- 
1.2 ins. 


eter D between model and prototype leads one to question the propriety of 


WwW 
relating model and prototype as does the author (p. 806-4). Certainly in the 


second phase of mound growth, where the surface of the mound contacts water 
in the spreading basin and acquires the potential control of its surface, there 


is far greater probability of “water logging” of the surrounding ground for Ww 


< 1 that for a4 = 10. Where the spreading basin is adjacent to stream chan- 


nel there is also a greater probability of substantial return flow from spread- 
ing basin to the stream channel. 

Commenting on the latter possibility the author states that “no significant 
return flow from the spreading basin to the stream channel will materialize 
SO long as the gradient between the two is moderate,” (p. 806-8), suggesting 
a slope of 5 horizontal to 1 vertical as critical but not supporting this sugges - 
tion by theory or test. The only positive statement the author has made in 
this connection is that “to stop the flow from the spreading basin the water 


[. Cons. Engr., San Francisco, Calif. 


2. Baumann, Paul, Ground water movement controlled through spreading, 
Trans. A.S.C.E., 1952, Vol. 117, Figs. 18 to 21. 
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surface of the channel mound would have to rise to the same elevation as that 
in the spreading basin.” Author’s Fig. 3b is referred to in illustrating this 
statement. The writer has attempted a sketch development of the flow net in 
the region of an assumed river levee (writer’s Diag. 2). This is more de- 
tailed than is shown on the author’s Fig. 3 and is one which the writer con- 
siders to be more in accord with actual field experience as well as with re- 
sults of model tests for percolating flow under levees and dams. The 
percolating flow through the vertical section A-A', immediately back of the 
near river bank in a slice 1 ft. wide perpendicular to the river bank, can be 
computed from the equation Q=KAi, assuming K as 1000 Meinzer units, A = 


’ 
1 ft. x 92 ft. = 92 sq. ft. and i = So 


Q, = 1000 x (92 §_ _ 0.32 sec. ft. 


7.5x24x60x60~ 37 


A similar computation at the far bank of river in vertical section B-B' is 


1000 x 92 6' 
Qp = 7.5 x 24x 60x60 * 230° ~ °025 sec. ft. 

The indicated loss must be assumed to be seepage into the river channel, 

Qa - Qp = -286 sec. ft. 
eee 200 

This is =) 
which percentage would be greater for sections taken at lesser angle to the 
river. This computation illustrates that even for a gradient flatter than 5 to 
1, the seepage into the river may be “significant.” The author’s optimism 
regarding avoidance of significant return flow from spreading basin to chan- 
nel at gradients less than 5 to 1 would appear to be unfounded in the case of 
actual field conditions. 

Broadly viewing the author’s mathematical approach, as developed in his 
earlier paper, it may be pointed out that there is no law of physics but which 
is true only within certain limits. The author in developing his equations has 
omitted to investigate and define limits and in his former paper2 does not 
indicate within what ranges his physical phenomena are mathematically valid. 
For example, in setting up boundary conditions for unsteady flow he states in 
condition No. 2 that y = 0 at x = Lg for all values of t in which 0<Lg<oo. 

This disregards the case where seepage escapes from the bank of a ditch, 
river, or lake above free water surface, or the existence of escape gradient. 
It also assumes that the advancing margin of the mound has a feather edge 
rather than blunt nose. The rapid initial rise of water, as often observed in 
test holes following initial stream flow, would indicate a blunt nose. 

In the author’s condition No. 3 he states y =f (x) at t = 0°, apparently as- 
suming that the profile of the water table is a continuous linear function of 
distance from spreading ground throughout its length at any and all times. 
This assumption is predicated upon the transmitting material being homo- 
geneous for any value of x. This is a condition never met with in nature. It 
also omits consideration of discontinuity due to intervening horizontal water | 
surface in open ditch or river. Such condition is of common occurrence. 

The “fitting” of the “basic parabola” observed in model tests with the 


= 87% of the percolating flow at the near bank of the river, 


2. Op. cit., Equation 4, page 1027. 
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abolic type of equation developed mathematically for the mound profile for 
eady flow does not justify the extension of such equations beyond the specific 
mits of the model test. Obviously the results of such tests apply only to the 
mditions established by the model. In the absence of a great variety of con- 
itions, which only a great number of tests made on many different models 
ould provide, there exists no foundation for generalization or correlation 
tith mathematical equations. 
The author points out that his basic equation for unsteady flow? has its 
nalogue in the flow of heat through a prismatic non-radiating bar. It is to 
e noted that for this condition to obtain for percolating flow through granular 
il, the soil stratum must be enclosed in sealed surfaces which, in a flow 
et, would be represented by flow lines rather than by equipotential lines. In 
fgetice, this condition is fulfilled for a percolating stream in homogeneous 
confined granular by the free water table and the impervious stratum at the 
ottom. When such a flow is intercepted by an open ditch, river or lake, how- 
wer, the upper limiting surface changes from the water table to the stream 
nk and bed. The stream bank and bed then becomes an equipotential line, 
us changing the analogue to heat transference through a radiating prismatic 
- Hence, at this point a discontinuity occurs in the mathematical equations 
lating to the percolating stream. Recognition of this condition is an aid in 
derstanding why a substantial part of the percolating water seeps into the 
tream channel rather than moves “across the channel,” as stated by the 
uthor (p. 806-7). 


ROBERT T. KNAPP,* M. ASCE.—Mr. Baumann is to be congratulated on 
1is excellent paper on a very important subject. Although it was undoubtedly 
timulated by his interest in the water problems of Southern California, the 
oncepts and methods it presents are applicable to all semi-arid and arid 
egions where water supply is the limiting factor in determining the size of 
opulation that can be supported by the area. 

The writer was particularly interested in the Manhattan Beach experiments 
n the development of a fresh water barrier to prevent salt water intrusion in 
epleted groundwater basins. While following the action of the fresh water 
arrier as shown in the author’s Fig. 6a, b, and c, the writer observed a 
jinor discrepancy in Fig. 6c which seemed to have some interesting impli- 
ations. The author states that Fig. 6c was drawn with the aid of a crystal 
all. This figure represents the idealized ultimate condition when equilibrium 
2comes established. It seems to the writer that when the author was gazing 
rough the crystal ball at the ultimate position of the free surface of the 
ndward saline wedge, there must have been a minor irregularity in the 
rystal, because the free surface of the salt water appears to slope in the 
rong direction. : 

It seems reasonable to assume that in a density current flow in a porous 
edium, there will be a vanishingly small friction at the interface between 
e two liquids. In the first place, the relative velocity will be extremely 
nall. In the second place, the porous medium is fixed in space, and the 
ear that restricts the flow of the overlying layer is that between the liquid 

this layer and the porous medium. Hence, in the case of salt water in- 
usion, it can be assumed that, for equilibrium conditions, pressures will 


Op. cit., Equation 4, page 1027. 
Prof., Hydr. Engr., California Inst. of Technology, Pasadena, Calif. 
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always be constant along any horizontal plane and will be equal to the weight 
of a unit column of sea water whose height is the distance from that horizon- 
tal plane to sea level. The slope of the interface will therefore be deter- | 
mined solely by the pressure from the overlying layer. 

In the case of the landward saline wedge cut off from the ocean by the re- 
charging flow, pressures must also be constant along horizontal planes, but 
the pressure level will not be determined by the depth below the surface of 
the ocean, but solely by the pressure in the overlying fresh water. The con-. 
figuration of the interface between this isolated body of salt water when it 
finally reaches equilibrium will depend upon its geographical location at 
equilibrium. It will, of course, run downhill to the lowest available spot and 
there form a lake. In Fig. 6c it looks as if this lake was formed within the 
recharge zone, that is, where the fresh water mound is still sloping landward 
If this is the case, the interface would slope seaward, as shown in the writer’ 
Fig. 1. If this salt water lake should happen to come to rest under a region 
in which there were pumping wells, then the interface would include salt wa- 
ter mounds whose peaks would lie under the pumping wells, as shown in the 
writer’s Fig. 2. 

This principle of zero shear at the interface makes it very easy to calcu- 
late the possible distance of the salt water intrusion from the known ground- 
water elevations since the hydrostatic pressure at the interface must just 
equal the depth below sea level of the interface. The slope of the interface is 
always in the opposite direction to that of the groundwater surface. The in- 
terfacial slope is always the steeper, since the ratio of the two slopes is 
equal to the ratio of the density of the fresh water to the difference in density 
between the salt and fresh water. If the relative density of the salt water is 
assumed to have the nominal value of 1.025 with respect to the fresh water, 
the downward slope of the interface will be forty times the groundwater slope 
As stated previously, these considerations apply only to equilibrium condi- 
tions, and hence represent only the maximum salt water intrusion that would 
occur for a given set of continuing groundwater elevations. 
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Discussion of 
“EXTENDING STREAMFLOW DATA” 


by W. B. Langbein and C. H. Hardison 
(Proc. Paper 826) 


R. D. GooprIcu,! M. ASCE.—The paper on “Extending Streamflow Data” 
s one that the writer finds not only interesting but very timely and useful in 
onnection with the engineering investigations of the Upper Colorado River 
ommission. A great deal of its work is covered exactly by the title of this 
aper. The record of the authors wide experience in this field and their many 
uggestions for the proper application of mathematical and graphical correla- 
ion procedures should be especially useful wherever short term records 
just necessarily be utilized extensively as in the water supply studies for 
1is Commission. Not the least valuable feature of this paper is the authors 
iscussion of the limitations which are inherent in any method which can, at 
est, give only approximate results. 

One very important problem which has arisen in this Commission’s studies 
3 that of fitting the short term records at newly established stations at im- 
ortant points on a stream into the water supply picture at the earliest prac- 
cal date, and, at the same time, to know how reliable the estimated quantities 
ill be when the records at the new station are extended back for perhaps 20 
ears, 

On the Dolores River in Colorado, it has been found desirable to discon- 
nue the gaging station at Gateway because of changes in the stream channel 
sriously affecting both gage and control. This station is a few miles east of 
i¢ Colorado- Utah State line, and a new station has been established down- 
ream at about an equal distance west of the state boundary. After numerous 
ials and considerable study, it was found that when concurrent monthly re- 
rds in 1000 acre-foot units for three water years were plotted on log paper, 
the annual discharges were plotted with 10,000 acre-foot units, the points 
ll along the same general trend lines as the monthly records. Furthermore, 
-obably because of the fact that the two stations are on the same stream only 
} miles apart, the trend of the plotted points was nearly along a 45° line. It, 
erefore, appeared possible to compute an equation for the regression curve 
hich would express the relation of monthly discharges at these stations with 
Satisfactory degree of reliability. Could not the same identical equation be 
ed to give annual values of discharge simply by taking ten times the monthly 
it for the unit of annual flow? 

To test this possibility, the discharge records at two locations on the 
jlorado River about 75 miles apart were treated in the same manner. The 
onthly records for the water years 1951-2, 1952-3 and 1953-4 at Glenwood 
rings on the Colorado River were added to those on Roaring Fork for the 
tals at the upper station and the concurrent records at the station near 
1meo on the Colorado River were used for the record station for which 


Eng. Consultant, Upper Colorado River Comm., Grand Junction, Colo. 
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estimates of annual discharge were to be carried back for a total period of 21 
years, that is to 1934, when the Cameo station was started. This gave a 
series of 39 items including the three annual discharges. These records in- 
clude one year of high flow, one of very low flow and one of nearly average 
discharge. The standard error of estimate for this series, using discharges 
and not their logarithms, was 25.1 units of flow or 10.8%. In order to have a 
series that was as homogenous as possible, three of the monthly records were 
omitted for the next computation. One of these records had a departure of 
one standard error, one a departure of two and the third a departure of more 
than three standard errors. This first trial indicated beyond question that the 
use of the logarithms of the discharge would give a distribution that, as the 
authors state, would be more nearly normal than if the discharges themselves 
are used. The following equations were derived. Logarithms used are to 
base 10. 


log Y = 0.95989 log X + 0.13765 


¥ = 1.373 x°-960 


Based on the 33 monthly records and the 3 annual ones, the adjusted coeffi- 
cient of correlation was r = 0.9976 and the standard error was S = + 0.03161 
log units or + 7.55% and - 7.00%. When the actual annual estimates were 
computed and compared with the 21 years records available at the Cameo 
station, the results turned out better than indicated above. The mean absolute 
error was only 4.8% indicating a standard error for this period of + 6%. In 21 
years, 6 or 7 times one might expect differences of 6%or more. However in 
2 years, 1934, a low year, and in 1942, a high year, the error near 6% 
(-5.7%), and five times the difference between recorded and computed flows 
was less than 1%. More trials and study will be necessary before one may 
safely adopt this use of a single correlation equation for both monthly and 
annual discharges, but under similar conditions to those described the pro- 
cedure outlined seems to have sufficient promise to justify much further 
study. 

The writer has been very much interested in the study of this paper and 
has found much of value which can be used in his own work and little to criti- 
cize. Personally, whenever possible we prefer to use mathematical pro- 

cesses with graphical methods for illustration and checks. The accompany- 
ing diagram gives a comparison of the records with the computed discharges 
obtained in the above test. The authors are to be highly commended. The 
writers own less technical studies tend to confirm their conclusions. 


CARROLL F. MERRIAM! and E.'T. SCHULEEN2 M. ASCE.—Those who 
have had extensive experience in studying the hydrology of the Susquehanna 
River will approve the dominant thought expressed by the authors: that by 
means of extending short term records on the basis of correlation to long 
term base records, the knowledge of streamflow can be advanced in the most 
economical manner. The proposed use of roving satellite stations is the 
sensible solution to the critical demand for more and more basic hydrologic 
data. The engineering profession should encourage the U. S. Geological Sur- 
vey in furthering such a program. 


1. Prospect Harbor, Maine, formerly with Pennsylvania Water & Power Co. 
2. Assistant to Manager, Safe Harbor Water & Power Corp., Conestoga, Pa. 
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For a number of years the Pennsylvania Water & Power Company? and the 
Safe Harbor Water Power Corporation had been developing methods for analy- 
sis of past records of streamflow at gaging stations and hydroelectric power 
stations in the Susquehanna basin. It is important that one of these proce- 
dures is very similar to what the authors call “equal per cent duration,” and 
which they point out gives results, as far as can be shown in Fig. 2, identical 
with the correlation method. 

The principal differences in practice were in the use of daily instead of 
monthly data and the use of stage instead of the logarithm of discharge. The 
obvious objection to daily data, the tremendous amount of labor involved, was 
overcome by business machine methods. Punch cards, for which many other 
statistical applications have been found, were used to record the daily dis- 
charges at 21 stations for the past 24.years. Periods of time were selected 
for which records at the stations to be compared were continuous and un- 
affected by ice or other abnormalities. The sorting machine listed the dis- 
charges in parallel columns and in ascending order of magnitude. Inasmuch 
as there were the same number of entries in each column, those appearing on 
the same line had equal frequency of occurrence. In this phase discharge was 
used for the reason that the given data were so expressed, but conversion was 
made to stage on the lists by reference to the rating curve in effect at the time 
of observation. Plotting the stages shown on each line, one against the other, 
resulted in a relationship curve of equal probability of occurrence, similar to 
that shown at the extreme left of Fig. 2. This relationship could then be used 
to estimate stages and discharges at the station with a short term record by 
reference to those at the station with the longer record. Actually it is used 
to estimate the stage that is most likely to occur at a down stream station on 
the basis of an observed stage at an upstream station. 

Stage was used, not for the sake of normalizing the decided aiew in the 
discharge frequency, which the authors accomplished through use of a log 
scale of discharge, but rather because the work was of such a nature that 
shifts that had been made in rating curves from time to time would have 
caused substantial variations in the results. 

The authors seem somewhat surprised that the several methods give al- 
most identical results. Their conclusions appear to indicate that what the 
writers, in the development of stage relation curves, had termed “conjugate 
stages” or “conjugate discharges,” and which might be defined as stages or 
discharges having the same frequency of occurrence, do have greater reality 
than originally supposed. This concept should prove very useful in solving 
flood routing problems. 

If one were interested only in long term mean discharge, a much easier 
procedure than that proposed by the authors is the use of the double mass 
curve. The accumulated monthly inches of runoff for the station being in- 
vestigated are plotted against the accumulated runoff at a base station, or 
better yet, against the mass curve of the basin runoff pattern, which, in the 
case of the Susquehanna, consisted of the average of approximately 100 gages 
weighted by the square root of the drainage area. The slope of the double 
mass curve gives the relation of the mean runoff at the station to that of the 
basin. This method has been extensively used in studies of the Susquehanna 
to derive useful estimates from many records altogether too short or frag- 
mentary to be depended upon otherwise. In this way useful results have been 


3. Now merged into Pennsylvania Power & Light Company. 
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salvaged from about 50 gages since discontinued, many of them with only a 
few years of record. The reasonableness of these estimates has been veri- 
fied by comparing runoffs from gaged areas with the residual runoff from un- 
gaged areas on a “family tree” diagram.* The success of the methods lead 
the writers to conclude that the policy of roving satellite stations is sound 
and certainly should be adopted. 

It should be emphasized that this is not an end in itself, but rather the be- 
ginning of a long road toward greater achievements in the field of applied 
hydrology. Past experience has also shown that while the U. S. Geological 
Survey should be encouraged to forge ahead, private industry has much to 
gain by placing a shoulder to the wheel. 


WILLARD M. SNYDER,* A.M. ASCE.—A good, objective treatment of an 
example of the use of statistics in hydrology is presented. Certainly carefully 
selected drainage areas should be capable of sufficiently good correlation so 
that base stations could serve as “predictors” for satellite sites with rather 
closely controlled variances of prediction. The similarity in characteristics 
of adjacent streams has been used for such varied purposes as interpolation 
for missing records, derivation of synthetic unit hydrographs, and prepara- 
tion of regional maps of streamflow parameters. 

The authors invite discussion on particular points of their methodology 
when they state “------ the details of the techniques are far from completely 
developed.” There can be no serious question of the statement that statistics 
should be the tool of hydrology, rather than vice versa. This idea should be 
extended to state that all phases of methematical and numerical methods are 
but tools. The basic science must provide the physical logic. 

The graphical statistical correlations are particularly useful when the 
problem to be analyzed combines small samples with very evident correlation 
of a few variables. But when samples become large, and when causation is a 
combination of many separate parameters and iterative solutions required for 
graphic multiple correlation become agonizingly tedious. Under the latter 
conditions more economical analyses can probably be accomplished using a 
product-accumulating desk calculator and following the conventional least- 
Squares techniques for multiple correlation. It is acknowledged at this point, 
that techniques for selection of proper mathematical models for hydrologic 
analyses are also far from completely developed. 

It appears necessary to distinguish between statistical tests based on fre- 
quency distributions of uncorrelated data, and tests based on the residual er- 
rors in data after adjustment by regression analysis. Highly skewed fre- 
quency distributions must be “normalized” for testing, and a logarithmic 
transformation apparently accomplishes this purpose. But it is not self evi- 
dent that the frequency distributions must be normalized for regression analy- 
sis. It is more important** that the residual errors be normally distributed 
With equal variance along the line of regression. Least-squares fitting by 


4, Progress Report On The Analysis of Rainfall Data, Carroll F. Merriam, 
Transactions of the American Geophysical Union, 1938, p. 532. 

* Civ. Engr. (Hydr.) Tennessee Valley Authority, Hydraulic Data Branch, 

_ Hydrology Section, Knoxville, Tenn. 

** Anderson, R. L., and Bancroft, T. A. Statistical Theory in Research. 
McGraw-Hill, 1952. p. 153. 
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definition produces-minimum residual variance, no matter what the distribu- 
tion of the variates, and approximate statistical tests may be based on these 
results. The seriousness of the approximation is difficult to determine and 
in critical instances should be referred to mathematical statisticians. 

The authors use the expression “standard error of estimate” for the square 
root of the variance of residuals. This nomenclature is possibly misleading. 
A better expression is “standard deviation of residuals.” It is, of course, in- 
herent in the theory of least-squares that the residual variance be uniformly 
distributed along the line of regression. When fitting is done by logarithms, 
the fitting does produce a standard deviation of residuals which is a constant 
percentage of the dependent variable after transformation back to the original 
units. 

It does not seem correct to pass over the assumption that “the mean line 
is correctly located.” The analyst cannot know whether or not this is true. 
Consequently, he cannot know whether approximately two-thirds of future 
occurrences will be within a uniform band of plus and minus one standard 
deviation from the fitted line. When the regression is used for prediction of 
future events, the variance of the prediction is the residual variance plus the 
variance of the fitted line.. 

The variance for a single predicted value of Y' for a particular value X' is 


and the variance for the average of all predicted values of Y' for the particu- 
lar value X' is 


t —,2 
Sag) Dane tae 
= 2 
S (X-X) 
where s2 is the variance of the residuals 


X is the mean of X 


and §(X - X)2 is the sum of squares about X. 

From these formulas it can be seen that the variance of prediction depends 
upon the particular value of X' on which the prediction is based. Even for 
X' = X (prediction at the mean value) the variance of individual predictions 
is not the same as the variance of residuals, but depends upon the size of the 
sample. The variance for the average of all predictions does become s2/n 
for X' = X. The “engineering” preference for errors as percentages is thus 
not justified for variances of prediction but only for variances of residuals of 
fitting. The authors point out the break in variance between original record ~ 
and predictions. 

The distinction between the variance residual to fitting and the variance of 
prediction should be maintained for the “unregressed” model. 

The use of the geometric mean to obtain a unique line of relation between 
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two variables has been mentioned by other authors.! The slope of the fitted 
line by this method is the geometric mean of the two slopes obtained by mini- 
mizing the residual sum of squares along first the vertical and then the hori- 
zontal axes. Minimizing the product of the vertical and horizontal residuals 
produces the same result. It is not clear that the graphical method of obtain- 
ing the medians of vertical and horizontal strips would result in the same line. 
Any conclusions based on analogy between the two methods would have to be 
considered approximations if the lines are different, since the distribution of 
residual errors, and hence residual variance, will vary with the fitted line. 


Practical 
I. Gumbel, E. J. “Statistical Theory of Extreme Values and Some 
Applications. ” National Bureau of Standards Applied Mathematics Series 


No. 33 p. 16. 
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Discussion of 
“LUNAR-CYCLE MEASUREMENT OF ESTUARINE FLOWS” 


by Irvin M. Ingerson 
(Proc. Paper 836) 


FRANKLIN C. CRAIG,! A.M. ASCE.—This paper is an important contribu- 
tion to a subject that presents many complex problems. Much of the planning 
for the future development of the Central Valley area depends upon a knowl- 
edge of how much flow must be maintained in the lower Sacramento River 
channel. The measurement of flow discussed by Mr. Ingerson is therefore of 
considerable economic importance. 

Of prime interest is the high degree of accuracy that must be attained if 
results are to be truly indicative of net outflow. In a tidal cycle measure- 
ment accidental errors tend to become negligible because of the large number 
of observations that are made, assuming of course that exactly the same 
points of observation are maintained and other precautions are taken to as- 
Sure accurate work. Systematic errors which may apply unequally to flood 
and ebb tide flows offer the greatest possibility of producing erroneous re- 
sults. The assumed figures of Table 1, which are commensurate to actual 
volumes of flow, illustrate the importance of avoiding errors of this type. 
The table is based throughout on hypothetical flood flows of 200,000 cfs and 
ebb flows of 205,000 cfs. 


Craig.on Ingerson, Table 1 -- Comparative Error Possibilities 


Measured Assumed Measured Assumed Indicated Error 
Average Percent Average Percent Net in Net 
Flood Flow in Error Ebb Flow in Error Outflow Outflow 
cfs cfs cfs efs 
200,000 c@) 205,000 0 5,000 ie} 
204,900 +2 209, 100 +2 5,100 100 
220,000 +10 225,500 +10 5, 500 500 
204,000 +2 205,000 0 1,000 4,000 
204,000 +2 200, 900 -2 -3,100 8,100 
196,000 -2 209,100 +2 13,100 8,100 


Fairly large errors can be tolerated if they apply equally to flood and ebb 
flows but even very small errors can be disastrous if they are of varying 
magnitude and especially of opposing sign. 

Possible sources of systematic errors which might not apply equally to 
flood and ebb flows are (1) variable density effects caused by incomplete mix- 
ing of salt and fresh water, and (2) wind effects. 

The average of point velocities at .2 and .8 depth does not always repre- 
sent the true mean velocity in the vertical. This has been found to be especi- 
ally true in some tidal streams where the heavier salt water tends to move 


1. Hyd. Engr., Geological Survey, U. S. Dept. of the Interior, Sacramento, 
Calif. 
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along the bottom quite independently of the lighter fresh water on top. Flood 
and ebb tide flows are sometimes affected differently by variations in density. 
A limited number of velocity observations were made at each tenth of the 
depth, during the course of the measurement that has been described by Mr. 
Ingerson. These special observations suggest that the mean of the .2 and .8 
depth observations may be adequate at the measuring section that was used, 
at least under the conditions of tide, flow, and wind, at which the comparisons 
were made. Salt and fresh water mixing action is apparently more complete 
here than it is in some tidal streams. 

Some lateral shifting or wandering of currents in a direction not perpen- 
dicular to the measuring section occurs during periods of low velocities. 

Such errors as might be introduced from this source would have a tendency 
to be compensating during periods of flood and ebb flow. 

Little is known concerning the magnitude of error that may be introduced 
by wind. Abnormal movement of meter cups induced by vertical boat move- 
ment, distortion of vertical velocity relationships, shifting currents, and 
change in the normal mixing action of fresh and salt water are all possible 
effects of wind. Wind movement often varies; hence, errors of unequal mag- 
nitude may occur during flood and ebb flow periods. In the absence of definite 
knowledge concerning wind effects, it would seem to be good practice to place 
the most confidence in observations made during periods of minimum wind 
movement. If winds occur during a measurement, it would be good practice 
to extend the period of measurement well into a period of calm because of the 
effects of wind on channel storage. 

One of the purposes of the measurement of flow that has been described in 
the paper is to determine the amount of flow which may accrete to the stream 
from unmeasurable ground-water sources. The amount of these accretions 
may vary seasonally, or from year to year because of change in ground water 
conditions or consumptive use by crops. For this reason a continuous record 
of flow is desirable, and consideration should be given to possible methods of 
rating the channel. 

First consideration might be given to a proposal that the salinity gradient 
between two ends of a selected reach may be used as a factor in determining 
outflow. If the salinity gradient could be thus utilized, it would furnish the 
most economical method of obtaining continuous records of flow because it 
would eliminate the need for frequent expensive current meter measurements. 
Incomplete mixing of fresh and salt water would seem to present a major 
source of error in determining the salinity gradient. This would be especially 
true if the measurements of salinity bear a different relation to the true 
salinity during flood and ebb flow periods. 

The nature of a tidal stream is somewhat different from that of an ordinary 
stream of mild slope entering a lake. The ordinary stream has a water sur- 
face profile that is concave upward within the limits of backwater effect and a 
definite relationship exists between stage, slope and flow. In a tidal stream, 
however, the profile of the mean water surface tends to arch upward by an 
amount which varies with the range of tide levels, the volume of flow, the 
level of the ocean, and other factors. This variable arching wage? enters inte ; 
the problem of rating a tidal stream. 

The modern theory of tidal action in the ocean is that it consists of the 
rocking up and down of water in large natural basins across virtually tideless 
nodes. Tidal action is said to develop when the natural period of oscillation 
within the basin coincides with that of the tide-producing forces; principally 


ASCE DISCUSSION 955-43 


the moon and the sun. Tidal action at San Francisco (writer’s Fig. 1) is com- 
plex because the rocking action takes place along two axes. A tidal pattern 
results that is a mixture of the diurnal and semi-diurnal types of tides. 
These two basic tidal patterns are illustrated by tide graphs for Pensacola, 
Florida and Hampton Roads, Virginia. The mixed tidal pattern is character- 
ized by succeeding tides that frequently vary considerably in magnitude and 
that do not usually oscillate about mean sea level. These characteristics, of 
course, add to the complexities of streamflow measurement and rating prob- 
lems. 

When the ocean surface rocks upward, a wave of translation2 is said to be 
induced within the estuary which is propagated upstream until its energy be- 
comes depleted by friction and by the force of the downflowing current. Like- 
wise when the ocean rocks downward, a negative wave of translationlis induced 
which is also propagated upstream. However, the energy of the negative wave 
becomes depleted much more rapidly than that of the positive wave because 
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Craig on Ingerson, Fig. 1—Tidal Patterns 


the negative wave travels in shallower water and is therefore subject to 
greater frictional effects. Also, because more water must be moved under 
the influence of the negative wave, it has more work to do. It is the unequal 
depletion of the energy of the positive and negative waves that causes the 


1. Applied Hydrology, by Linsley, Kohler, and Paulhus; McGraw Hill Book 


Co., New York, 1949. 
2. “Flow of Water in Tidal Canals” Col. Earl I. Brown, Transactions 


American Society of Civil Engineers, 1932. 
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mean daily water surface profile of the tidal stream to tend to arch upward, 
thus creating a head to achieve a balance of wave force. 

The unequal rate of depletion of the energy of positive and negative waves 
creates rather special rating problems as is evidenced by the operation of a 
gaging station on the Sacramento River at Sacramento. Here the stream is 
subject to a range in tidal action of about two feet but not to a reversal in the 
direction of flow. A reach of channel extending twenty miles downstream 
from Sacramento to Snodgrass Slough has been used for rating purposes. In 
this reach the net effect of tidal action over a period of a complete tidal cycle 
is much the same as that which would.be created from any source producing 
variable backwater effects, except that the forces which create abnormal 
stages are active within the rating reach and may create even more abnor- 
mality in stage at the upstream end of the rating reach than they do at the 
downstream end of the reach. This is in marked contrast to the behavior of 
ordinary backwater effects. The result is that the fall through the reach is 
not a dependable measure of the effective slope at any point in the reach. 

This places a definite limitation on the type of rating relationship which may 
be used. 

Considerable success has been had in rating the Sacramento River at Sac- 
ramento by a rating fall method, a standard method used in rating many back- 
water-affected streams. This method makes use of a ratio between measured 
fall and rating fall, and a ratio between actual discharge and a rating dis- 
charge, in accounting for the effect of variable abnormalities in stage intro- 
duced by tidal action. Experience in the use of the method under a wide range 
of conditions is needed before its dependability for use on tide affected 
streams can be definitely established. 

Considerable success has also been had in rating the Sacramento River at 
Sacramento by methods which develop corrections either to abnormal fall or 
to abnormal stage. These methods have as their basis the depletion of wave 
energy as measured by the depletion in range of stage (between crests of 
positive and succeeding negative waves), from the downstream to the upstream 
end of the rating reach. The corrections to fall are applied to a unit-fall rat- 
ing extended downward as a straight line and the corrections to stage are 
applied to a base rating curve. 

The problem of rating the lower Sacramento River channel is complex be- 
cause of the physical characteristics of the channel and because the flow is in 
two directions. There are no reaches in the lower Sacramento River channel 
that are well adapted to the methods that have been used to rate the Sacra- 
mento River at Sacramento. With the exception of a few short, somewhat 
constricted reaches (author’s Fig. 1), the lower channel is very wide and is 
in the nature of a bay rather than a river. Since little fall develops during a 
portion of the tidal cycle, the minimum length reach that could be used if fall 
were a factor in the relationship, would be about fifteen miles unless special 
stage recording equipment were developed. In a reach of this length, a large 
portion of the flow would be continually entering or leaving channel storage, 
hence it is evident that special methods, which consider changes in channel 
storage, would be required to rate this channel. 

Also of concern is the effect of wind. Strong upstream winds prevail dur- 
ing much of the time, but light or strong, upstream or downstream winds may 
occur within a matter of hours. The fetch is about 30 miles and waves sever- 
al feet high often develop. Of more concern than surface waves are small 
changes in mean stage and rates and distribution of flow which are caused by 
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Craig on Ingerson, Fig. 2—Changes in Tide Pattern due 
to wind effects and normal depletion. 

wind. Here again the danger is that flood and ebb flow relationships will be 
affected differently and especially by winds of varying magnitude. The im- 
portance of very small changes in normal relationships is evident when it is 
considered that the entire range of flow, from zero to about 400,000 cfs takes 
place within a range of stage of only about 5 feet. A net outflow of 5,000 cfs 
would be represented by a change in stage of only about a tenth of a foot. 

Some idea of the magnitude of change in stage that is the result of wind can 
be gained from a study of the change in pattern of tidal action as it moves up- 
stream. Shown in the writer’s Fig. 2 are recorded tide crests at San Fran- 
cisco and gage-height graphs at Collinsville and Sacramento. Lines have been . 
drawn connecting wave crests in order to emphasize the pattern. At San 
Francisco the line connecting the crests of the negative waves makes a 
broad downward sweeping arc. At Collinsville this line reaches an unnatural 
looking peak on July 12 which appears to be associated with wind action. The 
most pronounced effects can be noticed when the wind quickly subsides, as on 
July 20. It appears that wind may affect normal river levels by as much as 
several tenths of a foot. Since wind effects depend on both duration and mag- 
nitude, they are difficult to evaluate. Unless the rating method properly con- 
siders the effects of wind, it would be of questionable value except under mild 
wind conditions. 

Since changes in upstream channel storage are continually in progress, 
flow in the lower Sacramento River must be computed through at least a fort- 
nightly and preferably a lunar monthly period in order to arrive at a depend- 
able value of net outflow. 

Quite probably, dependable records of daily flow might be obtained by the 
use of a number of continuously operated current meters suspended from pile 
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structures spaced at intervals across the stream. The velocities at the points 
at which the meters were suspended would be related to the average velocity 
of the stream as predetermined from regular measurements. Whether a few 
or many current meters and recorders would be required might be deter- 
mined from an analysis of the series of measurements made in September 
1954 and described in the paper. If a large number of meters were required, 
the system, of course, would be very costly. It should, however, produce 
records that would be dependable except during periods of very heavy winds. 

The limitations of the stationary meter in measuring extremely low velo- 
cities and the difficulty in determining direction of flow are recognized. A 
well balanced Price type meter in excellent mechanical conditon is depend- 
able in measurement of velocities of 0.3 foot per second. The moving boat 
method as described in the paper has shown that, on an average, the pattern 
of velocities in a tidal cycle, from a velocity of 0.3 foot per second through 
the period of slack-water to a velocity of 0.3 foot per second in the opposite 
direction, is very nearly a straight line. If the current meter were suspended 
from a rigid structure, there should be no appreciable error in adjusting ob- 
served velocities to this pattern during the relatively short periods when very 
low velocities prevail. This procedure would also eliminate need for concern 
regarding direction of flow. The use of a fixed meter has its place in the 
measurement of tidal flows if its limitations are borne in mind. 

The instrumentation research engineers are making rapid strides toward 
perfecting new techniques and equipment for measuring extremely low water 
velocities. One device being tested by the Geological Survey seems to offer 
great promise for measuring point-velocity by the rate of heat loss from a 
probe to the water flowing past it. This device is hoped to be free of some 
of the difficulties experienced with earlier versions of the hot-wire anemome- 
ter. Other developments of ultrasonic wave velocity measuring devices! may, 
in the future, permit recording the instantaneous mean water velocity on a 
line (sound wave path) across rivers. Such new devices may be of material 
assistance some day soon in measuring complex tide-affected streams like 
the Sacramento River. 

The moving boat method and the special methods of analysis of data that 
have been presented in Mr. Ingerson’s paper are ingenious solutions to prob- 
lems not often encountered in ordinary stream gaging operations. These 
methods may prove to be of value in other special situations such as the mea- 
surement of streams of very low velocities and especially in the measure- 
ment of flow of other tidal streams. 

In view of the new data obtained in the measurement described in the paper, . 
it is believed that additional observations and studies in line with the avenues 
of approach presented in the foregoing discussion might point the way to the 
development of a reliable rating of the net outflow from the Sacramento-San 
Joaquin Delta. 


L. J. TISON.2—A first very interesting method presented in this paper is 
the new “moving boat method” of using a standard current-meter to measure 
the low velocities. 

A second remarkable part is the description of the “Delta Outflow 


1. Swengel, R. C., Hess, W. B., Waldorf, S. K., “The Ultrasonic Measurement 
of Hydraulic Turbine Discharge”; A.S.M.E. paper No. 54 A 54. 
2.’ Prof., Univ. of Ghent, Belgium. 
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Measurement” during a lunar-cycle, between times of equal tidal-volume. 
By choosing-times of equal tidal-volume, said the author, the item of tidal- 
volume goes out of the hydrologic equation and this equation is reduced to the 
terms of accurately measured flows and to only one unknown: the net use of 
water in the Delta area. 

A theory presented in 1908 by a Belgian engineer Van Brabandt, suggests 
another method for the determination of the net use of water in the Delta. 


Theorem of Van Brabandt 


The volume between two surfaces, the first one joining the heights of the 
slack-water reversal on the ebb in the different cross sections upstream the 
cross section where the flood flow is to be established, the second one joining 
the heights of the following slack-water reversal on the flow, is equal the 
flood flow modified by the volume of the lateral flow. 

The demonstration presented by Van Brabandt was a very long one and 
was based on geometrical considerations. We present the following demon- 
stration. 


AB is the bottom of the river and CD is a simultaneous tidal line between two 
cross sections @ en w : @ is the fixed last cross section of the river, while 


the cross section w varies with 
b+5e 
The tidal-volume change — se equal the algebraic summation of the 


At 
volume flowing laterally « and of the volume flowing through the cross 
brat bear 


Frat br Ok beat trdr 
Niet eae UGE Tes ) a (1) 
b ¢ [a t 


hence: 
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ov) TA, 20 Dr _, 0 
Frag 1 te” OOS BEDS 
The volume flowing through the fixed cross section & is not altered by the 
variation of “ and A equal zero. 
Hence equation (2) reduces to: 
ON, po, + eel = 0 (3) 


E045 Bt OeO3 


On the other side the flood flow to be determined in an arbitrary cross sec- 
tion “4 is: 


e 
Fy = q, ab (4) 
¥ 
te is the time of the slack-water reversal on the ebb and t; is the time of the 
following slack-water reversal on the flood, in the cross-section w ; q, is the 


instantaneous flow in the same cross-section. 
We obtain, taking the derivative of (4): 


G bt 
le BS 2 i dt « | dk + G9, 4), (5) 
t 


The two last terms are the values of the instantaneous flows in the cross sec- 
tion w for the times of the slack-water reversal on the flow and of the slack- 
water reversal on the ebb respectively. These values therefore equal zero 
and we write: 


b. t, rk. 
M2 ae | bea Fooln dt (6) 


while 


Var : 


Owing to the equation (3), we may write: 
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age + | SAP - [tia o 


High Water 
4s of s/ack- water 


Heid 


Fig. 2 shows different characteristic curves for a tidal river. The cross sec- 
tion L is the first cross section presenting slack-water reversals. Upstream 
this cross section, there is no more flood flow. 

Integrating the equation (8) between the cross-sections a en L , we obtain: 


+) + (3 [ea f'n [Be 4 


the flood flow indeed equal zero in the cross section L. 


oi 
Yd 
a cau is the tidal-volume change, between the cross sections 


s ens + ds, from the time of the slack-water reversal on the ebb in the sec- 
tion s and the time of the slack-water reversal on the flood in the same sec- 


tion. 
It is therefore represented by the hatched part of the fig. 3 where the 
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curves (1) and (2) are the simultaneous tidal lines for tj and te respectively. 


A 
spa v 
F20b is consequently the volume between two surfaces, the firs 

by 


one joining the heights of the slack-water reversals on the ebb, the second 
one joining the heights of the slack-water reversals on the flood. 


(s 
Se F, i) 
It is easy to see that [o a4 t dt is the lateral flow. In the problem 
C] by 

of the Delta at the confluence of the Sacramento River and the San Joaquin 
River, this lateral flow is the part of the net use in the Delta area during a 
flood period and between the sea and the first cross-section L presenting 
slack-water reversals. 

This part of the net use in the Delta is therefore the difference between 
flood flow and the volume of Van Brabandt. 

The result of this method cannot be compared with the result of the author: 
this last result is an average for a total tidal-cycle and for the whole Delta. 

However we feel that our method presents some interest. 
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Discussion of 
“INTEGRATING THE EQUATION OF GRADUALLY VARIED FLOW” 


by Ven Te Chow 
(Proc. Paper 838) 


CLINT J. KEIFER,! A.M. ASCE and HENRY HSIEN CHU,2 A.M. ASCE.— 
Integrating the equation of gradually varied flow has been attempted by many 
investigators and success has been limited only to certain approximations or 
for certain specific type of channels. The author is to be complimented for 
advancing the methodology in the fact that only one set of tables need to be 
used to solve both the friction loss and velocity head change terms. However, 
the writers experience some disappointment in trying the author’s method for 
the closed conduit case. When the discharge through a closed conduit is 
greater by a few percent than K/So for the full section, which should not be 
considered as an infrequent occurrence, the normal depth y, becomes non- 
existent and rendering the final formula, Equation (38) meaningless. Even 
within the realm of its applicability, the writers attempted to test the steadi- 
ness of exponents N and J which should remain reasonably unchanged for any 
considered range of proportional depth y/D. Figure 6 shows a plot of N, M, 
and J versus y/D for the circular section. (The two broken curves, m and 
N/m will be explained later). The N and M curve was already presented by 
the author in Figures 3 and 4, but for a smaller range. The writers are 
rather surprised at the rapid variation of the exponent J for any y/D beyond 
0.4. The J value approaches positive or negative infinity as the proportional 
depth nears 0.65 from below or above respectively. 

Lee and others(17) presented a very similar method in which they arrived 
at the following equation, using the author’s notations 


Y Vn 7 pm N 41 
X=- 2 [u-F(u, N) + x5, F(u" I + CONSTANT (41) 

where m is another hydraulic exponent defined by 
Ameen (42) 


which is similar to Equation (37) except in the third term, the varied flow 
function, N/m has been chosen as the exponent instead of J. It can be shown 
that 


m= aye (43) 


[. Senior Sewer Designing Engr., Bureau of Eng., Dept. of Public Works, 


Chicago, Ml. ’ 
2. Civ. Engr. IV, Sewer Planning Div., Bureau of Eng., Dept. of Public Works, 


Chicago, Ill. 
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Fig.6. CURVES OF HYDRAULIC EXPONENTS 


FOR CIRCULAR SECTIONS 
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A plot of m and N/m versus y/D are also shown in Figure 6 as dotted curves. 
It seems to the writers the N/m curve possesses better steadiness than the J 
curve and Lee’s method should be good for y/D up to 0.8. However, it still 
occurs to the writers that the closed conduits actually present a different 
problem which has not been solved satisfactorily for the region near crown 
by any method presuming constant hydraulic exponents. The writers pre- 
sented a different approach(19) devised for closed conduit primarily, but the 
principles hold true for open channels also. 

As far as open channel cases are considered, the writers wish to see a 
graph showing the variation of J with respect to y/b using z as parameter. 
That will provide a better guide to those in the profession who intend to adopt 
the author’s method and want to know where they stand on the matter of 
accuracy. 


REFERENCES 


17. “Gradually Varied Flow in Uniform Channel on Mild Slope” Ming Lee, 
Harold E. Babbitt and E. Robert Baumann, University of Illinois Engineer- 
ing Experiment Station Bulletin Series No. 404, Vol. 50, No. 28, Nov. 1952. 


19. “Backwater Function by Numerical Integration” Clint J. Keifer and Henry 
Hsien Chu, American Society of Civil Engineers, Transactions 1955. 


CORRECTIONS 
u 7% a d Waa ” 
P. 2 - Equation (3) Denominator at right side / +4 (4) should be /*+ ay! 29) 
“" » 


P, 3 - Equation (13) Right side ° & "should be - oa 


s 


P, 9 - Line 21 from top “From Figure 3” should be “From Figure 5.” 


ROBERT Y. D. cuun,! J.M. ASCE.—Dr. Ven Te Chow has advanced the 
analysis of integrating the gradually varied flow equation for prismatic chan- 
nels a step forward toward simplicity and accuracy. But most important to a 
practicing engineer, Dr. Chow has presented a new direct method of analysis 
through the use of simple charts and one table. The method is accurate 
enough to be applied to practical problems involving channel shapes most 
frequently encountered, and can be used rapidly enough to be time saving. 

Once a graduate student of Dr. Chow, the writer fully appreciates the 
masterful manner in which this broad subject is presented in a concise, or- 
derly, and simple form. However, the writer would like to correct an error 
found in the general differential equation of gradually varied flow (Equation 
3), and to discuss the nature of the profile of the water surface by utilizing 
the equation for gradually varied flow (Equation 14). Equations 3 and 14, 
from Proceedings Paper 838, are repeated here for convenience. 


Oye as eee (3) 
ay _ - 


[. Asst. Hydr. Engr., Div. of Water Resources, Los Angeles, Calif. 
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rs 

£95 98> yooe He (14) 
dx 1-26 
ee 


In Equation 12, Dr. Chow reveals that either the second term in the 
denominator of Equation 3, or Equation 12 may be incorrect. The second 
term in the denominator of Equation 3 is incorrect, and will be corrected 
herein. 


FIG 1 GRADUALLY VARIED FLOW IN A VERY SHORT REACH dx 


Figure 1, as shown on page 838-28, is reproduced above with a few minor 
changes. Using Bernoulli’s Theorem, the upstream and the downstream sec- 
tions, as shown in Figure 1, were equated. After canceling, rearranging and 
factoring, the Equation 3-Corrected is obtained. 


Ay Se etic 4B x ¥y 
2g dy 


dy _ So - 5S 
(be o> ue (Ee) 
dy 


Nature of Profile 


With only minor computation, the equation of gradually varied flow (Equa- 
tion 14) can be a useful tool in determining the nature of the profile of the 
water surface of a gradually varied flow condition. It is a frequent occur- 
rence that the engineer would desire only the direction of the slope of the 
water surface. Also, if the engineer knows the nature of the profile, he can 
approach the solution, with the method outlined by Dr. Chow, more intelli- 


gently. Two general conditions of flow, subcritical and supercritical, will be 
discussed hereinafter. 
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Subcritical Flow 


When subcritical flow occurs, the following relationship exists: 


Y>Yec 
ZZe 


However, the depth of flow can either be greater or less than the normal 
depth. When the depth of flow is greater than the normal depth: 


Y>Yn,thenK>K, 


Substituting in Equation 14: dy 7) = a negative sign 
dx (+) 
The negative value shows that the slope of the water surface is rising to the 
left, assuming the usual sign convention, and the direction of flow to be in the 
same direction as shown in Figure 1. 
When the depth of flow is less than the normal depth: 


y< Yn » thenK< K, 


Substituting in Equation 14: dy _  (-) _ positive sign 

dx (+) 
The positive sign shows that the slope of the water surface is lowering to the 
eft. 


Supercritical Flow 


In the case of supercritical flow, the same procedure can be carried out as 
n the case of subcritical flow, and the following relationship can be shown: 


With Supercritical flow: y<Y¥e 
Z< Ze 
With depth of flow less than normal depth: y < yn, then K< Ky 
dy (-) viescncea 
— =--+— = a negative sign 
dx (-) : 


.’.The slope of water surface is rising to the left 
Nith depth of flow greater than normal depth: 
y>yYn,thenk > ky 


dx 


.’.The slope of water surface is lowering to the 
left 


dy =- a =a positive sign 
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: 
Shown below is Table A giving a summary of the discussion on the nature of | 
profile. | 

TABLE A 


SUMMARY TABLE ON NATURE OF PROFILE, 
UTILIZING THE BQUATION, 


1 - KG 
= 
eae» PV mak ae 
dx 7 1 - Ze 
zn 
Subcritical Flow : Supercritical Flow 
YP>V¥e3 F7& 3 Y<¥e3 22 ZB 
Y > 3 : YaVn ; yo oon y<Yn 
K 7 Ky ee Ky So K+ ky 2 K<ky 
aye yer io. way eyes _ -e 
ax = cigs GY Asa = ue eae at 
Slope rising; Slope lower-; Slope lower , Slope rising 
to left . ing toleft , ing to left , to left 


The Division of Water Resources, State of California, under the direction 
of Harvey O. Banks, State Engineer, and Max Bookman, Principal Hydraulic 
Engineer, Engineer-in-Charge of the Division’s Southern California office, 
encourages its engineers to participate in preparation of technical papers 
and in its discussions. This environment is fully appreciated. Also, the 
assistance of Jack J. Coe, Associate Hydraulic Engineer and Robert F. 
Clawson, Associate Hydraulic Engineer, both of the Division of Water Re- 
sources, for reviewing and commenting on this discussion is acknowledged 
with thanks. 


MASASHI HOM-MA.!—The integration method of the equation of gradually 
varied flow has been improved step by step by many authors. Now the prob- 
lem seems to be almost perfected by Mr. Chow. Although the adaptability of 
the Manning formula to the flow in steep channels still remains a problem fo 
future investigations, the other assumptions adopted by the author seem to be 
reasonable and have not much to be added. 

Regarding the author’s method, a point of theoretical interest may be 
added. The author’s differential equation of gradually varied flow may be 
written as: 


‘dx == 3n. (1-92 hve Mou al (1) 
du So 1-uN yn 1- uN 7 


1. Prof., Faculty of Eng., Tokyo Univ., Bunkyo-ku, Tokyo, Japan. 
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Mouret2 has pointed out that the critical slope S, is a function of water depth 
and this fact causes the water surface profile to show a characteristic of flow 
in steep channels or even in channels of likely mild slope (Fig. 1). This char- 


acteristic is that Eq. 1 will produce a point of inflection on the surface profile 
in the domain of y> yy >yc. 


Fig.l Movrefs charaderistic for the flowin wide 
rectangular Channel 


If Nis equal to M in Eq. 1, we will get S. = constant and such phenomenon 
cannot be explained analytically. The slope S¢ is defined as the value of Sy 
which is obtained by eliminating Q from the following equations: 


1 - [yp (@, 8) / yy] “=0 (2) 


= 1-[ye @D/y]™=0 ®) 


In these two equations Q is coupled with y as Q2/yN and Q2/yM respectively. 
Therefore, if Nis equal to M, y must be eliminated with Q. Then Sc becomes 
independent of y and in this case the inflection point, if any, can not be de- 
termined. 

‘Now if Nis not equal to M, there may be a possibility of explaining 
Mouret’s characteristic by the aid of numerical calculation. The condition 
that Eq. 1 has an inflection point in the domain of y > yn > Yc, 28 shown in 
Fig. 1, requires the following equation: 


(N - M) +MuN = n( 22M at the point of inflection (4) 
M Yc 
wu 


2. Mouret, “Hydraulique Generale,” 1927. The writer has not seen the 
original book but the outline of corrections by Mouret and other authors 
as presented by G. Formica in “Energia Elettrica,” No. 3, 1955, p. 249. 
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The writer has tried.an approximate calculation for the example presented by 
the author, where N = 3.65, M = 3.43 and yp/yc = 3.36/2.16. The result indi- 
cates that the Mouret characteristic occurs atu = .185. This condition of 
flow is doubtlessly out of the question for practical considerations. However, 
the writer is interested to see whether Mouret’s characteristic has any sig- 
nificance to the flow in a channel of circular section. This problem perhaps 
can be made clear by a numerical calculation. 


ALLAN NEWMAN, ! J.M. ASCE.—Since the publication in 1932 of the 
Bakhmeteff(13) method for integrating the equation of gradually varied flow, 
there have been various modifications of his procedure in an attempt to re- 
duce the labor involved in its application. The Engineer in attempting to 
evaluate these various methods is naturally led toward choosing one which 
has a sound theoretical basis, conforms to the physical evidences of the phe- 
nomenon it describes, and which is applied with a minimum of time consum- 
ing computation. It seems that on this basis the Bakhmeteff method supple- 
mented by the graphs of Kirpich, (2) must be chosen. Its advantages over the 
Von Seggern(18) method have been previously discussed by Kirpich and will 
not be repeated here. The author’s work has succeeded in eliminating the 
additional table required by the Von Seggern method, but at the expense of 
introducing certain computational complications. A comparison of the 
Bakhmeteff equation for gradually varied flow with that of the author illus- 
trates this. 

The Bahkmeteff equation is: 


L= 244-4) - (8) FS, N)-F4, NT 


where 0 (sigma) the critical slope and N the hydraulic exponent are obtained 


from the Kirpich graphs by computing 4 and (4486) bs 


The author’s equation is: 
L= 2 )(4-4)-[FU,N)- FY, ND] 
+ el") LF Gi) - Fu) )} 


where M, the exponent for the section factor, and N the hydraulic exponent is 


obtained graphically by computing + , J equals N and the value of Ye 


is obtained graphically by computing rasa 5 
It is clear that the complicated third term of the author’s equation may be . 
replaced by multiplying the second term by the simple factor (\- 22) - In fact 


1. Project Engr., Andrews, Clark and Buckley, Cons. Engrs., New York, 
N. Y. 
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z 
by noting that, @= as , it can be shown, by suitable arithmetic procedures, 


that the simple Bakhmeteff equation forms an identity with the author’s 
equation. 

While in some applications values of o must be obtained by steps, these 
would be no more numerous over a given reach, than those required for the 
evaluation of the author’s exponent M, inasmuch as an error in this factor is 
compounded by being in effect multiplied by itself since it also controls the 


M 
values of ts ; J and fF ( uv J) . Also, the value of N in the 
Bakhmeteff equation for practical purposes need only be evaluated to the 
nearest tenth which is more consistent with graphical accuracy, while the 
author abstracts the values of N, and M, to hundredths from his graphs, a 
procedure which is probably required by his method. Finally, while some of 
these difficulties are apparent in the author’s comparatively simple illustra- 
tive problem, solutions by the author’s method of the somewhat more com- 
plex problem of computing a discharge curve (Von Seggern(18) ex. 7) or of 
the problem of obtaining water surface elevations at specific distances from 
the control depth, would show how onerous the computations become, as com- 
pared with solution by the Bakhmeteff equation. 

Two aspects of the paper deserve commendation, 

First, the Bakhmeteff table of varied flow functions is now published as a 
Proceedings Separate making it more generally available, and secondly a 
brief survey of the various existing methods of computing water surface 
curves has been included, giving the casual reader a broader perspective into 
the subject background without requiring the investigation of numerous 
references. 


STEPONAS KOLUPAILA.!— Author is to be acclaimed for his ingenious ef- 
fort to find a correct and lucid solution for a complicate problem of nonuni- 
form flow in prismatic channels. An excellent idea to confine the computation 
to a single function was crowned with success. This varied-flow function was 
originated by B. A. Bakhmeteff; tables provided for these functions proved to 
be too short, thus the author put enormous amount of work in expansion of 
these tables. This is his particular achievement. 

The use of a single table is an answer to the problem, as suggested by the 
writer in his discussion to the paper by M. E. Von Seggern (Transactions 
ASCE, 115(1950), p. 95). 

It is regrettable, however, that author deliberately disregarded the fact of 
nonuniform velocity distribution across a current and neglected so-called 
Coriolis’ correction factor. So did Bakhmeteff, although he had agreed to 
change this in a new edition. 

Writer believes that this new excellent method will find wide appreciation 
among hydraulic engineers. 

The summary of existing methods (table 1) could be enlarged by several 
more, as by H. Walther (1903), N. N. Pavlovskii (1924), N. Westerberg (1924), 
I. I. Levi (1928), A. N. Rachmanov (1930), R. R. Chugaiev (1931), M. D. 
Chertousov (1934), P. W. Werner (1940), A. Vitols (1942), I. I. Agroskin 


1. Prof. of Civ. Eng., Univ. of Notre Dame, Notre Dame, Ind. 


955-60 HY 2 April, 1956 


(1944), R. B. Jansen (1951), T. G. Voinich (1953), etc. 

Writer would like to improve one historical detail. The tables for use of 
Bakhmeteff method were not lost in the turmoil of revolution and civil war: 
they were copied in hand-written form, until they had been added to a second 
Russian edition of the book in 1928. These tables had no typographic errors; 
exponent were different (2.0, 2.5, 3.0, 3.25, 3.5, 3.75, 4.0, 4.5, 5.0, and 5.5) 
and the values of 7 were more detailed. These tables had been reprinted 
many times in Russian textbooks and manuals (not mentioning the name of the 
author anymore; such plagiarism can be explained by amnesia). 


CORRECTIONS. — 


p. 838-7, line 3 from the top, “Let v = u2” should be “Let v = uN/J” 


p. 838-7, line 1 below Eq. (37), “in which u = y/y,, v = uJ” should be “in 
which u = y/yp, v = uN/J.” 


p. 838-10, line 8 from the top, “and v = uJ” should be “and v = uN/J_» 


p. 838-10, the 3rd and 5th columns of v and F(v,J) in the table under step 
4) of the example should be corrected as 


Vv F(v,J) 
1.625 0.213 
1.015 1.293 

-1.080 


p. 838-10, the equation under step 6) of the example should be corrected 
as 


L = 2,100 [0.476 - (-0.877) + (0.180) (-1.080) ] = 2,430 ft. 
p. 838-13, line 11 from the bottom, “v = uJ” should be “v = uN/J,” 
p. 838-25, line 1 from the top, “(p. 10 or 12)” should be “(p. 10 of 12).” 
p. 838-27, line 1 from the top, “(p. 12 or 12)” should be “(p. 12 of 12).” 
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Discussion of 
“RAINFALL DEPTH- DURATION RELATIONSHIPS” 


by Herbert M. Corn 
(Proc. Paper 840) 


D. M. HERSHFIELD,! and W. T. WILSON,2 A.M. ASCE.—The author’s 
problem evidently was to synthesize the rainfall-intensity regime for a re- 
mote region lacking short-duration data. In the absence of local data, United 
States data were used (Yarnell’s tabulation of each station’s most intense 
rainstorm as of 1933). 

A question that immediately arises is: how analogous is the rainfall re- 
gime of the remote area to that of the area from which the author selected 
his 60 stations? The United States storms studied were mostly summer 
storms and, as the author states, were generally convective. A convenient 
measure of the degree to which convective storms typify a region’s climate 
is the incidence of thunderstorms. In the central United States, the region 
from which the 60 storms apparently were selected, the mean annual number 
of thunderstorm days varies from about 30 to 70, with an average of about 
50.° Examination of climatic records for French Morocco shows the mean 
annual number of thunderstorm days to vary from about 2 to 20, with an aver- 
age of about 8. 

The moisture source region for intense storms in the central United States 
is the South Atlantic and Gulf of Mexico, where surface dewpoints approach- 
ing 80°F ensure a high moisture charge. In French Morocco, as in California, 
the moisture source region is the relatively cold water to the west of the 
coast, where dewpoints rarely exceed 60°F—a much lower moisture charge. 

The 2-year 24-hour precipitation in the central United States varies from 
about 2 to 5 inches, whereas in French Morocco less than 10 percent of the 
stations report a 2-year 24-hour value in excess of 2 inches. Figure 7, which 
shows the relation of hourly to 24-hour rainfall as.a function of thunderstorm 
2xperience, is based on northern United States stations for estimating short- 
juration rainfall in remote cold regions. 

Another question that deserves consideration is: what basis is there for 
xxtending the generalized depth-duration curve beyond the range of the data— 
say, from 80 minutes to 12 or 24 hours? The curve shown in figures 3 and 4 
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is based on short-duration rainfall; in spite of the ‘Introduction’s’ implication 
of its application to estimating short-duration rainfall from 12-hourly data 
from remote areas, and in spite of the ‘indications’ suggested in the ‘Conclu- 
sion,’ the writers question any basis for extrapolating the curve to durations 
of 12- or more hours. 

The mass curve of rainfall for most days of appreciable rain, or even for 
most storms, is typified by periods of relatively intense rainfall, with little 
or no rainfall the rest of the time. These periods of intense rainfall are 
known to meteorologists as “burst,” and as Corn has shown by his analysis of 
Yarnell’s data they tend to have a typical and continuous distribution of rain- 
fall increments. The writers see no reason to assume that the “within burst” 
distribution would be anything like the “among burst” distribution. Among the 
bursts, according to analysis of mass curves of daily rainfall of, say 2 inches 
or more, the periods of no rainfall typically occupy nearly half the day. Ex- 
ceptions to this include hurricane rainfall. Attempts to generalize on the 
typical distribution would have to recognize seasonal, regional and other dif- 
ferences. 

The curve of figures 3 and 4 says, as we understand it, that the ratio of 5- 
to 120-minute rainfall is the same as the ratio of the 1- to 24-hour rainfall, 
and that this ratio has general application. The ratio of the 1- to 24-hour 
rainfall from this curve is about 10 percent, whereas this ratio varies from 
about 30 percent to 60 percent in the central United States and about 10 per- 
cent to 70 percent for the entire United States. Figure 8 shows a regional 
difference in the 10- to 60-minute rainfall between coastal and interior zones 
of western United States.° This regional difference is related to average 
hourly rainfall amounts. Figure 9° is based on Yarnell’s and other data and 
shows that the ratio of 10- to 60-minute rainfall averages about 0.55 for total 
hourly rainfall of 1.0 inch to 2.0 inches, and about 0.30 for total hourly rain- 
fall of 4.0 inches to 5.0 inches. 

To illustrate the limitation of extrapolating such a curve beyond the range 
of data, we can use the Lexington, Kentucky storm of July 3, 1931, which is 
described in figure 5 and table 2. The suggested method has been applied to 
a total duration of 50 minutes, and is shown to work well. It fails when ap- 
plied to the 24-hour period, which has a total depth of 3.57 inches for 24 
hours. An hour is roughly 4 percent of a day; entering figure 4 with a dura- 
tion of 4 percent, the corresponding rainfall is about 10 percent of the daily 
depth. Ten percent of the daily total, 3.57 inches, is .36 inch, whereas the 
observed value is about 8 times that much. This difficulty can be avoided in 
part if the duration of actual rainfall for the reporting period of 24 hours (or 
12 hours) is known. This information is seldom available, however, except 
at stations having continuous observation of the weather. Also, the sequence 


5. U. S. Weather Bureau, Dept. of Commerce. Weather Bureau Technical 
Paper No. 24. . Rainfall intensities for local drainage design in the United 
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periods. Part I: West of 115th meridian. Washington, D. C. 1953, 

Figure 6. 
6. U. S. Weather Bureau, Dept. of Commerce, and Corps of Engineers, War 


Department; Hydrometeorological Report No. 5, Thunderstorm Rainfall, 


Vicksburg, Mississippi, 1947, Figure 141. 
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of rainfall increments is rarely available at stations where a synthesis of the 
regime would be pertinent, 

Another feature of this curve should be noted. The coordinates of figure 2 
are Cumulative precipitation and cumulative time, both expressed in percent 
of total; the points are distributed quite closely around a simple curve on 
semi-logarithmic paper. This good fit will also be the case when two series 
of random numbers, with the same operations performed on them, are plotted 
on either logarithmic or semi-logarithmic paper. 

The author briefly discusses the lack of similarity in the measurements 
for stations of rather close proximity and concludes that these differences 
are a result of poor observational practices. Convective precipitation is 
naturally spotty and it is not unusual for the heaviest precipitation to be con- 
centrated in a small area. Whatever dissimilarity in gage measurements 
appears may well be a function of convective rainfall pattern rather than an 
indication of the quality of the observations. 

An examination of the Burlington, Vermont, maximum annual rainfall 
values for durations of 5 minutes to 24 hours for 50 years of record reveals 
that it is not unusual for 75 percent of the rainfall in these storms to occur 
in about 10 percent of the storm time. What the author refers to as 
“peculiar” is simply a characteristic of rainfall regime for Burlington. In 
fact, this situation occurs, on the average, about once every three years. 

In conclusion, the writers heartily endorse Mr. Corn’s statement that more 
research is needed, and his hope that continued cooperation of hydrologists 
and meteorologists will produce the desired answers. 
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GEOLOGY OF SOME AMERICAN ESTUARINE HARBORS 
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ABSTRACT 


The essential requirements of harbors are protection from wave action, 
space for maneuvering of ships, and adequate depth of water. The principal 
causes of natural bodies of water suitable for harbors are (1) rise of sea 
level which causes the sea to inundate previously-formed valleys or low 
places; (2) deformation of the land, causing the land to sink so that the sea 
occupies the depressions; (3) differential erosion, which causes the soft rocks 
to be washed away, leaving low areas where the streams flow and valleys 
form, so that when the sea level rises or the land sinks, water can come into 
the low-lying areas; (4) near-shore deposition of bars which form barriers 
to bays, thus protecting them from waves; (5) glacial action which either 
scours out a basin which subsequently fills with water or deposits morainal 
and other material, forming islands or bars that protect the harbor from 
waves; and (6) mouths of large and deep rivers. 

Examples of typical harbors developed by these processes are Boston, 
formed by glacial action and depression of the land; New York, Philadelphia, 
Norfolk, Savannah, Charleston, and San Diego, formed by drowning of rivers 
sither by rise of sea level or lowering of the land; New Orleans and Portland, 
lear the mouths of large rivers; Galveston, by the formation of barrier bars 
und widening of the sides by wind action; and San Francisco and Seattle by 
leformation and lowering of the earth’s crust. All these harbors except New 
yrleans and Portland show the effects of differential erosion. 

Factors that affect silting and shoaling of harbors are: distribution of 
surrents; quantity, size-distribution, and mineral composition of sediment; 
‘mount and rate at which water is brought down by streams; salinity distri- 
ition; type of clay minerals; and amount of organic matter and sewage in 
he water. 


jote: Discussion open until September 1, 1956. Paper 956 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
82, No. HY 2, April, 1956. 
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INTRODUCTION 


The control of harbors is a difficult problem. Two essential requirements 
of harbors are: (1) location in a place protected from wave action, and (2) ade- 
quate depth of water. The distribution of depth of water in natural harbors 
represents equilibrium conditions with the environmental factors affecting the 
area. If man modifies the harbor, the water and sediments attempt to reach 
equilibrium in light of the new conditions, frequently much to his discomfort, 
as witnessed by the more than 100 million yards of material that must be 
dredged annually from the harbors in the United States to keep them open. 

The Committee on Tidal Hydraulics of the Corps of Engineers, under the 
leadership of Clarence F. Wicker and J. B. Tiffany, for several years has 
been studying the environmental processes that affect depth of water in har- 
bors. These factors have been well-expressed by Wicker and Rosenzweig 
(1950) in the first published report of this committee. The present sympo- 
sium, of which this paper is a part, summarizes the progress that has been 
made by the committee since 1950. This paper describes the geologic causes 
of harbors and the environmental factors that affect the character of estuarine 
deposits in the United States. It also contains a brief description of the geo- 
logic cause of typical American harbors. 


Geologic Factors Affecting Harbors 


The principal geologic factors affecting estuaries or harbors are; (1) the 
rate of change of sea level, (2) the deformation of the land, (3) the climate or 
rate of change of the climate, (4) the nature and distribution of the bed rock 
in the area, and (5) the inherited configuration of the land resulting from its 
previous geologic history. The soil, which is one of the principal sources of 
the sediments deposited in the harbor, is also a result of these environmental 
factors, as Jenny (1950) has shown so clearly. Soil depends primarily upon 
the climate, bed rock, and topography or slope of the land. Man, in his efforts 
to cultivate the land, upsets the processes which produced the previous equi- 
librium or near equilibrium. The chief difference caused by man is accelera- 
tion of erosion and greater supply of sediment to the harbor, which results in 
shoaling problems of concern to the community. Man also, in deepening and 
straightening channels and in building breakwaters and jetties, interferes with 
the natural circulation, thus resulting in the deposition and scour of sediment 
in unwanted places. 

The geologic factors that influence the deposition and character of sedi- 
ments are essentially the same factors as the hydraulic engineer recognizes. 
The geologist, however, thinks primarily in terms of the present state, and of 
_ the rate of change, of the environment. Four fundamental factors control en- 
vironment. These factors are climate, configuration of land, nature of bed- 
rock, and past geologic history. The sediment and soils in individual environ- 
ments strive to attain, or have already attained, equilibrium with their en- 
vironment. If geologic events, or man, cause the environment to change, the — 
old equilibrium is upset and the area tries to come into equilibrium with the © 
new conditions. Thus, in areas in which geologic conditions or climate are, — 
or have been recently, changing at a relatively rapid rate, the sediments and 
depth of water are likely to be different than they were formerly. | 

As the factors causing shoaling in harbors are discussed by other papers — 
in the present symposium, they will not be described at length here. The 
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principal factors are current distribution; quantity, size distribution, and 
mineral composition of sediment; amount and rate at which water is brought 
down by streams and in by tides; the distribution of Salinity; and the addition 
of organic matter in the form of sewage to the harbor. The type of clay min- 
erals is a very important factor affecting deposition, but it is a result of a 
combination of the basic environmental factors just mentioned. The type of 
clay brought down by the stream materially influences the sediments deposited 
in the harbor. The presence of water of different salinity or ionic concentra- 
tion in the harbor causes changes in the clay minerals (Grim, 1953). 

As clays have only just begun to be studied intensively by geologists, exact 
conclusions as to what happens or may happen have not yet been reached. 
Certainly ions are exchanged in the clay molecules, causing changes in the 
physical and chemical behavior of the sediments. New clay minerals also 
may be formed. Shepard and his associates (1953, 1954) have found that on 
the continental shelf off Southern Texas, chlorite and glauconite take the place 
of common clay minerals brought down by the rivers. Recent work also has 
shown that organic matter forms compounds with the clay minerals (Grim, 
1950, pp. 265-269). That is, the organic matter ties into the structure of the 
clay and changes significantly its properties. Such reactions between organic 
matter and clay, aided probably by the effect of bacteria, perhaps are one of 
the causes of the serious silting below Philadelphia. In this area more than 
15 million cubic yards of material must be removed annually. Harbor engi- 
heers would do well to give serious attention to the problem of clay mineralogy 
with respect to silting of harbors. 


Effect of Change of Sea Level 


The principal geologic cause of harbors in the United States is the recent 
world-wide rise of sea level (Marmer, 1948); (Russell, 1952); (Shepard, 1948). 
Geologists are generally agreed that, not so long ago geologically speaking, 
sea level was lowered at least 200 feet and probably as much as 300 feet, and 
perhaps more, all over the world. The cause of this lowering of sea level is 
ascribed to the removal of water from the ocean to form the great glaciers 
that covered the polar regions of the world during the past ice age. These 
ice caps began to melt at a rapid rate 10,000 to 20,000 years, or perhaps 
longer ago. This melting has resulted in a gradual rise of sea level. The 
rise depended upon the vagaries of climate, and has proceeded at a variable 
rate. According to Disney (1955), sea level along the coasts of the United 
States has risen about half a foot in the last 50 years. It is not at all certain, 
however, that the rate was so rapid before that time. The distribution of 
ancient cities around the world suggests that sea level may not have changed 
appreciably during the last 2000 years. Nevertheless, there is little doubt 
but that sometime before the time of Christ sea level had risen materially 
from its previous low stage. As a result, the sea invaded the land and drowned 
the valleys of streams that flowed out to the ocean during glacial time. As 
sea level rose, silt was deposited in the drowned valleys. In places, particu- 
larly on the Gulf Coast, silting has proceeded almost as rapidly as rising of 
the sea (Thompson, 1955); in other places, such as Chesapeake and Delaware 
bays, sedimentation has not kept pace with rising level, with the result that 
long estuaries have formed. 
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Deformation of Land 


In addition to rise of sea level, the land itself in places has been deformed 
in relatively recent geologic time. The coast line near Seattle, San Francisco, 
and Boston has been depressed and broken by earth-deforming movements. 
This deformation has proceeded much more rapidly than the addition of sedi- 
ment, with the result that deep bays now exist in these regions. 

In other areas, such as the Atlantic and Gulf coasts, the land has been 
slowly tilted or depressed seaward. This tilting of the land resulted in down- 
cutting of the valleys in such a way that the present streams flow at lower 
elevation with respect to the adjoining land than in earlier geologic time. In 
places the axis of tilt is landward from shore with the result that the land 
seaward from the axis of tilt has been depressed below sea level. As this 
tilting began a long time ago, areas formerly occupied by sea have since been 
filled with sediment. An example of such warping is the Mississippi embay- 
ment shown in Fig. 1. This embayment has since been blanketed with sedi- 
ments of Tertiary and Quaternary age. 

In the areas formerly occupied by ice, north of a line extending northwest- 
ward from New York City to the Great Lakes and beyond, the land has been 
elevated as much as 700 feet since the retreat of the glaciers. The southern- 
most limit of tilting and the extent of the ice are shown upon Fig. 1 (Leverett 
and Taylor, 1915); (Flint, 1948). This tilting has resulted in depth changes 
along the New England Coast. The uplift caused by this tilting tends to lessen 
the effect of the previous deformation of the coast in that area. 

The warping of the continent has taken place in steps. Long periods of 
quiet have interrupted the deformation of the land. During these periods of 
quiet the streams have moved back and forth across their valleys and carved 
more or less broad flood plains. When the land began to move upward again, 
the streams cut new trenches in the old flood plains. The end result isa 
series of terraces. Most of these terraces are above Sea level, but Chesa- 
peake and Potomac bays have a pronounced terrace whose surface is 10 to 15 
feet beneath the present water level. Either this terrace was cut during the 
ice age when Sea level was lower, or the land has been tilted to such an extent 
that the former terrace was depressed below sea level. 


Differential Erosion 


The hardness of the underlying bed rock also influences the processes of 
erosion and the resulting configuration of land. If the rocks are hard they do 
not erode; if they are soft they wear away readily. As a result, valleys form 
where the rocks are soft. Hills or headlands are produced where the rocks 
are hard. Consequently, harbors or estuaries generally are found where the 
rocks are soft. If the rocks and sediments that form the borders of the es- 
tuaries or harbors are poorly consolidated or cemented, such as on the Gulf 
Coast, they erode readily by waves. As a result the harbors increase in size. 
The banks of Galveston Bay are estimated to be eroding at the rate of three 
feet a year (U.S. Corps of Engineers, 1942, p. 61). The borders of Chesapeake 
and other bays similarly are being eroded by the wind and waves, with the 
result that these bays too are increasing in size. 

The problem then arises—what happens to the silt that is eroded from the 
banks? Is it transported farther out into the bay and deposited in shoals, or 
is it carried out to sea? According to Thompson (1953), many bays in the 
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Gulf Coast maintain a constant depth depending upon the wave size, but Shepard 
(1953) cites data to show that other bays are silting up. Thompson points out 
‘that the larger bays are deeper than the smaller bays because they have 
greater fetch for the generation of waves. Asa consequence, wave action 
extends to a greater depth and stirs up the bottom. The depth of these bays 
thus represents an equilibrium condition in which any excess sediment that 
normally would be deposited is agitated by wave action and kept in suspension 
until such time as it is transported out to the ocean. Shepard estimates that 
the Trinity River in Texas transports five times as much sediment through 
Galveston Bay as is deposited in the bay. That is, more than 80 per cent of 
the sediment by-passes the bay and eventually reaches the sea. 

The problem, however, is not so simple as just outlined. Much depends 
upon the size of the particles that are transported to the bay or are already 
on the bottom of the bay. The sediments in the Gulf-Coast bays are fine- 
grained clays which easily stay in suspension once they are agitated. The 
Sand and silt brought down by the Trinity River into Galveston Bay, as Shepard 
has shown, have shoaled the water at the entrance of the river into the bay by 
as much as 5 feet in the past 100 years. It is also conceivable that a river 
could bring down so much fine silt and clay that the waves could not transport 
away all that was brought down, even though the material was readily trans- 
portable by the water in the bay. 


Effect of Currents 


Currents also affect silting in harbors. If currents are directed seaward, 
More sediment can be removed than if they are in some other direction. 
Longshore currents may produce barrier bars, which cause bays to form 
behind them; as for example, Pamlico Sound on the Atlantic Coast and Gal- 
veston Bay on the Gulf Coast. These bars generally produce smooth arcuate 
islands or spits parallel to the coast, leaving marshland and estuaries behind 
them with more or less narrow outlets. Bars are more common on drowned 
than on rising coasts. The California Coast has few barrier bars compared 
with the Atlantic Coast. 


Glacial Action 


Continental glaciers also profoundly modify harbors. The ice continually 
moves outward from centers of accumulation. As the ice migrates it picks 
up rock and silt. This debris is deposited when the ice melts. At the front 
of the glaciers, determined by the position where the rate of melting equals 
the rate of advance of the ice, large masses of debris called terminal moraines 
are formed. The rock material that is continually brought along by the gla- 
ciers is dumped as the ice melts. These terminal moraines, when located 
along the coast, may form islands which cause good harbors. Boston Harbor 
is an example. rsp 

Also when the glacier eventually melts, all the debris which is contained 
within it is left behind upon the ground as a blanket. This morainal material 
varies in composition from huge blocks many feet in diameter to silt and clay. 
Boston Harbor and Puget Sound were covered with ice which dropped a great 
variety of material on the bottom of the water. In Boston Harbor great varia- 
ion, from silt and clay to gravel, is found throughout the harbor. Also the 
jepth of water varies, depending upon the amount of debris deposited. So 
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much material has been laid down in some places that shoals or isles were 
formed. Puget Sound is so deep that the glacial material deposited in it has 
caused few problems of navigation. 


Geologic Description of Harbors 


Boston Harbor 

Boston Harbor is of glacial and deformational origin (Currier, 1952; 
Shaler, 1893). The New England Coast has been broken and depressed be- 
neath the surface in past geologic time. Following this action, the area was 
covered by thick sheets of ice, which when they melted left behind a more or 
less thick mass of morainal material on the ground and off the coast. Most 
of the sediment now found in Boston Harbor consists of this morainal mate- 
rial. Sediment size ranges from clay to huge boulders. Much of the clay 
presumably represents recent deposition, but most of the harbor bottom is 
just what the glacier happened to leave. The depth of water is variable, owing 
to the vagaries of glacial deposition. 

The subsurface sediments, as indicated in Fig. 2, consist of sand, silt, and 
clay of variable thickness, generally 20 to 30 feet. In places these sediments 
rest upon crystalline bedrock, in other places they rest upon blue clay, as a 
rule 40 to 60 feet thick. This clay locally has been oxidized to form yellow 
clay. Interspersed in some of the deposits are peat and organic clay. 


New York Harbor 

New York Harbor is essentially a drowned river. The Hudson River ex- 
tends seaward into a deep submarine canyon that crosses the continental 
shelf in a trench (Veatch and Smith, 1939). The canyon can be traced to a 
depth of more than 10,000 feet below sea level. In the vicinity of New York 
City, the presence of this canyon represents a former submergence of the 
land. Following that period of lowering, a series of silts and sandy silts 
were deposited. Subsequently, during the last glacial period when sea level 
was again lowered, the river cut a valley in these sands and silts. Since that 
time this glacial valley has been filled with sand and silt brought down the 
river. The depth of the ancient gorge of the Hudson River at New York is at 
least 280 feet (Fig. 2), and drill holes up the river have gone more than 400 
feet without striking bedrock (Fluhr, 1953). Long Island is covered by the 
end moraine of a big ice sheet. Long Island Sound occupies the low ground 
behind the moraine. Manhattan is composed of crystalline rocks. 


Philadelphia Harbor 

The harbors south of New York are limited landward by what is called the 
Fall Line. A line extending more than 200 miles southward from New York 
through Trenton, Philadelphia, Baltimore and Washington separates the soft 
Tertiary and Cretaceous sediments of the Atlantic Coastal Plain from the 
crystalline and other rocks of the Piedmont area to the west. Falls in the 
river at the contact between the crystalline and sedimentary rocks terminate 
navigation. The trace of this “Fall Line” and the distribution of crystalline 
and sedimentary rocks along the Atlantic Coast is shown upon Figure 1. : 

Philadelphia is near the headwaters of Delaware Bay. This bay represents 
the drowned valley of the Delaware River. Depth to sedimentary bedrock in 
places is more than 200 feet, and depth to the crystalline rocks is several 
hundred feet more. The old valley clearly was formed when the land stood 
considerably higher than now. Since the time of formation of the valley, the 
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land has sunk and the valley filled with sediment. Subsequently, the valley 
Was again deepened, during the glacial lowering of sea level. It has since 
been largely filled with material brought down by the Delaware and Schuylkill 
rivers. The silting problem in the harbor is complicated by the addition of 
finely powdered coal from the anthracite mines up the Schuylkill River and by 
the presence of organic clay deposits, which cause a serious dredging prob- 
lem downstream from Philadelphia. 


Norfolk Harbor 

Norfolk Harbor is at the mouth of the Chesapeake and Potomac estuaries. 
These estuaries are drowned valleys which have not yet filled with sediment 
Since the glacial lowering of sea level. The sediments are mainly silt and 
clay. 


Savannah and Charleston Harbors 

The rivers on the southern Atlantic Coastal Plain have not been depressed 
below sea level as much as the northern rivers. The harbors at Savannah 
and Charleston appear more to be the result of post-glacial rise of sea level 
than of downwarping of the Atlantic Coast, though downwarping cannot be posi- 
tively excluded. The problems in these harbors are mainly those caused by 
man when he tries to change nature’s equilibrium by deepening the water. 
Before man came, a natural depth of 10 to 15 feet prevailed. When these 
depths were exceeded, nature tended to fill in the dredged areas in places 
where currents were weak. 

The sediments consist of silt, clay and silty sand. The river brings down 
much silt and clay; a considerable proportion of the sand is brought in from 
the ocean by tidal currents, as is attested by the presence of hornblende in 
the beach sands and in the sand in the harbor, compared with its scarcity or 
absence in the sediments brought down by the streams entering the harbor 
(Schultz, 1952). The distribution of hornblende indicates that much of the 
Sand is transported 5 miles, and some of it as much as 10 miles, upstream 
from the ocean. As shown by Figure 2, the subsurface deposits consist of 10 
to 60 feet of sand, gravel, and shells, presumably in part representing ma- 
terial filled in during the post-glacial rise in sea level. Below this are sev- 
eral hundred feet of blue clay, shale, sandstone, and other sediments of Ter- 
tiary and Cretaceous ages. Similar subsurface conditions are found at 
Savannah. 


New Orleans Harbor 

The harbor at New Orleans is on the Mississippi River. As the harbor 
lies one hundred miles upstream from the mouth of the river, it is hardly an 
estuarine harbor, though salt water is now wedging a long way up the river. 
Under normal conditions the water ranges in depth between 35 and 100 feet 
from New Orleans to the principal mouths at Southwest and South passes. 
The volume of water coming down the river and going out Southwest Pass is 
sufficient to keep the sediment in suspension to considerable depths. The 
depth of water in the other channels is less. The chief problem of maintain- 
ing the channel to specified depths is to maintain adequate flow in the river 
and to control the outlets by proper jetties (Holle, 1951). 

The geologic history of the delta has been well described by Russell (1952), 
Morgan (1952), Fish (1952), and Scruton (1955). The sediments on the delta 
are very fine sand, silt, and clay. According to Holle (1951, p. 124), the 
river transports some 500 million tons of sediment annually. This sediment 
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is two per cent very fine sand (62 to 125 microns); 48 per cent silt (4 to 62 
microns); and 50 per cent clay (less than 4 microns). As in other harbors, 
silting is a more serious problem in the small channels than in the large 
channels. At New Orleans the subsurface sediments consist of 100 or more 
feet of sand, silt, and clay, below which are Pleistocene deposits, in part of 
oxidized clay (Fig. 2). These clays have been downwarped in relatively re- 
cent geologic time. 


Galveston Harbor 

Galveston Bay is typical of the barrier bays along the Gulf Coast. It is an 
elliptical bay, 20 miles in maximum diameter, protected from the ocean by 
a barrier bar composed of fine sand. The Trinity River, one of the major 
Texas streams, flows into the northwest part of the bay. The average depth 
of the bay is 7 to 9 feet. According to Shepard (1953) the bay has shoaled 
about two feet in the past century, and as much as five feet on the delta of the 
Trinity River. Shepard also states that more than 80 per cent of the silt 
brought into the bay is transported out to sea, and less than 20 per cent is 
deposited on the bottom. The sediments of the bay consist of mud and sandy 
silt to a depth of 50 to 100 feet, below which is firm clay of Pleistocene age. 
The surface of the Pleistocene clay slopes seaward gently. It is encountered 
at greater depth at the south than at the north end of the bay. 


San Diego Harbor 

The harbors on the West Coast are different from those on the East Coast. 
The United States west of the 105th meridian, as shown in Fig. 1, is essen- 
tially mountainous; whereas much of the East and Gulf Coasts are bordered 
by a more or less gently sloping plain. The coastline in California, Oregon, 
and Washington in general is bold and rugged. Natural harbors are few and 
far between; most of them are found where drowned valleys exist. The sedi- 
ments tend to be coarser than those in the harbors on the Atlantic Coast, per- 
haps owing to the greater relief and more arid climate on the West Coast. 
This inference, however, is subject to confirmation by detailed investigation 
of the problem. 

San Diego Bay is essentially a drowned bay that has not filled completely 
since the glacial lowering of sea level. It is not evident as to whether the 
drowning is the result of post-glacial rise of sea level or whether it repre- 
sents deformation or depression of the land. The sediments are sand and 
Silty clay to a depth of at least 100 feet (Fig. 3). 


San Francisco Harbor 

San Francisco Bay has a complicated geologic setting. It has been one of 
the most active areas of deformation of the crust of the earth in the United 
States for millions of years. Its present form results from a period of cru- 
stal deformation prior to Mid-Pleistocene, during which time a series of 
hills and valleys comparable to those now seen in the City of San Francisco 
were formed (Louderback, 1953). At this time the Sacramento River flowed 
into the ocean through Golden Gate (Fig. 4). Subsequently, the area now oc- | 
cupied by San Francisco Bay was downwarped and the Berkeley Hills elevated 
During this elevation of the Berkeley Hills, the Sacramento River maintained 
itself through Carquinez Strait. The elevation of the hills and the downwarp- 
ing of the bay went on at intervals during which several formations of sand, 
sandy clay, and clay were deposited. At least five formations have been 
recognized (Louderback, 1953); (Trask and Rolston, 1951). The oldest three 
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[ these formations consist of alternating sequences of clay, sandy clay, and 
and. The upper two formations are sand and soft silty clay. They are as- 
pciated with valleys carved in the older and firmer sand and clay during late 
lacial and post-glacial time. 

Recent work by the writer indicates that the lower three formations, which 
re older than the valleys formed during the last period of lowering of sea 
evel, contain a series of beds in definite order, beginning with waxy, fine- 
rained clay at the base, followed in turn by beds of sandy clay, and fine to 
iedium-grained sand. Each bed in this sequence is 5 to 10 feet in thickness. 
his three-member sequence of beds is repeated at least five times, and has 
een traced in bore holes over a horizontal distance of more than five miles 
n the east side of the bay. Seemingly, these sequences are associated with 
limatic changes, which perhaps could represent various phases of glacial 
tages of the Pleistocene era. The waxy clays at the base of each sequence 
ive the impression of having been formed during relatively humid periods of 
limate. Sand grains in these waxy clays seem to have been dissolved and 
isintegrated during the process of weathering, whereas the sand in the over- 
ring sandy clays has not been so changed. 

Following the deposition of these tripartite sequences of strata, the sea 
as lowered and a series of broad valieys, as much as 1000 feet in width, 
<tending to a depth of more than 200 feet, were cut inthem. An interesting 
ature of this drainage system is its course south and west of Yerba Buena 
land and over the rock sill between Yerba Buena Island and San Francisco, 
1 which the present bridge stands (Fig. 4). This valley cutting of the late 
‘e age progressed in steps, as data from bore holes indicate the presence 
two terraces on the sides of the old valleys. The stream which drained 
ie south end of the bay during early Pleistocene time went around the east 
id north sides of Yerba Buena Island. It did not extend across the sill, be- 
use the sill at that time formed a hill at least 100 feet high. 

The deposits on the surface of the bay consist of mud, sand, and gravel. 
he distribution of types of sediment on the present bay bottom is shown on 
ig. 5. This figure is based upon data presented by U. S. Coast and Geodetic 
irvey charts. Silty clay covers most of the bottom in shallow water and in 
aces protected from currents. On the east side of the bay this clay extends 
_places as much as three miles from shore, in water less than ten feet deep. 
also fills the valleys carved in the older firm clays. In the areas of these 
cient valleys, the mud may be as much as 100 feet thick. On the San Fran- 
sco side of the bay, the clay likewise is thick in the sites of the old valleys. 
ypical bore-hole profiles are shown in Fig. 3. 

The bottom of the bay in the center of Golden Gate is 360 feet deep and 
msists of large blocks of rocks. On either side of the deepest part of the 
ite, the texture of the deposits becomes progressively finer as the depth of 
iter lessens. First is a zone of fine gravel, which grades into coarse sand, 
th to the eastward and to the westward. Outside the Golden Gate is a 
escent-shaped bar, composed of fine to medium sand. This bar rises to 
thin 20 to 30 feet of the surface. It obviously represents an equilibrium 
ndition between tidal currents and sediment supply. Fine sand of about 
0 microns diameter is found in the bottom of the channel northward through 
2 upper bay and through the series of bays as far as the outlets of the San 
aquin and Sacramento rivers. 
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Portland Harbor 

Portland Harbor is located at the junction of the Willamette and Columbia 
rivers, 100 miles above the mouth of the Columbia. A serious problem in 
connection with this harbor is the formation of bars in the Columbia River 
downstream from Portland. Most of the deposits brought down as bed load 
in the river are fine sand. The average diameter of this sand, according to 
Hodge (1934) is about 200 microns. Hodge states that sand drifts both to the 
north and the south at the mouth of the river, but the greatest drift is to the 
north. 


Seattle Harbor 

Seattle Harbor is located on Puget Sound. This sound represents a part 
of the earth’s crust that has been sharply depressed in relatively recent geo 
logic time. The water is mostly 500 to 800 feet deep, with :a maximum dept 
of 930 feet. The deposits consist mainly of soft, silty clay, with sand where 
sandy beaches are present. Puget Sound is connected with the ocean by Ad- 
miralty Inlet and the Strait of San Juan de Fuca. The water in Admiralty Inl 
is shallower than in Puget Sound, and in one place the inlet is crossed by a 
sill on which the water is only 240 feet deep. Puget Sound was covered by ic 
during the Pleistocene Age. Depth to bedrock varies greatly in the vicinity 
of Puget Sound, as shown by the borings illustrated in Fig. 3. The subsurfac 
deposits encountered in these bore holes consist of sand, silt, and gravel lef 
behind by the glacier, or by streams which reworked the glacial deposits. 
According to Fleming (1954), sea level was at least 250 feet lower in Puget 
Sound during the late glacial period than it is now. 


SUMMARY AND CONCLUSIONS 


The necessary requisites of harbors are protection from waves, sufficie1 
depth of water, and space for maneuvering ships. These requirements are 
met in several ways geologically. The most common is inundation of river 
mouths by rising sea level, as represented by harbors at New York, Phila- 
delphia, Norfolk, Charleston, and San Diego. Another process is subsidence 
and deformation of the land so that the sea occupies low areas. Examples 
are San Francisco, Seattle, and Boston. Another factor affecting harbor for 
mation is differential erosion of the land either before or after the sea cams 
in. Examples of differential erosion before the advent of the sea are Delaw: 
and Chesapeake bays. Galveston Bay is an example of erosion by wave ac- 
tion after inundation. Bars deposited across the mouths of rivers, as at 
Galveston, also may lead to the formation of harbors. Deposits left by gla- 
ciers may form islands and leave deposits that protect bays from the open 
Sea, aS at Boston. Large and deep rivers, such as the Mississippi and Co- 
lumbia, carry so much water that they keep their channels open to consider: 
able depth for a long distance, thus allowing the development of inland ports 
such as New Orleans and Portland, Oregon. The depth of water in harbors 
is controlled by the amount of rise of sea level or sinking of the land and by 
the deposition of silt in the harbor. Factors that affect silting are distribu- 
tion and strength of currents; quantity, size-distribution and mineral com- 
position of sediment; amount and rate at which water is brought down by 
streams; distribution of salinity; type of clay minerals; and es of orga 
matter and sewage in the water. 
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Table 1 
List of Figures 
Geologic features of the United States. Glacial data after Flint (1948); 


other data from geologic map of the United States issued by the U.S. 
Geological Survey. 


Logs of typical bore holes in harbors on the East Coast and the Gulf 
of Mexico. From U. S. Corps of Engineers (1952). 


Logs of typical bore holes on the West Coast of the United States. 
From U. S. Corps of Engineers (1952). 


Ancient valley patterns, San Francisco Bay, California. 


Distribution of sediments, San Francisco Bay, California. 


April, 19 


HY 2 


956-14 


S2 vim Z1v9s 


SU0IsIO man aN 


wouvosos 


54 


wossejsoy? 


LNIWAVEGWIZ 


FIO, ION 


“er 
} 
Preesepeliad (} . 


BZ 


“gy \ 


wossogh . 
u 


“43002 -O00S NO!I1VA313 


UOsSOAjOD 


06910 vos 


Sut ' 


096/IU019 UOS 


Puosss0g 


eB) 


——~__ 40085 
ae \ *@ ~ 


SES 


SCE 


SAVANNAH NEW ORLEANS 


CHARLESTON 


PHILADELPHIA 


W_YORK 


AUVs 


4S NOGBNOS-IwNV>D 
W4ildSOHW ALInWHD 
TwilidSOH VA 


NW330 1¥-N3A3U9S 1d 
OD B3MOd 3373 AVS 


OD W3IMOd 33913 AVS 


@Y¥OD Y3dVd 
ONY OVE NOINN 


433F0Nd ONISNOH 
390!1u8 ¥3IdOOD 


3001ue w3e002 | fj 


aN mas 
ENS 


373M 4330 


w33uD AVMOT IY 
wv3IN 113M 2S ON 


PN 'N3anv29 ON T714EE SRE SONNE 
CN NZaMWD £ ON [see els SRS SSS 


AVG WUVMIN 


| SYONY3AOD 
40 


3N_QA 0001 |. | 
ALID NOINN fr 


“wy NOSONH 


SMOUYYN 


290148 

BIA DASAW 
T3NNN4 8345 3HDBOO 
TINNNL 


O84 IW HABON | 


18OguIv NV9O7 
TINNNA JissvOL 


TRASK 


4333 NI ¥31VM MOT NW3W M0738 Hid30 


8 g 8 
Nu on 3 
z ae 
Ri sa3 S3zoee8 
| ¥ |se8tesa 
: > Szeesz3 
NQNEEEESN 
SSSSAY ASS 
ISNEENSEN a 
3 Peed 
- Fa) w uJ 
2 Bebleea) xb s 
< rale< ru 
sslssele2| 2 2 
.SONWS w 
[| svaaHs™ RSS aNo4sanin 
NR AAAS 3NOLS3WI 809 —> 
EEE NQQQQQ HAN BNOAS3MIT 119 ——e 
RSE NRN FAQS 3NOASaMIT 0101 —e 
= = > 
od es 
s 3$ 
F a* , 
tn : z 
PASSE z 8 
NJ 0 $ 4 


w 
E>] 
ay 
o 


907 ON 


FILL 


(é 370vHS) 
w20u F AVID 


W/SMELLS 


> 


HARD CLA’ 
HARDPAN 


LEGEND 


”20N 
NvdduvH 


& 
1334 NI 434VM MOT NV3N M0738 H1d30 


ai 
s 
' 


447 
ses 


4 


« 
[} 
z 
#4 
« 
ae See 3 
> we wl 
+o Ir 
oz £2 22 gn 
ZS 52 << << 
40 54 dS Sz 
” Vo nS 
iti 
i 


(62) 
1269 


269 


SAND 


a 
> 
< 
« 
oO 


o| GRAVEL 
SAND & 


e 

5 

= 
FEROS 


MUD OR 
MUCK 


LOGS OF BORINGS 
ATLANTIC AND GULF COAST HARBORS 


H 
695 


CONSISTENCY 
F-FIRM H-HARD M- MEDIUM 


SOFT 


ie 


956-15 


” 
iS 
i 
a 
fa} 
BS 
3 


144 


LIMESTONE 


25 JUNE 52 


2136 


956-16 HY 2 April, 195 


aad er feats MOT NV aN sy Ad Hieao 


SHALE, ROCK OR 
26 JUNE 52 


(He) HAROPAN 


Jeera eae AS 
ee EE SSN 


Hei a A ss: oo 
lvat: ie NE p 
ane RE PE: ia ERE 
REL as é . 
ove Mat ogee iebateed* g* e* 
SST A 


j 


4 


CONSISTENCY 
3-SOFT F-FIRM H-HARO M- MEDIUM 


LEGEND 
GRAVEL 
GRAVEL 


FIGURE 3 
LOGS OF BORINGS 


PACIFIC COAST HARBORS 


Sy 


a “Shue” 


ASTORIA 


3NOLSONVS 


COOS BAY 


iencene P— MEBEIRSESTS Sa) STVHS ¥ aNossonvs 
“Gnonmon PMB SE CONST ES SSSA awnes 1a05 


| eee PRIN EWR 
| AH po INNS SRS fowwnens ee 


a0 sonr bo SNNINGEEEEE SNE SE 
a SSS 
| mamas b WAVES Ee 


le suainon DP | 


° 
: 
< 
q 
uw 
z 
< 
7) 


SuvIUs AVIO/M NOMVRYOT OVORd NVS 
INOT SNAIVIE WILVM WNIGSVD QOVeIATIS yd 


ao BEN Sy nt Simba eae 3208 
79 $f NN ee i 
pee 


Lait dpe dh 
\ FE SEX RR QS RRR ES SFO SEY | 


LOS ANGELES 


é 

34 01 ON JEN g 
$9 0N BESSY 

5 2) 0N BRBY | 


SAN_DIEGO 


123d NIWAZLVM MOT NVZN MOT28 H1ids0 


956-17 


TRASK 


ASCE 


wSt ee! ia 4 q DE SE e2Zl Ov e221 


2 


SAO 10190)6 A\J03 ame ame om 
SAajOA (0190) 2407 qx 


VINUOSITVD ‘AVE OOSIONVYS NYS. 
ee SAZTIVA IWI0V19 31V7 


Ning 
idle) 
Ny 
& 
4 


puowyds 


® 3uNdis 


y 


April, 1956 


HY 2 


956-18 


puos A418 Sete meee 
@ puos aula E24 a 


puos Ayi1s 
- 8 4S Apuos = 


us @ hor ES 


fe VINYOSITVD ‘AVE ODSIDNVES nvs 
'LN Ine Pe) 
‘"s| SLNZNIGZS 
$ 3uNdI4 


i Paper 957 HY 2 


JOURNAL 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 
ee 
SSS 


USE OF ZONING PRINCIPLES IN FLOOD PLAIN REGULATION* 
Joseph I. Perrey, ! M. ASCE 


SYNOPSIS 


Regulation of the use of flood plains should provide for protection of the 
floodway required for carrying flood waters and for the prevention of flood 
damage through unwise use of those parts of the flood plain not required for 
the floodway. 

Zoning ordinances for the protection of the water-carrying capacity of 
floodways should prohibit all new construction or alterations of existing 
structures that may reduce the floodway cross-section or otherwise limit the 
Capacity of the floodway. The construction of bridges, utilities, or other 
necessary structures in floodways should be permitted only after careful 
Study by competent engineers has determined that the proposed structure will 
not unduly restrict the capacity of nor adversely affect the efficiency of the 
floodway. 

The objective of regulation of occupance of the flood plain is to reduce or 
prevent damage resulting from flooding. This objective may be reached by 
the use of zoning ordinances, building ordinances, and subdivision control. 

Good administration of zoning ordinances or other types of flood plain 
regulation requires that floodways and flood plains be clearly defined on maps 
So there can be no misunderstanding of the area subject to regulation. Ele- 
vations of flood heights or other controlling elevations should be stated in 
ordinances and such indefinite terms as “flood height” or “high water mark” 
Should be avoided. 

Quotations from several ordinances and laws, now in use in various cities, 
counties and states, furnish guidance for those who wish to draft ordinances 
for the regulation of flood plains in their own communities. 


INTRODUCTION 


The need for some regulation of the use of our flood plains to reduce the 
ever mounting toll of lives and damages taken by floods has been recognized 
Note: Discussion open until September 1, 1956. Paper 957 is part of the copyrighted 
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for many years. However, little has been done to achieve any substantial de- 
gree of control of flood plain development until quite recently. Even now the 
efforts in this field are minute and scattered when compared with the country’s 
overall flood control needs. 

The vast damage resulting from the recent New England floods generated 
by hurricane Diane has given added emphasis and publicity to the need for 
flood plain regulation, particularly on the smaller streams where the control 
of floods by reservoirs is difficult. 

The Engineering News-Record, on September 1, 1955, in editorializing on 
the effect of hurricane Diane on future flood control planning stated: 


“From a more optimistic viewpoint, - - there are things that can be done 
to minimize the destructiveness of floods in small river valleys if the af- 
fected communities will stand for the radical changes that are necessary 
and accept a reasonable share of the costs that are involved. For example, 
to provide a safe channel capacity, many buildings may have to be removed 
or set back, bridges raised and lengthened and the channel bottom cleaned 
of debris accumulated over many years. Possibly a parallel flood channel 
may be required, or flood walls or dikes may have to be built. 


“Outside the towns, pumping plants, power stations and factories that 
cannot be moved out of the flood plain can be diked. Inside or outside, 
dwellings may have to be prohibited in areas known to be subject to dan- 
gerous flooding. And finally, if sites for flood detention reservoirs are 
available, the necessary dams can be built.” 


The editorial points up two necessary actions of importance to planners. 
The redeveloping phase, in which the old mistakes of encroaching on the 
stream and its floodway are corrected, and the development phase, in which 
prevention of future damage is obtained by proper regulation of the use of the 
flood plain. 

The corrective action needed in New England is no different from that 
needed in hundreds of other areas. In all such areas there is a need for 
comprehensive planning for flood control which will take into account all 
phases of community development. 


Comprehensive Planning 


The authority for local officials and planning commissions to develop and 
adopt master plans for the orderly development of governmental units and 
their environs is granted by state laws. The object of such legislation is to 
improve the health, safety, convenience and welfare of the citizens and to 
provide for future development of communities in a manner that will promote 
res general welfare and provide for the efficient and economical use of public 
unds. ; 

The state planning laws contain provisions for master plans to include 
comprehensive surveys and studies of existing conditions and recommenda- 
tions setting forth plans for the development, redevelopment, improvement, — 
extension and revision of physical situations. Flood control and prevention — 
are usually as prominent in the state enabling acts as Sanitation, utilities, 
transportation, education and many other factors that are generally given 
prominent places in community plans. 
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Comprehensive planning for flood control at the local level should provide 
for coordination of all means of flood control that may be applicable to the 
local situation. Improved channels, local protection works, retention reser- 


_voirs, limited use of the flood plain, and structural adjustments to floods 


should all be considered and fitted into their proper place with respect to the 
master plan for community or area development. Planning commissions 
should make special studies of the flood problems with the help of competent 
consultants, and the resulting plan should coordinate both the development 
and redevelopment phases. 

The master plan, however, is only an overall framework of desirable future 
developments and needs to be implemented by more detailed plans and regu- 
lations for accomplishment. Zoning ordinances, building ordinances and sub- 
division control are effective instruments for putting the master plan into 
operation. 

Zoning is an extension of the police powers of the state to local officials, 
to be used as an integral part of planning to the end that public health, safety 
(from fire, flood and other danger), comfort, morals, convenience and general 
welfare may be promoted. Zoning ordinances must be reasonable and court 
decisions indicate that the factor of reasonableness, which the law requires, 
is strongly supported by a comprehensive, expertly-prepared master plan. 


Floodway Requirements 


In any plan of development the channel capacity requirements of the flood- 
way must be given full recognition. This is essential to any sound plan that 
will permit the greatest use of the resources of the flood plain with the least 
damage. 

At this point it might be helpful to future clarity to define the terms flood 
plain and floodway. 

The flood plain is the area adjoining a river or stream which has been or 
hereafter may be covered by flood water. 

A floodway includes the channel of a river or stream and those portions 
of the flood plains adjoining the channel which are reasonably required to 
efficiently carry and discharge the flood water or flood flow of the river or 


stream. 


The capacity requirements of the floodway are determined by nature and 
cannot be ignored without disastrous results. Encroachments on a floodway 
by any property owner creates a nuisance that is detrimental to other prop- 
erty owners and the general welfare of the community. 

Therefore, it becomes necessary to have two types of regulation of flood 
plains to prevent unwise land use. One serves to preserve and protect the 
floodway and the other to protect users of the flood plains. Both of these 
types apply in the development and redevelopment phases of planning. 


Community Development Planning 


The prevention of future flood damage resulting from occupance of the 
floodway or flood plain can best be accomplished through a master plan for 
community development which has given adequate consideration to flooding 
and to floodway requirements. Such a plan also will provide for adequate 


control of the floodway and flood plains. 


957-4 HY 2 April, 1956 
Control of the Floodway 


The objective of floodway control is to preserve adequate space for carry- 
ing flood waters without increasing flood heights and to limit land use which 
may unduly reduce the efficiency and capacity of the floodway. 

The control or regulation, sometimes known as channel capacity regula- 
tion, as obtained through laws or ordinances that restrict or prohibit en- 
croachments on streams by dumping waste along the banks, filling the parts 
of the flood plain required for the floodway, or construction of bridges or 
other structures which decrease the channel cross-sectional area or other- 
wise restrict the floodway. 

This type of regulation is frequently contained in state flood control laws 
which are worded in general terms. For example, the Indiana Flood Control 
Act (1945) states: 


“It shall be unlawful to erect, use or maintain any structure in or on any 
floodway as a permanent abode or place of residence, or to erect, make, 
use or maintain any structure, obstruction, deposit, or excavation in or on 
any floodway - - -, which will adversely affect the efficiency of or unduly 
restrict the capacity of the floodway.” 


The states of Iowa, New Jersey, New York, Pennsylvania and others have 
somewhat similar laws to protect floodways. 

Although state laws of the type quoted above are fairly effective in con- 
trolling the design and construction of bridges and location of large industrial 
plants, they are frequently inadequate in regulating the development of new 
subdivisions and controlling the filling of flood plains and dumping of waste 
along streams. The failure in this respect is largely administrative and is 
due to inadequate personnel and facilities for policing the state and to the lack 
of defined floodways throughout the state for the guidance of flood plain users. 

The problem of encroachment on floodway and stream channels is one re- 
quiring constant vigilance and control at the local level. State and Federal 
agencies cannot and should not be expected to reach into every corner of the 
nation to police the actions of individual citizens who through ignorance or 
intent engage in objectionable encroachments. 

At the local level, the city, county or regional development plan should 
provide a floodway zone or district which sets aside the area required for the 
floodway. All occupance of the floodway should be prohibited except for those 
uses which will not adversely affect its efficiency or restrict its capacity. 

Some uses of floodways are definitely beneficial and may be permitted. 
Parkways and some agricultural uses, for example, will increase the effi- 
ciency of the floodway by keeping it cleared of undergrowth and other ob 
structions. 

Necessary structures such as bridges, should be permitted only when de- 
signed and constructed in accordance with plans for the improvement of the 
efficiency and capacity of the floodway, or-when adjustments are to be made 
in the floodway to compensate for any adverse effects of the construction. 
The determination of compliance with plans for floodway improvement or the - 
adequacy of compensating adjustments should be made by competent hydraulic 
engineers. 

Few local master plans and zoning ordinances are adequate in defining the 
limits of floodways and in providing for proper regulation or control. How- 
ever, the zoning ordinance (1951) of Milwaukee County, Wisconsin, provides 


ASCE PERREY 957-5 


an example of an ordinance in which a separate floodway district is established 
and regulated. The ordinance creates Upland, Valley and Channel Districts 
along streams, with the “upland” districts being on land above flood heights, 
the “valley” districts on the flood plains and the “channel” districts in the 
floodways. The ordinance for the “channel” districts contains the following 
provision: 


F-3—Channel District 

In the F-3 Channel District no building or structure, retaining wall or 
revetment wall, except bridges or dams, and no fillings --- or any kind of 
dump shall be placed. The boundaries of the F-3 Channel District shall 
be such lines as shall be officially determined from time to time by ordi- 
nance of the County Board ---. 


The boundaries of a floodway district should be clearly and precisely de- 
fined on maps of sufficiently large scale that there will be no misunderstand- 
ing of their locations. The ordinance cannot be properly administered unless 
the location of the zone of prohibited structures is definitely known. 

The determination of the floodway requirements and boundaries is a tech- 
nical task that involves a study of past and probable future flood discharges, 
characteristics of the floodway, and possible modifications of the floodway. 
History of past floods is frequently insufficient to furnish a safe guide for 
floodway capacity requirements and future greater floods must be considered. 

The technical experience for making such studies will be lacking on most 
planning boards and on most city engineering staffs. However, most states 
now have water resources agencies with responsibilities for making studies 
of floods and other water problems. Help in defining floodways usually can 
be obtained from such agencies. For detailed planning of improved channels 
or other works to increase floodway capacity it may be necessary to employ 
consultants with experience in this field. 

In administering exceptions from the floodway zoning ordinance for neces- 
Sary structures it may be advantageous for planning boards of small commu- 
nities to require approval of state agencies as a prior condition to granting 
the exception, particularly if state laws prohibit encroachment on floodways 
and if adequate technical review of the project is given by the state agencies. 


Control of the Flood Plain 


The objective of control of the flood plain outside of the floodway should be 
the protection of the lives and property of the users. This protection may be 
obtained by preventing unwise use of the flood plain by limiting activities to 
those that will not be unduly damaged by flooding. The means of control are 
obtained through zoning ordinances, building ordinances and subdivision 
control. 

Zoning ordinances employing land use regulations are a common devise 
for controlling flood plains. Many cities and counties now have such ordi- 
nances, of which the zoning ordinance (1953) of Kokomo, Indiana is an example: 


FP-1—Flood Plain District. 

This district is created to protect the public health and to reduce the 
financial burdens imposed on the community, its governmental units and 
its individuals by frequent and periodic floods and the overflow of lands. 
The boundaries of this district have been determined from data in the files 
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of the Indiana Flood Control and Water Resources Commission and as such 
lands become protected from the risk of overflow they shall be removed 
from the FP-1 District and be reclassified into the appropriate Use Dis- 
trict in Section 13 of this Ordinance. 


1.. Uses Permitted: 

a. Any uSe permitted in an A-1 and R-1 zone as long as the use does 
not require the erection of a structure intended for residential purposes 
nor shall any existing building so located and used for residential purposes 
be enlarged or structurally altered. 

b. Public parks, playgrounds, and recreational areas so long as a 
structure intended for residential purposes is not erected. 


2. Uses Prohibited: 
a. All uses, other than those permitted above are prohibited. 


The permitted uses in the A-1 and R-1 districts mentioned in the Kokomo 
ordinance included general agricultural use, public parks, playgrounds, rec- 
reational areas, and buildings other than residences. Hence, the only pro- 
hibited buildings are residences. 

Ordinances regulating only residences are too limited in-scope and do not 
provide protection against the unwise use of flood plains by industries, utili- 
ties and others who may be subject to severe flood damage, and whose damage 
may have a greater general economic effect on a community than damages to 
residences. 

The zoning ordinance (1949) of Azusa, California is somewhat broader in 
its restrictions than the Kokomo ordinance and has the following provisions: 


Water Conservation Districts 

Statement of Policy. There are certain areas of the City which under 
present conditions are not suited for permanent occupancy or residence 
by persons for the reason that they are subject to the periodical spreading 
of water and other hazards. Therefore, for the public interest, health, 
comfort, convenience, preservation of the public welfare, the City Council 
does hereby create districts within which it shall be unlawful to erect, 
construct, alter, or maintain structures of any type. 

Use Prohibited. All buildings or structures, mining and/or removal of 


rock, sand, and gravel, or manufacturing any products from the same for 
commercial purposes. 


The body of the Azusa ordinance contains a clear statement of the purpose 
of the ordinance and relates it to public safety and welfare, which is neces- 
sary for all such ordinances. 


The Zoning Ordinance (1950) of Shelbyville, Indiana, is stated in somewhat 
shorter terms: 


Article VI—Conservancy Use 
Section 13. C-1 Conservancy District ; 
The following regulations shall apply in all C-1 districts. 


A. Uses Permitted. 

1. Customary agricultural operations provided that no odor—or dust- 
producing substance or use shall be permitted within two hundred (200) 
feet of a residence or business district boundary. 

2. Municipal recreation uses. 


ASCE PERREY 957-7 


3. No building or structure designed or intended for permanent use or 
occupancy other than fences will be permitted. 


It should be noticed that the Kokomo ordinance termed the district a flood 
plain district and the other two ordinances used the terms conservation or 
conservancy. The Shelbyville ordinance contains no reference to flooding in 
any part of the ordinance although it was written to prevent unwise use of 
flood plains. The avoidance of reference to floods or flooding either in the 
title or body of a flood plain zoning ordinance defeats the primary purposes 
of protecting against flood damage by making the purpose of the act vague or 
even unknown and encourages relaxation in enforcement. Such practices are 
sometimes excused on the theory that favorable public relations cannot be 
maintained if a community’s unfavorable or dangerous features are acknowl- 
edged. However, in a test, the court’s attitude toward a flood plain zoning 
ordinance will be more favorable if the purpose of the ordinance is clearly 
stated and if public safety and welfare are stressed as a justification. 

Building ordinances can sometimes be used as a means of reducing poten- 
tial flood damage in flood plains. Such ordinances may require that the first 
floor of all buildings be placed above flood heights, may prohibit basements, 
or may require walls to be constructed to withstand water pressure. 

Such ordinances may be a part of building code requirements or more 
commonly may be a part of flood plain zoning ordinances. 

The Milwaukee County zoning ordinance (1951) makes use of structural 
adjustments in the portion of the ordinance creating “valley districts” (flood 
plain districts). 


F-2—Valley District: 

In the F-2 Valley District no building or structure will be erected, and 
no existing building or structure shall be moved unless the ground upon 
which said building or structure is to be erected and ten (10) feet beyond 
the limits of said building or structure shall, prior to or at the time of 
construction, be raised to such level that the main floor of said building 
or structure shall be not less than three (3) feet above the high water level, 
as shown on the maps above referred to. No basement floor or other floor 
shall be constructed below or at lower elevation than the main floor. 


Prohibitions against building below the high water mark are to be found in 
zoning or building ordinances of other communities but frequently they are 
buried among other regulations and are not recognized as flood plain regula- 
tions. Also, the elevation of the high water mark to be used in each locality 
must be given in the ordinance or be shown on zoning maps, as in the Mil- 
waukee County ordinance, if the ordinance is to be effective. 

Subdivision control is a very effective way of regulating the development 
of flood plain areas because it controls the very beginning of urbanization. 
The control should begin with the approval of the plat plan and should con- 
tinue through the development and construction stages. No construction 
should be permitted to start until adequate flood protection is provided for 
in the plan and made a part of the overall development. 

The methods of control consist of restrictions on the use of flood plains 
for new subdivisions, requiring structural adjustments, or requiring the use 
of flood protection works. All of these methods have been used successfully. 

The restrictions on use method is employed in the Kokomo, Indiana ordi- 
nance (1953) for the control of the subdivision of land: 
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1. No land shall be subdivided for residential use unless adequate access 
to the land over improved streets or thoroughfares exists or will be pro- 
vided by the subdivider from said land to a thoroughfare maintained by the 
governmental unit; or if such land is considered by the Commission to be 
unsuitable for such use by reason of flooding or improper drainage, (under- 
scoring added), objectionable earth and rock formation, topography or any 
other feature harmful to the health and safety of possible residents and the 
community as a whole. 


The use of structural adjustments is employed by Hammond, Indiana in its 
subdivision control ordinance (1955) by establishing a minimum elevation for 
streets: 


Part VI0—Minimum Street Elevation. 

1. The minimum elevation of residential streets shall be 596 feet mean 
sea level. 

2. Any special adaptation of this elevation so as to allow new streets to 
meet the grades of existing streets which are below this minimum shall be 
submitted in plan to the Board for written approval. 


The control of residential buildings is obtained in another section of the 
Hammond ordinance which requires that “the site of any residential building 
shall be filled and/or graded so as to reach a minimum elevation of 1 foot 
above the top surface of the curb of the street facing the lot on which the site 
is located.” 

Flood protection works, such as levees, may be used to advantage in places 
and should be required by ordinance for urban areas if no other means of 
protection are provided. 

No subdivision should be permitted to develop in an area subject to flooding 
unless adequate protection against floods is provided by the subdivider. If 
this policy is followed the subdivision will not become an economic burden on 
the community at a later date when enraged residents, after a flood, demand 
that flood protection be provided. Further, it will tend to prevent or retard 
the development of areas subject to flooding when other areas can be developed 
more economically. 

The Hammond, Indiana zoning ordinance (1955) establishes Flood Plain 
Districts and permits construction of residences if protected by the owner. 
The Hammond ordinance reads in part: 


Chart II—Flood Plain District 


Uses Permitted 
(c) Single family residences under the following conditions: 

1. The owners of the land must provide and pay for a system of adequate 
flood protection works along the course of the Little Calumet river. This 
system of flood protection works --- must be submitted in plan to the Indi- 
ana Flood Control and Water Resources Commission and approved by that 
Commission --- before it is considered by City of Hammond authorities 
for action. 

2. After approval of such plans by the Flood Control and Water Re- 
sources Commission, they shall be examined by the City Engineer and 
approved by the Board of Public Works and Safety. 

4. The standards to be followed by the Board of Public Works and Safety 
for any dikes included in these works shall be as follows: Such dikes shall 
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be constructed of impervious material; minimum elevation 598 feet mean 
sea level; uppermost surface width ten (10) feet, graded level and paved 
with cement or slag capable of bearing vehicular traffic; slopes on sides 
1-1/2 to 1, the river side being protected by an 18 inch dump of one to six 
inch stone rip-rap, the outside slope being closely seeded in grass, 

5. The City of Hammond shall be given a permanent easement for such 
works and dikes and their rights-of-way for purposes of their maintenance 
and control. 


The provisions of this ordinance are intended to either lift the burden of 
providing flood protection from the city and place it on the subdivider, or 
discourage the development of lands subject to flooding. 

One subdivision at Speedway, a suburb of Indianapolis, Indiana, has been 
developed with a flood protection levee built by the subdivider. The impelling 
force that caused the subdivider to furnish the protection, in this case, was 
not the requirements of the city or county plan commissions nor flood plain 
ordinances, but the refusal of the Federal Housing Administration to guaran- 
tee loans on any of the houses to be built in the area until the flood protection 
works had been approved by the Indiana Flood Control and Water Resources 
Commission. 


Redevelopment Planning 


Unwise past use of floodways and flood plains has led to the growth of 
blighted areas in many communities. Such areas have a depressing effect on 
adjoining areas and are a deterent to wholesome community growth. The 
correction of such conditions calls for bold planning and a remolding of com- 
munity activity. The corrective measures frequently extend over a long 
period of years. 

The redevelopment should be a part of the community master plan and 
should include the removal from the flood plain of all activities that are seri- 
ously damaged by flooding, and the removal of all buildings and obstructions 
from the floodway. However, through the use of flood protection works and 
structural adjustments many desirable and wholesome uses may be provided 
for on the flood plain. 

The master plan of the City of Cincinnati furnishes an excellent example 
of a river front redevelopment plan. The redevelopment area adjoins the 
central business section and has a most favorable location in relation to trans- 
portation and business. The plan calls for the removal of present decadent 
buildings and replacing them with parkways, stadium, heliport, administration 
and service center, convention and recreation center, and apartment buildings. 

The stadium will be built to withstand flooding without serious damage, the 
heliport deck will be above flood levels and have parking space underneath. 
All buildings will be built on the higher ground in the area with first floors 
above flood heights. Floods will cause temporary disruptions of some uses, 
such as parking, but no serious damage will result. 


SUMMARY 


The floodway requirements of rivers and streams as well as the human 
occupancy of the flood plains must be considered if we are to obtain the 
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greatest use of the resources of our flood plains with a minimum of damage. 
It is of primary importance that adequate capacity be preserved in the flood- 
way and that no encroachments be permitted. Flood plains, on the other hand, 
may be occupied under limited conditions. 

Floodway ordinances should permit encroachment only for very limited 
purposes, such as bridges, and then only when structural adjustments are 
made to the conditions of the floodway and no adverse effects on the channel 
capacity of the floodway will result. Flood plain ordinances, on the other 
hand, may permit more general use but such use should be conditioned on the 
prevention of flood damage through the employment of non-damagable types 
of installations or activities, structural adjustments or protective works. 

The limits of floodways should be based on thorough hydraulic studies of 
each local area. These limits, which will be determined by selection of a 
design flood for flood plain zoning purposes, will have a decided economic 
effect on the users of the flood plain and should be given careful considera- 
tion. Such studies of the magnitude and frequency of floods and channel ca- 
pacity requirements are technical in nature and should be made by competent 
hydraulic engineers. 

Many zoning ordinances place special emphasis on restrictions on resi- 
dential use of flood plains. This possibly results from a mistaken idea that 
protection against danger to lives and damage to homes should be a principle 
consideration. This idea frequently has led to the illogical conclusion that 
areas unsuitable for residential purposes should be satisfactory for industry 
and should be so zoned. Protection for all community activities should be the 
primary consideration, including industry, transportation, utilities, and rec- 
reation. Flood damage to poorly located industries frequently is much greater 
and has a more far reaching effect on the economy of a community than does 
damage to residences. 

The control of subdivision development for flood protection purposes should 
be at the point of approval of the plan for the subdivision. If the area to be 
developed is subject to flooding the original plan should not be approved until 
protection against flooding is provided. Once work has been started on the 
construction of streets and houses the problem of control becomes increas- 
ingly difficult as time goes on. Developers and other users of land in a flood 
plain frequently exert strong pressure against regulations that limit their use 
of land and such pressure will be intensified if restrictions are not enforced 
until after a large investment has been made in a partly completed improve- 
ment. 

Good administration of zoning ordinances requires good tools and good in- 
Spection procedures. The tools are adequate zoning maps with floodway and 
flood plain districts clearly defined, pertinent elevations shown on maps or 
given in ordinances, and sufficient approved bench marks in each district for 
use in establishing correct elevations. Inspections should be made in the de- 
velopment area before, during and after construction to determine the degree 
of compliance with regulations and to obtain a basis for corrective action if 
required. 

If planning commissions and zoning boards will accept comprehensive flood — 
control planning implemented by flood plain zoning as being as necessary to 
the health and economic welfare of a community as planning for adequate 
streets, Sewers, schools, shopping centers, industrial districts and other 
community improvements, then much can be done to limit and reduce the in- 
convenience, discomfort and damage resulting from floods. 
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Task Forces (continued) 


Flood Committee 
Economic Effects of Floods on Agriculture. 
Effect of Floods on Industry. 
Effect of Floods on Transportation. 
Maintenance and Operation of Flood Control Works. 
Use and control of Flood Mains. 


Research Committee 
Open Channel Aerated Flow. 
Cavitation in Hydraulic Structures. 


Sedimentation Committee 
Manual on Sedimentation 
Rates of Reservoir Sedimentation 
Distribution of Sedimentation in Reservoirs. 


The Editor regrets that space limitations prevent listing the names of 
our members serving on the Task Forces. Please refer to the Official 
Register 1956 for their names. 


The Tidal Hydraulics Committee has not appointed any Task Forces as 
yet since they are the newest addition to the Hydraulics Division. Chairman 
Wicker writes “Its area of interest is not merely the rise and fall of the tide 
and the flow in tidal waterways, but also the transportation and deposition of 
sediments under this complicated regimen, the effects of salinity intrusion 
on the tidal flows, and many allied subjects. At present, it is planning to 
broaden the interest in the subject and to expand the fund of knowledge by 
sponsoring papers at forthcoming Society Conventions.” 


Each year one member retires from a committee. Under normal condi- 
tions they have served your Division for five years. The success of a Divi- 
sion is in direct proportion to the amount of time and energy these men 
devote to the Society. Our Division is among the most active group, and it 
is our aim to maintain that position. 

We will miss George Hickox, C. W. Thomas, G. R. Schneider, R. K. 
Linsley, F. B. Laverty, and A. S. Fry. There is a deep feeling of apprecia- 
tion for their contributions to our Division. 


Director L. E. Rydell has been appointed as Contact Member for the 
Hydraulic Division. 
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Glenn E. Hands, President, Kansas City Section reports that the Hydrau- 
lic Conference was a huge success, 347 attended with representatives from 
17 states and 4 foreign countries. The papers were well received by the 
audience, reflecting careful preparation and presentation. 


Harvey O. Banks, Acting State Engineer of California, is Chairman of 
the Task Force on Spillway Design Floods. 

D. W. Van Tuyl is Chairman of the Task Force on Hydrologic Data. 

William F. Guyton, Chairman, Hydrology Committee, plans to hold a 
meeting of the committee at an early date. He reports the two Task Forces 
are getting their work underway and it appears they will have a successful 
year. 


The Committee on Snow, Ice and Permafrost originally a Joint committee 
was assigned to the Hydraulics Division. During the past two years your 
Executive Committee has questioned the advisability of sponsoring this 
group. Action was taken at the Berkeley meeting; the society committee on 
Technical Activities has been requested to transfer this group to the Soil 
Mechanics Division. 

A. F. Ghiglione, Commissioner of Roads for Alaska, Chairman of this 
Committee, concurs with the action of your Executive Committee. 


PLAN NOW TO ATTEND THE DIVISION MEETING 


Location. 

University of Wisconsin at Madison, Wisconsin. 
Date. 

August 22, 23, and 24, 1956. 
Technical Sessions. 

Design 

Floods 

Hydrology 

Research 

Sedimentation 
Host Section. 

Wisconsin. 


Professor Arno Lenz, President of the Section and a member of our 
Division, is directing the planning. 

Housing will be available at the dormitories for a package rate of $16.00 
per person for Tuesday, Wednesday, and Thursday nights and 3 meals on 


Wednesday and Thursday, and 1 meal on Friday. 
Tentative plans include an inspection trip to the Allis-Chalmers plant. 


A trip to a large brewery and also a night at the ball game. 
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Watch for further announcements in Civil Engineering. Details of ac- 
commodations at Hotels and Motels will be printed with the program. 


Hydraulic Division Sessions at Society Conventions 


To insure participation at all meetings, the Vice Chairman will be res- 
ponsible for the programs for the next calendar year. By inaugurating this 
procedure there is a better distribution of the work within the Executive 
Committee. 


The Technical Committees plan the programs often with papers of local 
interest as well as papers on subjects of general interest not only to Hy- 
draulic Engineers, but to the entire profession. With three Society, one 
Division, and a number of Regional meetings each year these committees 
give generously of their time to insure the success of the meeting. It is 
necessary to plan ahead almost two years in order to give the author suf- 
ficient time to prepare his paper. The co-ordination of the program is dif- 
ficult with this long range planning, since many engineers as jobs and 
may even move to other countries. 


It is believed that some of our members have or could prepare papers 
that would qualify for publication if they understood the basic requirements. 
Presentation of a paper at a meeting is not a factor in determining its 
merits. Many contributions to engineering would not hold the attention of 
an audience because of the technical aspects, but would be of great value 
and should be published in the Transactions. 

Papers of local interest will draw an excellent audience, but may fail to 
qualify for publication. In spite of the fact that authors obviously spend a 
great amount of time preparing their papers, we regret that it is not pos- 
sible to publish all papers. 

The following was prepared by the Chairman of the Publications Commit- 
tee to assist authors in qualifying for a review of their paper by a board 
when submitted for publication. 


Requirements Before Papers are Reviewed by 


Division Committee on Publications 


1. Titles should not exceed a 50-character-and space length. 


2. Type manuscript, original and two copies, double-spaced on one side 
of 83" by 11" paper. Illustrations, tables, etc., must accompany each copy 
of the manuscript. Papers exceeding 18 000 andy in length require special 
approval from the Society. 
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3. Forward the original and two copies to Mr. Harold T. Larsen, Man- 
ager, Technical Publications, ASCE, 33 West 39th Street, New York 18, New 
York. He will forward two copies to Mr. H. G. Dewey, Jr., Chairman, Hy- 
draulics Division, Committee on Publications, c/o San Francisco District, 
Corps of Engineers, U. S. Army, 180 New Montgomery Street, San Francis- 
co, California. 


4. No line of mathematics should be longer than 61". Equations in the 
manuscript should be lettered in black India ink, with capital letters 3/16" 
high and all other symbols and characters in proportions dictated by stan- 
dard drafting practice. Original typed equations are acceptable, but they 
will be reduced to 69% of their original size. 


5. Tables must be typed (originals, not copies) on one side of 81" by 11" 
paper within an invisible frame 61" by 104". If tables are small, more than 
one should be included in the frame. Captions should be included within the 
frame. Reductions will be 69% of original size. 


6. Illustrations should be drawn with black India ink on one side of 81" 
by 11" paper within an invisible frame64"by 101". If some illustrations are 
small, more than one should be included in the frame. Captions must also 
be included within the frame. Glossy prints of photographs may be furnish- 
ed when manuscript is first submitted. Reductions will be 69% of original 
size. 


7. Submit with manuscript on a separate sheet, a 50-word abstract, 
typed double-spaced on one side of 84" by 11" paper. 


8. Papers prepared for delivery on the program at technical meetings 
senerally do not meet the above mechanical requirements. Accordingly, 
before submitting such papers to Mr. Larsen, appropriate revisions should 
be made. It is to be noted that delivering a paper at technical sessions does 
not imply such paper will necessarily be accepted later for publication in 
Proceedings. Each paper will be evaluated on its own merits by the Com- 


mittee on Publications. 


Presentation of Technical Papers 


Chairman Tiffany sent the Editor a copy of a paper entitled “Toward 
Better Presentation of Technical Papers” published by the American 
Society for Testing Materials. With their permission we will quote from 
this article. 

When you present a technical paper, do you interest and enlighten your 
4udience, or do you irritate and stupify it? Unfortunately the Society ex- 
yeriences too many presentations of the latter type.**** After all the 
fforts by the author, the reviewers, the Committee, and the Editorial Staff, 
in unintelligible oral presentation is a grevious disappointment but one 
which is experienced only too often. 

The relative merit of manuscript papers is not reflected in the way they 
ire presented at the technical sessions of the Society. The best paper can 
ye made a worthless mishmash by an author who fumbles his way through 
ts presentation. On the other hand, the author of a paper which barely has 
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made its way through the reviewing screen can steal the show by a combina-— 
tion of alertness, stage presence, and detailed preparation. *** 

An author can do a great deal to improve his presentation without as- 
sistance. He can make a thoughtful selection of the major points in the 
manuscript and then arrange them in logical, outline order for speaking. 
These should be coordinated with visual aids, such as slides and charts.**** 

The author should exercise special care in the preparation of slides or 
other visual aids, and should make certain that such material can be seen 
easily from all parts of the audience.**** 

The Papers Committee is convinced that anyone who has the ability to 
prepare an acceptable technical paper also can present it effectively if only 
he takes the time to prepare his oral presentation. 


Dead Lines For Papers 


Title and Author ---6 months before meeting 
Submit Paper ---4 months before meeting 


This lead time is necessary in order to make local arrangements, have the 
program printed in Civil Engineering and have the paper reviewed by the 
Technical Committee arranging the session. 


Papers of interest to Hydraulic Engineers often appear in the Journals 
of other Societies. 
“Erosion Resistance of Concrete in Hydraulic Structures.” 
Journal of American Concrete Institute Nov. 1955 
Forum on Hydraulics includes Water Hammer and Cavitation 
Mechanical Engineering Dec. 1955 


Membership in Two Divisions 


In the January issue of Civil Engineering page 125 you will find the form 
for registration in a second division of your choice. It appears desirable 
to call this to the attention of your friends and associates if they are inter- 
ested in joining two divisions. 

The Board of Direction is to be congratulated in making this possible as 
many of our members have more than one interest in the field of Civil 
Engineering. 


Editing of the News Letter 


In past years the Secretary of our Division was also the editor of our 
News Letter. Tospreadthe work among the Executive Committee members, 
it was decided that the Past Chairman would take this assignment. Since 
this bulletin’s purpose is to keep you informed about our activities the Editor 
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ild appreciate hearing from you as to what you would like included in 
his paper. 

If you have any suggestions please write to the Editor. If you would like 
} become more active in Division work let us know about it. We can not 
romise to put everyone on a committee or task force on short notice, but 
st of members that are willing to work would be of great value when 
lanning for the future. 


: = T. J. Corwin, Jr., Editor 


PROCEEDINGS PAPERS 


the technical papers published in the past year are identified by number below. Technical- 
sion sponsorship is indicated by an abbreviation at the end of each Paper Number, the 
bels referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering 
Hanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
tary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and 
ping (SU), and Waterways and Harbors (WW) divisions. Papers sponsored by the Board of 
setion are identified by the symbols (BD). For titles and order coupons, refer to the appro- 
te issue of “Civil Engineering." Beginning with Volume 82 (January 1956) papers were 
ished in Journals of the various Technical Divisions. To locate papers in the Journals, the 
bols after the paper numbers are followed bya numeral designating the issue of a particular 
“nal in which the paper appeared. For example, Paper 861 is identified as 861 (SM1) which 
cates that the paper is contained in issue 1 of the Journal of the Soil Mechanics and Founda- 
3 Division. 
VOLUME 81 (1955) 


UL: 659(ST), 660(ST), 661(ST)°, 662(ST), 663(ST), 664(ST)°, 665(HY)°, 666(HY), 667(HY), 
tly 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 
78 s 


= 679(ST), 680(ST), 681(ST), 682(ST)°, 683(ST), 684(ST), 685(SA), 686(SA), 687(SA), 688(SA), 
89(SA)°, 690(EM), 691(EM), 692(EM), 693(EM), 694(EM), 695(EM), 696(PO), 697(PO), 698(SA), 
99(PO)©, 700(PO), 701(ST)°. 


4; 702(HW), 703(HW), 704(HWw)°, 705(IR), 706(IR), 707(IR), 708(IR), 709(HY)°, 710(CP), 
"UA(CP), 712(CP), 713(CP)°, 714(HY), 715(HY), 716(HY), 717(HY), 718(SM)°, 719(HY)S, 
2O(AT), 721(AT), 722(SU), 723(WW), 724(WW), 725(WW), 726(wW)°, 727(WW), 728(IR), 
29(IR), 730(SU)*, 731(SU). 


x: 732(ST), 733(ST), 734(ST), 735(ST), 736(ST), 737(PO), 738(PO), -739(PO), 740(PO); 
’41(PO), 742(PO), 743(HY), 744(HY), 745(HY), 746(HY), 747(HY), 748(HY)°, 749(SA), 750(SA), 
'SL(SA), |752(SA)®, 753(SM), 754(SM), 755(SM), 756(SM), 757(SM), 758(CO)°, 759(SM)°, 
reo(ww)°, 


: 761(BD), 762(ST), 763(ST), 764(ST), 765(ST)°, 766(CP), 767(CP), 768(CP), 769(CP), 
NCP), 771(EM), 772(EM), 773(SA), 774(EM), 775(EM), 776(EM)°, 777(AT), 778(AT), 
(SA), 780(SA), 781(SA), 782(SA)°, 783(HW), 784(HW), 785(CP), 786(ST). 


TEMBER: 787(PO), 788(IR), 789(HY), 790(HY), 791(HY), 792(HY), 793(HY), 794(HY)°, 
95(EM), 796(EM), 797(EM), 798(EM), 799(EM)®, 800(WW), 801(WW), 802(WW), 803(WW), 
eww), 805(WW), 806(HY), 807(PO)°, 808(IR)°. 


OBER: 809 (ST), 810(HW)°, 811 (ST), 812(ST)°, 813 (ST)°, 814(EM), 815(EM), 816(EM), 
\17(EM), 818(EM), 819(EM)°, 820(SA), 821(SA), 822(SA)°, 823(HW), 824(HW). 


$MBE 825(ST), 826(HY), 827(ST), 828(ST), 829(ST), 830(sT), 831(ST)°, 832(CP), 
3(CP), ateeh fir 836(HY), 837(HY), 838(HY), 839(HY), 840(HY), 841(HY)°. 


BER: 84 843(SM)°, 844(SU), 845(SU)°, 846(SA), 847(SA), 848(SA)°, 849(ST)°, 
s ees, cen, 853(ST), 854(CO), 855(CO), 856(CO)°, 857(SU), 858(BD), 859(BD), 


oo 
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sie , 862(SM1), 863(EM1), 864(SM1), 865(SM1), 866(SM1), 867(SM1), 868(HW1) 
: shia ay pees 872(HW1), 873(HW1), 874(HW1), 875(HW1), 876(EM1)°, 877 
87 se . 


* 


JAF ; Y1), 881(HY1)°, 882(HY1), 883(HY1), 884(IR}), 885(SAJ), 886(CP1}, 
a1(S/ : ie bot ae doclaats 891(SA1), 892(SA1), 893(CP1), 894(CP1), 895(PO1) 
96(PO1), 897(PO1), 898(PO1), 899(PO1), 900(PO1), 901(PO1), 902(AT1)°, 903(IR1)°, 904 


ot + Y a 


928(8A2), 929(SA2), 930(SA2), 931(8A2), 932(SAz)°, 
ww), 937(wWw2), 938(Ww2), 939(WWw2), ), 
(ema) 248 1 DeEMa}e ), 


je! iy . 
oh Sig See ¢ 
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